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SUMMARY 


This  report  documents  version  8  of  the  OMEGA  10  program  and  version  3  of  the 
OMEGA  1 1  program.  These  programs  were  developed  to  prepare  flight  and  ground  run-up  noise 
data  for  input  to  NOISEMAP,  the  Air  Force  community  noise  exposure  prediction  program.  All 
routines  in  each  program  are  documented  at  a  level  useful  to  a  programmer  working  with  the 
code  or  a  reader  interested  in  a  general  overview  of  what  happens  within  a  specific  subroutine. 

Both  the  OMEGA  10  and  the  OMEGA  1 1  programs  input  reference  aircraft  noise  data 
from  the  NOISEFILE  database.  This  NOISEFILE  database  contains  reference  datasets  for  both 
flyover  and  ground  run-up  operations  for  almost  all  current  military  aircraft.  There  are  typically 
three  to  six  flyover  and  ground  run-up  datasets  for  each  aircraft  where  each  dataset  defines  a 
different  engine  power  setting.  The  flyover  datasets  contain  seven  mean  single  event  measures 
plus  the  mean  sound  pressure  level  spectrum  of  peak  perceived  noise  level  (PNLM),  all 
normalized  to  a  minimum  slant  range  of  1000  feet,  a  surface  temperature  of  59°F  and  a  surface 
relative  humidity  of  70%.  Each  flyover  dataset  is  also  normalized  to  a  specific  reference 
airspeed  which  varies  depending  on  the  type  of  aircraft  and  the  power  condition.  The  ground 
run-up  datasets  contain  19  farfield  sound  pressure  level  spectra  measured  at  10  degree  increments 
at  a  fixed  radial  distance  around  one  side  of  the  aircraft.  These  ground  run-up  data  are 
normalized  to  a  fixed  radial  distance  of  250  feet  and  to  standard  day  meteorological  conditions 
(59°F,  29.92  inches  Hg,  70%  relative  humidity).  The  format  and  content  of  these  normalized 
reference  flight  and  ground  run-up  data  are  defined  in  Appendices  E  and  G,  respectively. 

The  OMEGA  10  program  inputs  reference  flyover  datasets  from  the  NOISEFILE 
database  for  a  specific  aircraft,  and  extrapolates  the  reference  sound  pressure  level  (SPL)  data 
from  the  reference  minimum  slant  range  (1000  feet)  to  22  profile  distances  (200  to  25000  feet), 
computes  the  single  event  measures  at  these  distances,  and  then  applies  a  set  of  rules  to 
extrapolate  or  interpolate  these  single  event  versus  distance  data  to  produce  distance  profiles  for 
up  to  seven  single  event  noise  measures  at  the  requested  power  setting,  airspeed,  temperature  and 
relative  humidity.-  The  seven  single  event  measures  are  A-weighted  overall  sound  level  (ALX), 
tone-corrected  A-weighted  overall  sound  level  (ALTX),  perceived  noise  level  (PNLX),  tone- 
corrected  perceived  noise  level  (PNLTX),  sound  exposure  level  (SELX),  tone-corrected  sound 
exposure  level  (SELTX)  and  effective  perceived  noise  level  (EPNLX).  In  the  print  mode,  the 
profile  data  for  all  seven  measures  are  always  computed  and  printed  and,  when  requested  by  the 
IPU  flag,  the  SELX,  ALX  and  EPNLX  data  are  written  to  the  card  image  file.  In  the  no-print 
mode,  which  is  designed  primarily  to  prepare  data  for  input  to  the  NOISEMAP  program,  only  the 
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one  or  more  SELX,  ALX,  or  EPNLX  measures  identified  on  the  code  sheet  are  written  to  the 
card  image  file.  The  content  and  format  of  these  flight  noise  profile  datasets  in  the  card  image 
file  are  given  in  Appendix  F. 

The  set  of  rules  referenced  in  the  above  paragraph  describe  the  procedure  required  to 
interpolate  and/or  extrapolate  the  OMEGA  10  profile  data  from  the  given  reference  data.  These 
rules  are  defined  in  the  documentation  for  subroutine  SETUPD6. 

The  OMEGA  1 1  program  inputs  reference  ground  run-up  datasets  from  the  NOISEFILE 
database  for  a  specific  aircraft,  extrapolates  these  SPL  spectra  from  the  reference  distance  (250 
feet)  to  each  of  the  22  profile  distances,  computes  the  AL,  ALT,  PNL  and  PNLT  single  event 
measures  for  each  spectrum  at  each  distance,  and  then  interpolates  these  reference  data  to 
generate  similar  distance  profiles  for  ALX,  ALTX,  PNLX,  and  PNLTX  at  the  requested 
temperature,  relative  humidity,  barometric  pressure  and  aircraft  engine  power  setting.  As  in  the 
OMEGA  10  program  above,  print  flags  are  defined  to  control  the  type  and  quantity  of  data 
printed  on  the  output  and  card  image  files.  The  format  and  content  of  these  ground  run-up  noise 
profile  data  are  described  in  Appendix  H. 

The  profile  datasets  written  in  the  card  image  files  as  described  in  the  above  programs  are 
the  primary  aircraft  noise  data  used  by  the  NOISEMAP  program  to  generate  noise  profiles 
around  a  specific  airbase.  With  these  programs  these  noise  data  can  be  tailored  to  the  specific 
weather  and  aircraft  flight  and  ground  run-up  power  conditions  for  that  airbase  which  will,  in 
turn,  improve  the  accuracy  of  the  noise  contours  produced  by  the  NOISEMAP  program. 

The  remaining  appendices  in  this  report  contain  the  following  items  for  both  the  OMEGA 
10  and  OMEGA  1 1  programs: 

(a)  a  complete  source  listing; 

(b)  the  program  code  sheet  and  the  procedure  for  setup  and  execution  of  the  program; 

(c)  a  complete  sample  run  including  a  listing  of  all  input  and  output  data. 

Appendix  K  is  a  user's  guide  for  the  MENU  108  program  which  can  be  used  to 
interactively  setup  and  execute  the  OMEGA  10.8  program.  This  MENU108  program  also 
enables  the  operator  to  select  specific  noise  measure  profiles  which  are  then  written  to  the  display 
for  easy  comparison. 


iv 


A  similar  MENUl  13  program  is  being  developed  to  setup  and  execute  the  OMEGA  1 1.3 
program.  MENUl  13  will  display  noise  measure  profiles  for  22  distances  for  selected  angles  as 
well  as  noise  measure  data  for  all  available  angles  for  selected  distances.  This  MENUl  13 
program  is  not  complete  at  this  time  and  thus  is  not  included  in  this  report.  When  this 
MENUl  13  program  is  available,  it  will  present  menus  which  are  very  similar  to  the  MENU  108 
menus  described  in  Appendix  K. 
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PREFACE 


This  updated  report  documents  the  OMEGA  10.8  and  OMEGA  1 1.3  computer  programs 
which  were  developed  to  prepare  aircraft  flyover  and  ground  run-up  data  for  input  to  the 
NOISEMAP  program.  This  work  was  performed  for  the  Noise  Effects  Branch  of  the  Armstrong 
Laboratories  at  Wright-Patterson  Air  Force  Base,  Ohio.  The  contract  monitor  for  this  effort  was 
Mr.  Jerry  Speakman. 

Special  thanks  are  due  to  Mr.  Jerry  Speakman  and  Mr.  Robert  A.  Lee  for  their  guidance 
and  assistance  in  this  effort. 
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INTRODUCTION 


This  report  contains  the  complete  documentation  of  the  OMEGA  10.8  and  OMEGA 
1 1.3  programs.  Included  in  Appendices  are  the  program  source  listings,  input  parameter 
code  sheets  with  standard  procedures  for  setup  and  execution  of  each  program,  format 
descriptions  for  all  input  and  output  records,  and  input  and  output  sample  test  data.  Both 
programs  were  originally  written  in  CDC  FORTRAN  Extended  (FORTRAN  IV)  for  the 
CDC  CYBER  74  and  CYBER  750  in  the  ASD  Computer  Center  at  Wright-Patterson  Air 
Force  Base.  Both  programs  have  been  modified  to  run  on  an  IBM  compatible  PC  with  a 
Lahey  FORTRAN  compiler.  With  several  minor  changes,  both  programs  can  also  be 
compiled  with  Microsoft  FORTRAN  Version  5.0. 

These  programs  were  developed  by  the  University  of  Dayton  Research  Institute 
(UDRI)  and  the  Air  Force  Aerospace  Medical  Research  Laboratory  (AFAMRL)  to  prepare 
flight  (OMEGA  10)  and  ground  run-up  (OMEGA  11)  noise  data  for  input  to  NOISEMAP, 
the  Air  Force  community  noise  exposure  prediction  program.  Both  programs  input 
reference  datasets  from  the  NOISEFELE  database  and  interpolate  or  extrapolate  the  profile 
measure  data  at  the  requested  power  settings  and  meteorological  conditions.  These  flight 
and  ground  run-up  reference  data  in  NOISEFTLE  were  computed  by  the  OMEGA  6  and  8 
programs,  respectively.  Some  of  the  analysis  procedures  used  in  these  programs  are 
described  in  AMRL-TR-73-107.® 

The  OMEGA  10  program  inputs  all  flight  noise  reference  datasets  from  NOISEFILE 
for  the  aircraft  being  analyzed.  This  flight  noise  NOISEFILE  database  is  opened  with  Unit 
7  in  the  program.  These  flight  reference  datasets  contain  seven  mean  single  event  measures 
plus  the  mean  sound  pressure  level  spectrum  of  perceived  noise  level  (PNLM),  all 
normalized  to  1000  feet  minimum  slant  range  at  standard  day  meteorological  conditions 
(59°F  and  70%  relative  humidity).  A  description  of  the  format  and  content  of  these 
reference  datasets  is  given  in  Appendix  E.  These  reference  data  are  used  to  compute  the 
ALX,  ALTX,  PNLX,  PNLTX,  EPNLX,  SELX,  and/or  SELTX  profile  data  at  the  reference 
power  setting  for  the  profile  output  meteorological  conditions.  These  reference  power 
profile  data  are  then  used  to  interpolate  and/or  extrapolate  the  profile  data  to  the  requested 
power  setting  and  airspeed.  These  final  profile  data  are  printed  in  tabular  form  on  the 
OUTPUT  file  and/or  written  to  the  file  opened  on  Unit  3  in  the  format  required  by  the 
NOISEMAP  program.  Since  only  the  ALX,  SELX,  or  EPNLX  data  are  used  by  the 
NOISEMAP  program,  only  these  profile  datasets  can  be  written  to  the  file  opened  on  Unit  3. 


1 


A  description  of  the  format  and  content  of  this  Unit  3  file  is  given  in  Appendix  F.  A 
complete  sample  problem  with  a  listing  of  all  input  and  output  data  is  given  in  Appendix  C. 

The  OMEGA  1 1  program  inputs  all  ground  run-up  noise  reference  datasets  for  the 
aircraft  being  analyzed  from  file  NOISEFILE  which  is  opened  on  Unit  7.  These  reference 
datasets  contain  spectral  data  at  10°  angular  increments  for  angles  0  to  180  degrees  at  a 
specific  power  setting  and  standard  reference  distance  (250  feet)  and  for  standard  day 
meteorological  conditions  (59°F,  70%  relative  humidity,  and  29.92  in.  Hg.).  The  format  and 
content  of  these  reference  datasets  are  described  in  Appendix  G.  These  reference  SPL 
spectra  are  then  used  to  extrapolate  spectral  data  for  the  same  angles  at  22  standard  profile 
distances  for  any  reasonable  meteorological  conditions,  all  at  the  reference  power  settings. 
These  extrapolated  spectral  data  are  used  to  derive  PNLX,  PNLTX,  ALX  and  ALTX  data  for 
each  angle  and  profile  distance.  These  PNLTX,  ALX  and  ALTX  data  for  two  power 
settings  are  then  used  to  interpolate  the  final  ground  run-up  noise  profile  data  required  by  the 
NOISEMAP  program;  they  are  written  to  a  file  on  Unit  2  in  the  format  described  in 
Appendix  H.  A  complete  sample  problem  with  a  listing  of  all  input  and  output  are  given  in 
Appendix  D. 

The  profile  datasets  written  to  files  on  Units  2  and  3  as  described  in  the  above 
programs  are  the  primary  aircraft  noise  data  used  by  the  NOISEMAP  program  to  generate 
noise  profiles  around  a  specific  airbase.  With  these  programs  these  noise  data  can  be 
tailored  to  the  specific  weather  and  aircraft  flight  and  ground  run-up  power  conditions  for 
that  airbase  which  will,  in  turn,  improve  the  accuracy  of  the  final  noise  contours. 


GENERAL  OVERVIEW  OF  THE  OMEGA  10  PROGRAM 


The  OMEGA  10  program,  hereafter  referred  to  as  simply  the  "program",  is  designed  to 
compute  descriptions  of  the  noise  of  an  aircraft  in  terms  of  A- weighted  overall  sound  level 
(ALX),  tone-corrected  A-weighted  overall  sound  level  (ALTX),  perceived  noise  level  (PNLX), 
tone-corrected  perceived  noise  level  (PNLT),  effective  perceived  noise  level  (EPNLX),  sound 
exposure  level  (SELX),  and  tone-corrected  sound  exposure  level  (SELTX)  as  a  function  of  slant 
distance  to  the  aircraft,  aircraft  power  setting,  and  aircraft  airspeed.  These  noise  measure  data 
(profile  datasets)  are  computed  for  aircraft  flyover  tests  as  outlined  in  AMRL-TR-73-107.(^^  The 
EPNLX,  SELX  and/or  ALX  profile  datasets  are  required  as  input  to  the  NOISEMAP  noise 
exposure  forecast  program. 

To  compute  the  above  integrated  noise  measures  at  the  requested  meteorological 
conditions  and  for  each  power  setting  (PSC)  and  airspeed  requested  on  the  code  sheet,  the 
program  first  reads  in  all  reference  file  datasets  for  aircraft  ACC.  A  description  of  these 
reference  datasets  is  given  in  Appendix  E.  Each  reference  SPL  spectrum  is  then  extrapolated  to 
each  of  the  22  standard  slant  distances  and  to  the  requested  temperature  and  relative  humidity 
from  which  PNLX,  PNLTX,  ALX  and  ALTX  noise  measures  are  determined.  These  four 
extrapolated  noise  measures  as  well  as  the  mean  PNLA,  PNLTA,  ALA,  and  ALTA  data 
computed  from  the  mean  reference  spectrum  and  the  mean  EPNLA,  SELA,  and  SELTA  data 
from  the  reference  dataset  are  used  to  compute  the  PNLX,  PNLTX,  ALX  and  ALTX  single  event 
measures  and  the  EPNLX,  SELX  and  SELTX  integrated  measures  at  each  distance,  power 
setting,  and  airspeed.  The  final  profile  noise  measures  are  then  interpolated  or  extrapolated  from 
these  single  event  noise  data  using  the  rules  established  for  the  specific  data  type. 

The  program  prints  plots  and/or  listings  of  almost  all  of  the  above  input  and  computed 
data.  The  following  is  a  summary  of  the  program  operations: 

(1)  The  program  reads  the  code  sheet  input  data  and  all  reference  datasets  for  aircraft 
ACC  and  initializes  numerous  test  variables. 

(2)  The  program  (subroutine  SETUPD6)  determines  which  reference  file  datasets  are 
required  to  compute  the  requested  profile  data  at  each  power  setting. 

(3)  The  cover  page  is  printed  when  PR>0. 

(4)  The  PNLA,  PNLTA,  ALA  and/or  ALTA  data  are  computed  from  each  reference 
spectrum  used  in  the  analysis. 
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(5)  A  plot  of  the  reference  spectrum  is  printed  on  Page  G  (IPR>0). 

(6)  The  PNLX  through  SELX  (7  measures)  single  event  noise  measures  are  computed 
from  the  22  extrapolated  SPL  spectra  for  air-to-ground  and  ground-to-ground 
propagation  for  the  reference  dataset  power  setting  and  program  reference  airspeed 
(250  knots).  Computer  listings  of  the  SPL  spectra  plus  these  seven  single  event 
measures  at  each  distance  and  plots  of  these  single  event  measures,  all  at  the 
reference  file  power  setting  and  airspeed,  are  printed  only  when  IPR>0  for  the 
NP=0  option. 

(7)  The  interpolation  and  extrapolation  rules  are  then  applied  to  compute  the  final 
profile  data  for  one  or  more  single  event  measures  at  each  power  setting  (PSC). 
Computer  listings  and  plots  of  these  profile  data  are  printed  only  when  IPR  equal 
one. 

(8)  These  EPNLX,  SELX  and/or  ALX  profile  data  are  written  to  a  file  on  Unit  3  when 
IPU  is  greater  than  zero.  A  complete  description  of  these  profile  datasets  is  given  in 
Appendix  F.  Since  the  remaining  PNLTX,  PNLX,  SELTX  and  ALTX  measures  are 
not  used  as  input  to  the  NOISEMAP  program,  they  are  not  written  to  the  file  on 
Units. 

(9)  The  final  summary  page  is  always  printed. 

(10)  After  all  data  are  analyzed  for  all  aircraft  or  when  an  error  occurs,  the  program 
terminates  with  a  call  to  EXIT  where  the  argument  is  one  when  an  error  has 
occurred  or  zero  when  there  is  no  error. 

Examples  of  all  the  above  plots  and  listings  are  presented  with  the  sample  problem  in 
Appendix  C.  When  IPR  is  less  than  one,  only  the  summary  page  plus  any  error  messages  are 
printed  on  the  OUTPUT  file.  For  IPR  greater  than  zero,  all  single  event  measures  are  always 
computed  (all  MEAS  are  set  equal  to  one),  while,  for  BPR  equal  to  zero,  only  the  requested 
measures  (MEAS>0)  are  computed  (default  is  SELX). 

The  CDC  FORTRAN  Extended  (FORTRAN  IV)  computer  language  was  used  for  the 
original  version  of  this  program.  The  common  and  subroutine  features  of  the  language  were  used 
extensively  throughout  the  program  to  save  computer  time  and  memory. 

Over  the  past  ten  years,  this  program  was  converted  from  CDC  FORTRAN  to 
FORTRAN  77  for  the  Perkin-Ehner  and  later  to  Microsoft  and  Lahey  FORTRAN  for  the  PC. 
Version  10.8  is  currently  being  compiled  with  the  Lahey  compiler. 
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The  following  sections  describe  the  detailed  tasks  accomplished  by  the  program.  It  is 
intended  to  document  the  procedures  within  each  subroutine  at  a  level  useful  to  either  a 
programmer  reading  this  while  working  with  the  code  or  a  reader  simply  interested  in  what 
happens  with  a  specific  subroutine.  The  algorithms  used  to  compute  the  noise  measures  are 
described  in  detail  in  the  individual  subroutines. 

Note  that  Appendix  K  contains  a  user's  guide  for  the  MENU  108  program.  The 
MENU  108  program  can  be  used  to  interactively  setup  and  execute  the  OMEGA  10.8  program  as 
well  as  perform  some  post-processing  of  the  noise  measure  profiles. 
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GENERAL  ORGANIZATION  OF  THE  OMEGA  10  PROGRAM 


The  general  organization  of  the  entire  program  is  shown  in  Figure  1 .  The  diagram 
indicates  access  to  the  various  routines  rather  than  program  flow;  for  example,  MAIN  calls 
subroutine  CDIST  which  in  turn  calls  subroutines  CPNL,  CPTC  and  CAL  and  subroutine  CPNL 
calls  function  FNOY.  In  general,  program  flow  is  from  left  to  in  this  diagram.  Function  ICV  is 
called  numerous  times  from  the  MAIN  routine  and  from  subroutines  OUTG,  OUTH,  OUTJ, 
DELTA6,  and  SETUPD6. 

Using  Figure  1  as  a  guide,  this  section  summarizes  in  very  general  terms  the  functions 
performed  by  the  entire  program.  This  is  meant  to  serve  as  an  introduction  for  the  reader  to  the 
functions  of  the  individual  subroutines. 

The  control  routine,  MAIN,  reads  the  job  control  record,  the  aircraft  code  sheet  record, 
and  the  operation  power  information  record  (or  records).  The  required  input  and  output  files  are 
opened  in  routines  CLFILE  or  RDFILE.  For  non-standard  temperature  and  humidity,  subroutine 
ALPH  is  called  to  compute  the  atmospheric  absorption  data.  Numerous  program  and  test 
variables  are  initialized  at  the  beginning  of  the  MAIN  routine. 

Next  the  program  reads  all  reference  datasets  for  aircraft  ACC  from  the  reference  file 
opened  on  Unit  7.  These  datasets  are  checked  for  errors  which  will  result  in  the  aircraft  data 
being  deleted  from  the  job.  At  the  end  of  the  dataset  input,  subroutine  SETUPD6  is  called  to 
determine  which  datasets  are  required  to  compute  the  profile  data  at  each  requested  output  power 
setting  (PSC). 

Subroutines  TITLE  and  BPA  are  called  (for  IPR>0)  to  print  the  cover  page  and  initialize 
plot  arrays  used  in  subroutines  OUTG  and  OUTJ. 

Subroutines  CPNL,  CPTC  and  CAL  are  called  to  compute  the  PNLA,  C  and  ALA  data 
for  each  reference  spectrum;  the  PNLTA  and  ALTA  data  are  then  computed  in  the  MAIN 
routine.  Subroutine  CPNL  uses  function  FNOY  to  compute  the  noy  data  for  each  SPL  level. 
These  data  as  well  as  a  plot  and  listing  of  the  reference  SPL  spectrum  are  printed  on  output  Page 
G  by  subroutine  OUTG  (IPR>0). 

Subroutine  CDIST  is  called  for  each  reference  power  setting  (reference  dataset)  and  each 
type  propagation  (air-to-ground  and  ground-to-ground)  to  extrapolate  the  SPL  reference 
spectrum  to  22  standard  profile  distances  at  a  selected  temperature  and  humidity,  compute 
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Figure  1.  General  Organization  of  the  OMEGA  10  Program. 


PNLX,  PNLTX,  ALX,  ALTX  and  C  for  each  extrapolated  spectrum,  compute  the  SELX,  SELTX 
and  EPNLX  measures  for  each  distance  at  the  reference  file  airspeed,  and  finally  compute  the 
smoothed  EPNLX,  SELTX,  PNLTX  and  ALTX  data.  These  final  EPNLX,  SELX  and  SELTX 
profile  data  are  adjusted  to  the  program  reference  airspeed  (250  knots).  Then  all  seven  profile 
measures  are  stored  in  arrays  PRDA  (air-to-ground)  and  PRDG  (ground-to-ground).  This 
subroutine  calls  subroutines  CPNL,  CPTC,  and  CAL  to  compute  the  above  PNLX,  C  and  ALX 
data.  Subroutine  CPNL  also  uses  function  FNOY  to  compute  the  noy  data  required  by  the  PNLX 
algorithm. 

When  the  program  is  in  the  print  mode  (IPR>0)  with  the  NP=0  option  in  effect,  listings  of 
the  extrapolated  SPL  spectra  and  the  single  event  data  (PNL  etc.)  for  each  distance  are  printed  by 
subroutine  OUTH.  The  SPL  spectra  are  printed  on  output  Pages  H  and  L  for  air-to-ground  and 
ground-to-ground,  respectively;  the  single  event  data  are,  likewise,  printed  on  Pages  I  and  M. 
Subroutine  OUTH  calls  subroutine  HEAD  to  print  each  page  header  block.  Also  plots  of  these 
same  single  event  data  are  printed  by  subroutine  OUTJ.  The  air-to-ground  PNLTX,  PNLX,  ALX 
and  ALTX  single  event  data  are  plotted  on  output  Page  J  and  the  SELTX,  SELX  and  EPNLX, 
single  event  data  on  Page  K.  The  corresponding  ground-to-ground  single  event  data  are  plotted 
on  Pages  N  and  O.  All  these  data  are  computed  at  the  reference  power  setting  and  airspeed. 

When  the  program  is  in  the  print  mode  but  with  NP>0,  the  final  profile  data  for  the  seven 
single  event  measures  for  each  requested  power  setting  (PSC)  and  airspeed  (VX)  are  computed 
by  subroutine  DELTA6.  These  data  are  computed  by  applying  the  a  set  of  rules  to  extrapolate  or 
interpolate  these  data  from  the  reference  data.  Once  again,  listings  and  plots  of  these  single  event 
measures  are  printed  by  subroutines  OUTH  and  OUTJ,  respectively. 

In  the  no  print  mode  (IPR=0)  for  all  NP,  the  program  only  computes  the  single  event 
measures  requested  by  the  MEAS  code  sheet  parameter  plus  all  other  data  required  to  compute 
this  requested  measure.  No  plots  or  listings  are  printed  and  subroutines  OUTH  and  OUTJ  are  not 
called. 


Next,  when  the  program  is  in  the  profile  dataset  mode  (IPU>0),  subroutine  PPFDAT  is 
called  to  write  one  or  more  of  the  EPNLX,  SELX  and  ALX  single  event  measures  to  the  profile 
dataset  file  (Unit  3)  in  the  standard  profile  dataset  format  (Appendix  F).  Finally,  the  summary 
page  is  printed  and  program  control  is  returned  to  label  10  to  begin  the  next  aircraft  analysis  or 
terminate  the  job.  Note  that  the  PNLTX,  PNLX,  SELTX,  and  ALTX  single  event  measures  can 
not  be  written  to  the  profile  dataset  file  on  Unit  3. 
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DEFINITIONS  OF  SYMBOLS  AND  TERMINOLOGY 
USED  IN  THE  OMEGA  10  PROGRAM 


The  symbols  defined  here  are  used  in  this  report  and/or  in  the  OMEGA  10  program 
source  listing.  Some  of  the  symbols  used  in  the  program  are  really  dummy  variables  used  in 
only  one  or  two  routines  and  redefined  in  each  routine;  most  of  these  symbols  are  not  included  in 
this  list  of  symbol  definitions.  Symbols  which  are  arrays  will  be  listed  with  their  array 
dimensions.  Variables  I,  J,  and  L  are  usually  (but  not  always)  used  as  array  subscripts  as 
follows: 


(1)  The  subscript  "I"  is  a  running  index  associated  with  any  one  spectrum  or  standard 
profile  distance. 

(2)  The  subscript  "J"  is  a  running  index  associated  with  any  one  band  in  the  set  of  1/3 
octave  frequency  bands. 

(3)  The  subscript  "L"  is  a  running  index  associated  with  any  one  power  setting. 
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SYMBOL 


DESCRIPTION 


AC 

ACC 

ALA(15) 

ALTA(15) 

ALTX(22) 

ALX(22) 

ATNC(31) 

ATNR(31) 

AW(31) 

BLK 

C0MD(15) 

CRI 

D2X(22) 

DATE 

DATEN 

DELN 

DG 

DRAG(3.15) 

EA(22,24) 

EPNLA(15) 

EPNLX(22) 


Aircraft  name  (10  characters  or  less), 

Aircraft  code  (3  characters  --  see  code  sheet). 

Mean  reference  A-weighted  overall  sound  level  in  dBA. 

Mean  reference  tone-corrected  A-weighted  overall  sound  level  in  dBA. 

Profile  tone-corrected  A-weighted  overall  sound  level  for  each  distance  for  the 
reference  power  setting  and  airspeed  in  dBA. 

Profile  A-weighted  overall  sound  level  for  each  distance  for  the  reference  file  power 
setting  and  airspeed  in  dBA. 

Atmospheric  absorption  coefficients  for  the  OMEGA  10  profile  output  (ITEMP  and 
ERHUM)  in  dB  per  1(K)0  feet. 

Atmospheric  absorption  coefficients  for  standard  day  conditions  (59°F  and  70%)  in 
dB  per  1000  feet. 

A- weighting  coefficients  in  dB. 

A  data  statement  variable  carried  in  OUTC  labeled  common  containing  a  blank 
Hollerith  character.  It  is  used  to  print  variable  format  data  and  initialize  numerous 
program  variables. 

Last  five  characters  of  each  reference  file  dataset  COMDECK  name. 

Comdeck  revision  identifier  (see  code  sheet). 

Factor  to  adjust  the  measure  data  to  the  correct  profile  slant  distance. 

Date  of  computer  run  (see  code  sheet). 

Date  of  computer  run  in  numeric  form  which  is  used  as  part  of  the  data  identification 
code  on  tab  plots  and  in  the  reference  dataset  (see  code  sheet). 

Noise  adjustment  factor  added  to  each  band  in  the  reference  spectrum  (dB). 

Ground-to-ground  propagation  adjustment  factor. 

Drag  configuration  read  from  the  reference  dataset  and  printed  on  output  Page  G. 

Excess  atmospheric  attenuation  in  dB  for  bands  17  to  40  and  distances  200  to  25000 
feet. 

Mean  reference  effective  perceived  noise  level  for  each  power  condition  in  EPNdB. 

Profile  effective  perceived  noise  level  in  EPNdB  for  each  distance  for  the  reference 
power  setting  and  airspeed. 
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SYMBOL 

ERRFLG 

ET(2) 

EXTMX 

FJ 

FL(24,5) 

FM(24,4) 

FREQ3(31) 

HO 

IBNH 

IBNL 

ID 

IMS(15) 

IORD(35) 

IPAGE 

IPR 

IPROP 

IPTC 

IPU 

IRD 

IREQ(15) 


DESCRIPTION 

Logical  variable  which  is  set  to  .TRUE,  when  there  is  an  error  in  an  OMEGA  10  run. 
This  is  used  by  the  NOISEMAP  MCM  to  detect  an  error. 

Engine  type  read  from  the  reference  dataset  and  printed  on  Page  G  and  on  the  second 
profile  dataset  comment  card. 

Maximum  permitted  extrapolation  from  the  reference  value  at  a  slant  distance  of 
1000  feet  for  the  first  air-to-ground  measure  (EXTMX=5.0  dB). 

Constant  used  in  the  perceived  noise  level  computations;  FJ=0.15  for  1/3  octave 
band  data. 

Data  statement  array  used  in  noy  computations  (function  FNOY)  containing  the  band 
sound  pressure  levels  in  dB  given  in  Table  4. 

Data  statement  array  used  in  noy  computations  (function  FNOY)  containing  the 
reciprocals  of  the  slopes  given  in  Table  4. 

Geometric  mean  and  lower  limiting  frequencies  required  to  compute  atmospheric 
absorption  coefficients  for  1/3  octave  data. 

Reference  minimum  slant  range  for  the  power  setting  being  computed. 

Largest  band  number  index  (31  corresponds  to  band  40). 

Initial  band  number  index  (1  corresponds  to  band  10). 

Variable  frequently  used  to  identify  the  SPL  data  to  be  used  in  the  specific 
computations. 

Reference  minimum  slant  range  in  feet  for  each  power  setting. 

Plot  ordinate  scale  values. 

Part  of  page  number  identification. 

Program  print  control  flag  (see  code  sheet). 

Propagation  type  code;  1  for  air-to-ground  and  2  for  ground-to-ground. 

Index  of  band  which  determines  the  tone  correction  for  the  reference  distance  profile 
spectrum. 

Program  profile  dataset  print  (file  on  Unit  3)  control  flag  (see  code  sheet). 

Profile  distance  index  which  corresponds  to  the  reference  minimum  slant  range. 

Program  flag  used  to  indicate  which  reference  data  are  required  to  compute  the 
requested  profile  output. 
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SYMBOL 

IREQC(7,15) 

IRHUM 

ISEQD(7) 

ISRC(31) 

ITEMP 

IV(15) 

rVER 

IVX 

LFLG(15) 

LIM 

MEAS(3) 

MM 

N 

NFFLYl 

NFFLY2 

NP 

NPM 

NR(15) 

0PC(15) 

0PCC(15) 

0PCD(15) 


DESCRIPTION 

Indices  of  reference  power  setting  and  two  slope  power  settings  required  to  compute 
the  profile  output  for  each  requested  power  setting.  See  subroutine  DELTA6  for  the 
complete  definition. 

Reference  relative  humidity  in  percent. 

Data  statement  array  in  subroutine  OUTH  containing  the  indices  of  the  measure  data 
in  array  SENX  in  the  sequence  required  for  printing. 

Integer  value  of  mean  reference  data. 

Reference  surface  temperature  (°F). 

Reference  aircraft  velocity  in  knots  for  each  power  setting. 

Program  version  code. 

Integer  value  of  the  reference  or  profile  aircraft  velocity  in  knots. 

Program  flag  used  to  flag  profile  data  when  the  extrapolation  limit  is  exceeded  or 
when  data  are  omitted  because  of  extrapolation  problems. 

Extrapolation  check  flag.  Extrapolation  is  checked  only  when  LIM=1. 

Program  control  variable  used  to  select  the  profile  measures  for  the  no-print  mode 
(see  code  sheet). 

Maximum  number  of  reference  power  settings  per  aircraft  (MM=15). 

Number  of  reference  file  power  settings  read  from  NOISEFELE  for  aircraft  ACC. 

Character  (*24)  variable  used  to  store  the  first  of  two  reference  dataset  files  name. 

Character  (*24)  variable  used  to  store  the  second  reference  dataset  file  name. 

Number  of  output  power  settings  for  which  profile  data  are  requested  for  this 
aircraft. 

Maximum  NP  permitted  by  the  program  (NPM=15). 

Number  of  runs  (measure  locations)  averaged  to  obtain  the  given  reference  data  for 
each  power  setting  (the  reference  data  were  computed  in  the  OMEGA  6  program). 

Operation  power  code  for  each  reference  power  setting. 

Operation  power  code  for  each  profile  output  power  setting. 

Data  statement  array  containing  the  default  profile  output  operation  power  codes. 


12 


SYMBOL 


DESCRIPTION 


0PCR(15) 

OTC 

P(2,15) 

PC(2,15) 

PNLA(15) 

PNLTA(15) 

PNLTX(22) 

PNLX(22) 

PRDA(22,15,7) 

PRDC(22,7,2) 

PRDC(22) 

PRDG(22,15,7) 

PRDI(22,15) 

PS(2.15) 

PSC(15) 

PSCF(15) 

PSIF(15) 

PSU 

PTC 

PV 

RV 


Operation  power  code  for  the  reference  data  from  which  the  OPCC  data  are 
computed  and  the  operation  description  taken. 

Operation  type  code. 

Power  setting  description  data  for  each  reference  power  setting. 

Power  setting  description  data  for  each  profile  output  power  setting. 

Mean  reference  perceived  noise  level  in  PNdB  for  each  power  setting. 

Mean  reference  tone-corrected  perceived  noise  level  in  PNdB  for  each  power  setting. 

Profile  tone-corrected  perceived  noise  level  for  the  reference  power  setting  and 
airspeed  in  PNdB. 

Profile  perceived  noise  level  for  the  reference  power  setting  and  airspeed  in  PNdB. 

Profile  data  for  EPNLX,  SELTX, ...,  ALX  for  air-to-ground  propagation  for  each 
reference  power  setting  and  adjusted  to  the  program  reference  airspeed  (RV). 

Final  EPNLX,  SELTX, ...,  ALX  profile  data  for  air-to-ground  and  ground-to-ground 
for  one  power  setting  (computed  by  subroutine  DELTA6). 

Array  used  in  subroutine  DELTA6  to  compute  the  profile  data  for  one  specific 
measure  and  type  propagation. 

Same  as  PRDA  above  except  for  ground-to-ground  data. 

Array  equivalent  (subroutine  DELTA6  argument)  to  specific  elements  of  array 
PRDA  or  PRDG.  Contains  all  available  profile  reference  data  for  one  measure  and 
one  type  propagation. 

Reference  power  setting  data  (value  and  units  in  character  format)  for  each  reference 
operation  power  code. 

Profile  output  power  setting  for  each  requested  operation  power  code  in  character 
format  (see  code  sheet). 

Same  as  PSC(15)  except  in  numeric  format. 

Reference  power  setting  in  numeric  format. 

Power  setting  units  for  this  aircraft. 

Tone  correction  in  dB. 

Profile  version  code  (see  code  sheet). 

Program  reference  airspeed  in  knots  (RV=250). 
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SYMBOL 

SELA(15) 

SELTA(15) 

SELTX(22) 

SELX(22) 

SENX(22,7) 
SOURCE(2.15) 
SPLA(  15,31) 
SPLX(22,31) 
SR(  15,40) 

SX(22) 

THETA(15) 

TYPE(15) 

VFCT 

VREF 

VX(15) 


DESCRIPTION 

Mean  reference  sound  exposure  level  in  dB  for  each  power  setting. 

Mean  reference  tone-corrected  sound  exposure  level  in  dB  for  each  power  setting. 

Profile  tone-corrected  sound  exposure  level  in  dB  for  a  specific  reference  power 
setting  and  airspeed. 

Profile  sound  exposure  level  in  dB  for  a  specific  reference  power  setting  and 
airspeed. 

Array  equivalent  to  arrays  PNLTX(22)  through  SELX(22)  in  blank  common. 
Identifies  the  program  and  date  from  each  reference  dataset. 

Mean  reference  sound  pressure  level  in  dB  for  each  power  condition. 

Profile  sound  pressure  level  in  dB. 

Mean  reference  data  for  each  power  condition;  equivalent  to  arrays  SPLA  through 
CA  in  blank  common  (see  subroutine  CDIST)  . 

Distance  data  in  feet  for  the  22  profile  distances. 

Mean  reference  angle  of  maximum  radiation  in  degrees  for  each  power  condition 
(directivity  angle). 

Character  variable  used  to  store  the  interpolation  code  (I,  N,  or  P)  for  each  reference 
dataset. 

Variable  used  to  adjust  aircraft  airspeed  from  reference  dataset  airspeed  to  program 
reference  airspeed  or  from  program  reference  to  final  profile  airspeed. 

Floating  point  value  of  the  aircraft  reference  airspeed  (IV). 

Aircraft  airspeed  in  knots  for  each  profile  power  condition. 
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DETAILED  DESCRIPTION  OF  THE  OMEGA  10  PROGRAM 


This  section  discusses  the  MAIN  routine  and  each  subroutine  in  the  OMEGA  10 
program  Procedures  within  most  routines  are  documented  at  a  level  useful  to  a  programmer 
reading  this  while  working  with  the  code  or  a  reader  interested  in  what  happens  within  a  specific 
subroutine.  Most  routines  contain  numerous  comments  which  should  be  very  helpful  in 
following  the  code. 

The  program  algorithms  and  I/O  are  discussed  in  the  routines  in  which  they  are  coded. 

The  program  code  sheet  in  Appendix  A  and  the  sample  problem  in  Appendix  C  are  referenced  to 
simplify  the  description  of  the  input  and  output.  The  more  complex  subroutines  are 
supplemented  by  flow  charts  drawn  from  the  point  of  view  of  function  performed  rather  than 
block  instructions. 

COMMON  VARIABLES 

Extensive  use  is  made  of  common  in  the  program  for  communications  between  the 
various  routines.  Many  of  the  storage  locations  in  blank  common  are  used  in  different  ways  or 
with  different  variable  names  throughout  the  program.  This  was  done  to  reduce  the  memory 
required  to  run  the  program  on  the  old  CDC  computers.  Several  of  the  large  arrays  are  included 
in  blank  common  rather  than  labeled  common  because  on  the  CDC  CYBER  computers  a  large 
blank  common  reduces  the  total  memory  required  to  load  and  execute  the  program.  The 
variables  used  in  labeled  common  are  usually  of  a  similar  type  and/or  used  in  many  of  the  same 
subroutines. 

The  variables  assigned  to  blank  and  labeled  common  in  the  MAIN  routine  and  the  total 
common  length  are  listed  in  Table  1.  The  subroutines  in  which  the  labeled  common  are  used  are 
listed  in  Table  2.  All  blank  and  labeled  common  are  included  in  the  MAIN  routine.  All  common 
variables  are  defined  in  the  complete  list  of  symbols  in  the  previous  section  of  this  report.  The 
blank  and  labeled  common  are  described  in  the  following  paragraphs. 

Blank  Common 

The  24,420  bytes  used  by  blank  common  in  the  main  routine  is  the  maximum  required  in 
any  routine  in  the  program.  Almost  all  routines  use  some  blank  common  but  only  a  few  require 
the  24,420  bytes.  The  variable  names  assigned  to  blank  common  vary  throughout  the  program. 
Most  variables  are  defined  to  communicate  with  several  subroutines  and  then  redefined  for  the 
next  series  of  routines. 
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TABLE  1 

MAIN  ROUTINE  VARIABLES  IN  BLANK  AND  LABELED  COMMON 


Length 
in  Bytes 


Labeled  Common 


Blank  Common 

HEADC 

COMPC 

OUTC 

CHAR 

ERRORS 

ffiNL 

IPAGE 

IV(15) 

ISC(9) 

AC 

ERRFLG 

EBNH 

IVX 

IMS(15) 

IORD(35) 

ACC 

L 

ITEMP 

PSIF(15) 

BLK 

SR(  15,40) 

IRHUM 

PSCF(15) 

CRI 

NR(15) 

rVER 

IREQC(7.15) 

DASH 

ISRC(31) 

OTC 

VX(15) 

DATC 

SPLX(22,31) 

SX(22) 

DATEN 

PNLTX(22) 

ATNC(31) 

DOT 

PNLX(22) 

ATNR(31) 

ET(2) 

AL'rX(22) 

DELN 

OPC(15) 

ALX(22) 

EPTC 

OPCC(15) 

EPNLX(22) 

EPROP 

OPCD(15) 

SELTX(22) 

MEAS(3) 

OPCR(15) 

SELX(22) 

ORD(43) 

PRDA(22,15,7) 

P(2,15) 

PRDG(22,15,7) 

PC(2.15) 

PP(84.3) 

PS(2,15) 

PSC(15) 

PSU 

PV 

X 

24420 

24 

1080 

176 

1331 

4 
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TABLE  2 

SUBROUTINES  CONTAINING  THE  LABELED  COMMON 


HHADC 

COMPC 

OUTC 

CHAR 

ERRORS 

HEAD 

CDIST 

IPA 

HEAD 

SETUPD6 

OUTG 

HEAD 

OUTG 

IPA 

CLFILE 

OUTJ 

OUTG 

OUTH 

OUTG 

RDFBLE 

PPFDAT 

OUTJ 

OUTJ 

CDIST 

SUMRY 

PPFDAT 

BLOCK  DATA 

OUTH 

tupg 

SETUPD6 

OUTJ 

BLOCK  DATA 

SUMRY 

PPFDAT 

BLOCK  DATA 

TITPG 

SETUPD6 

DELTA6 

SUMRY 

BLOCK  DATA 
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CHAR  Common 

This  CHAR  labeled  common  was  added  to  the  program  when  the  program  was  converted 
to  the  FORTRAN  77  and  rehosted  to  an  IBM  compatible  PC.  All  variables  in  this  labeled 
common  are  character  variables  which  were  previously  used  in  the  HEADC,  COMPC  and  OUTC 
labeled  common.  The  character  and  numeric  variables  were  separated  at  that  time  with  all 
character  variables  being  put  into  the  CHAR  labeled  common.  No  effort  was  made  to  group 
these  variables  in  any  other  way. 

COMPC  Common 

These  variables  are  a  combination  of  output  variables  and  variables  required  to  compute 
the  profile  data  or  both.  The  variable  names  are  the  same  throughout  the  program. 

ERRORS  Common 

The  error  flag  variable  stored  in  this  labeled  common  is  used  only  in  the  main  routine  and 
in  subroutines  RDFDLE,  CLFDLE  and  SETUPD6.  It  was  added  in  OMEGA  10,  Version  8  to  flag 
errors  as  required  in  the  NOISEMAP  MCM. 

HEADC  Common 

These  variables  are  primarily  output  variables  required  to  identify  the  data  being 
analyzed.  Most  of  these  data  are  printed  in  the  page  header  blocks,  the  plot  identification  blocks, 
and  the  output  COMDECK  comment  records.  The  variable  names  are  the  same  throughout  the 
program. 

OUTC  Common 

These  variables  are  required  to  set  up  and  print  the  listings  and  plots  in  subroutines 
OUTG,  OUT!  and  OUTH.  If  all  listing  output  are  deleted  (e.g.,  in  some  future  preprocessor 
version  of  this  program),  this  common  statement  will  no  longer  be  required.  These  storage 
locations  are  defined  the  same  throughout  the  program. 

MAIN  PROGRAM 

MAIN  is  the  executive  routine  for  the  entire  OMEGA  10  program.  Its  principal  function 
is  to  call  the  subroutines  required  to  compute  the  aircraft  analysis.  However,  it  also  reads  all  the 
aircraft  setup  (code  sheet)  data,  initializes  numerous  program  and  aircraft  variables,  and  reads 
and  checks  the  reference  file  (Unit  7)  input  data.  The  MAIN  routine  is  discussed  in  the  following 
paragraphs  using  the  program  listing  in  Appendix  I  and  the  flowchart  in  Figure  2  as  a  guide. 
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Method 

The  first  segment  of  the  MAIN  routine,  from  the  beginning  to  label  150,  performs  the 
following  data  input  and  initialization  steps: 

(1)  The  Lahey  FORTRAN  routine  GETCL  is  called  to  get  the  execution  statement 
command  line  arguments.  If  the  command  line  is  empty,  subroutine  RDFILE  is 
called  to  query  the  operator  for  the  setup  file  name  and  the  profile  dataset  file  name. 
When  the  command  line  is  empty,  the  list  output  is  always  written  to 
OMEGAIOS.OUT  and  the  reference  databases  are  always  assigned  to  file  names 
NTvI63FLYl  and  NM63FLY2.  These  list  and  reference  database  files  will  change 
when  the  program  and  database  versions  change.  If  the  command  line  is  not  empty, 
subroutine  CLFELE  is  called  to  parse  the  five  file  names  on  the  command  line.  The 
command  line  must  contain  all  five  I/O  file  names  as  described  in  the 
documentation  for  subroutine  CLFILE. 

(2)  Numerous  program  and  aircraft  variables  are  initialized.  The  SX  array  is  initialized 
to  the  22  standard  profile  distances: 

SX,  =antUog  feet 

where  I  is  the  standard  distance  index  (I  =  1  to  22). 

(3)  The  job  control  record  is  read  and  default  values  are  set  when  they  apply.  For  IPR 
equal  to  zero  or  blank  (no-print  mode),  the  IPU  flag  is  always  set  equal  to  one  and, 
if  all  MEAS  are  blank,  MEAS3  is  set  equal  to  one  (SELX).  For  EPR  equal  one 
(print  mode),  all  MEAS  flags  are  set  equal  to  one  (all  three  measures  are  always 
computed). 

(4)  The  aircraft  setup  records  for  aircraft  ACC  are  read  and  additional  aircraft  default 
parameters  are  set. 

(5)  The  atmospheric  absorption  coefficients  are  computed  for  non-standard  temperature 
and  relative  humidity.  The  coefficients  for  standard  conditions  (59°F  and  70%)  are 
stored  in  the  ATNR  data  statement  array. 

(6)  The  reference  dataset  file  is  opened  and  searched  for  all  reference  datasets  for 
aircraft  ACC.  All  data  for  aircraft  ACC  must  be  back-to-back  in  the  reference  file. 
A  maximum  of  fifteen  (MM)  are  read  and  stored  in  the  program.  If  several  datasets 
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with  the  same  operation  power  code  (OPC)  are  in  the  reference  file,  the  last  dataset  for 
this  OPC  is  stored  in  the  program.  These  reference  datasets  are  checked  for  errors  in 
aircraft  code,  operation  power  code,  and  record  sequence  number;  an  error  will  abort  this 
aircraft  analysis  (not  the  entire  job).  It  should  be  noted  here  that  the  current  version  of 
OMEGA  10  (Version  8)  defines  two  reference  dataset  files  in  variables  NFFLYl  and 
NFFLY2.  The  file  named  in  variable  NFFLYl  is  opened  for  military  aircraft  which 
includes  aircraft  codes  less  than  200  or  greater  than  349.  The  file  named  in  variable 
NFFLY2  is  opened  for  civil  aircraft  with  codes  from  200  to  349.  These  aircraft  code 
numbers  as  well  as  the  number  of  reference  dataset  files  may  change  in  future  versions  of 
this  program. 

After  all  input  data  are  read  and  checked,  the  program  calls  subroutine  SETUPD6  to 
apply  the  interpolation/extrapolation  rules  in  the  selection  of  the  one,  two,  or  three  reference 
datasets  required  to  compute  the  profile  output  dataset  for  each  power  setting  (PSC).  The  indices 
of  these  reference  data  for  each  PSC  are  stored  in  array  DREQC.  Also  flag  array  IREQ  is  set 
equal  to  one  when  the  corresponding  reference  dataset  is  required  to  compute  profile  data.  The 
IREQC  and  IREQ  arrays  are  defined  in  greater  detail  in  comment  statements  in  subroutines 
SETUPD6  and  DELTA6. 

For  IPR  greater  than  zero,  the  cover  page  is  printed  by  subroutine  TfTPG  and  several  plot 
arrays  used  to  print  the  plots  in  subroutines  OUTG  and  OUTJ  are  initialized  by  subroutine  BPA. 

The  label  600  loop  performs  the  following  operations  for  each  of  the  N  reference  dataset 
power  settings  (IREQ>0)  required  to  compute  the  requested  profile  noise  data  at  the  PSC  power 
settings: 

(1)  Subroutines  CPNL,  CPTC  and  CAL  are  called  to  compute  the  perceived  noise  level 
(PNLA),  tone  correction  (C  or  PTC),  and  A-weighted  overall  sound  level  (ALA) 
from  the  reference  SPL  spectrum  for  the  L^^  power  setting.  The  PNLA  data  are 
stored  in  SRl  33  and  the  ALA  are  in  SRl,35  The  tone-corrected  perceived  noise 
level  (PNLTA)  is: 

SRl.34=  SRl,33+PTC  PNdB 

and  the  tone-corrected  A-weighted  overall  sound  level  (ALTA)  is: 

SRl.36=  SRl.35+PTC.  dBA 
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(2)  Subroutine  OUTG  prints  (for  IPR>0)  a  plot  and  listing  of  the  SPL  reference 
spectrum  and  prints  the  data  computed  in  step  1  above,  the  EPNLA,  SELA,  SELTA 
and  THETA  from  the  reference  dataset,  and  numerous  additional  identification 
parameters,  all  for  the  reference  power  setting. 

(3)  The  540  loop  calls  subroutine  CDIST  to  compute  all  SPL  and  single  event  data  for 
air-to-ground  and  ground-to-ground  propagation.  Subroutine  CDIST  extrapolates 
the  reference  SPL  spectrum  to  each  profile  distance  and  then  computes  all  single 
event  data  for  each  distance.  Subroutine  OUTH  prints  listings  of  these  SPL  and 
single  event  data  and  subroutine  OUTJ  prints  plots  of  the  single  event  measure 
versus  distance  data  (both  only  for  the  IPR>0  for  the  NP=0  option).  The  data 
computed  and  printed  here  are  at  the  reference  power  setting  and  airspeed  (L^h 
power  setting). 

The  label  700  loop  performs  the  following  operations  for  each  of  the  NP  output  power 
settings  (PSC): 

(1)  Subroutine  DELTA6  is  called  to  compute  the  final  profile  data  for  all  single  event 
measures  when  IPR>0  and  for  each  requested  (MEAS>0)  measure  when  IPR=0. 

The  indices  of  the  reference  profile  data  required  to  interpolate  or  extrapolate  these 
profile  data  are  given  in  array  IREQC  which  was  evaluated  in  subroutine 
SETUPD6.  These  final  profile  data  are  stored  in  array  PRDC.  After  all  profile  data 
are  computed  for  the  L^h  power  setting  (end  of  loop  620),  the  LFLG  program  flag  is 
checked.  If  LFLG  is  less  than  zero,  no  profile  data  were  computed  for  this  L^^ 
power  setting  and  the  program  skips  down  to  the  end  of  the  700  loop. 

(2)  For  IPRR  equal  to  one,  subroutines  OUTH  and  OUTJ  are  called  to  print  and  plot 
these  final  EPNLX  through  ALTX  profile  data.  Note  that  whenever  profile  data  are 
printed  all  seven  measures  are  always  computed. 

(3)  For  IPU  greater  than  zero,  the  requested  (MEAS>0)  EPNLX,  ALX  and  SELX 
profile  measures  are  written  to  the  file  on  Unit  3  by  subroutine  PPFDAT.  In  the 
print  node  (IPR>0),  these  profile  data  are  usually  not  written  to  Unit  3.  A 
description  of  the  format  and  content  of  these  flight  noise  profile  datasets  is 
presented  in  Appendix  F. 

After  completing  the  label  700  loop  operations,  subroutine  SUMRY  is  called  to  print  the 
final  aircraft  summary  page  which  lists  numerous  test  identification  parameters,  identifies  all 
reference  datasets  read  for  this  aircraft,  and  lists  the  COMDECK  names  of  all  reference  datasets 
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required  to  interpolate  or  extrapolate  each  profile  dataset.  This  summary  page  is  always  printed 
and  is  not  controlled  by  the  BPR  print  flag.  Program  control  then  returned  to  statement  label  10 
to  perform  the  next  aircraft  analysis  or  to  terminate  the  job.  The  OMEGA  10.8  job  will  terminate 
with  a  call  to  EXTTfl)  if  an  error  has  occurred  (ERRET^G  =  TRUE)  or  a  call  to  EXIT(O)  when  the 
job  terminates  without  an  error  (ERRFLG  =  FALSE). 

SUBROUTINE  ALPH(REL,TEMP^BC) 

This  subroutine  is  called  from  the  MAIN  routine  to  compute  the  atmospheric  absorption 
coefficients  for  non-standard  output  temperature  (TEMP)  and  relative  humidity  (REL)  for 
frequency  bands  10  to  40.  The  coefficients  for  bands  10  to  16  are  always  zero.  The  coefficients 
for  standard  day  conditions  (59°F,  70  percent)  are  stored  in  the  MAIN  routine. 

Subroutine  Arguments 

The  subroutine  arguments  are  defined  as  follows: 

(1)  REL  and  TEMP  are  the  relative  humidity  and  temperature  (°F)  for  which  the  ABC 
data  are  computed. 

(2)  ABC  is  the  absorption  coefficient  array  defined  in  this  subroutine 

Data  Statement  Arrays 

Array  FREQS  contains  the  geometric  mean  and  lower  limiting  frequencies  for  1/3  octave 
data.  X  and  Y  contain  the  normalized  absolute  humidity  (hnormalized)  normalized 
molecular  absorption  coefficient  (Otnormalized)  data  given  in  SAE  ARP  866A.(^) 

Atmospheric  Absorption  Coefficient  Algorithm 

The  atmospheric  absorption  coefficients  for  1/3  octave  frequencies  are  computed  as 
described  below. 

(1)  The  absolute  humidity  statement  function  is: 

F(TEMP,REL)  =  {(a)(REL)}{antilog[(b)(TEMP)+(c)(TEMp2)+(d)(TEMP3)J} 

where 

a=  1.064764  X  10-2 
b  =  2.288074  x  10-2 
c  =  -9.589  X  10-5 
d  =  3.0  X  10-2 
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(2)  The  normalized  absolute  humidity  (hnormalized)  is  computed  as  a  function  of 
frequency. 

HN^FjTEMP.REL) 

fFREQY^ 

I  1010  J 

where 

FREQ  =  the  geometric  mean  or  lower  limiting  frequency  for  the  1/3  octave  frequency. 

(3)  A  quadratic  Aitkin  interpolation  function  is  used  to  compute  the  normalized 
molecular  absorption  coefficient  (ALN)  for  a  given  normalized  absolute  humidity 
(HN). 

(4)  Then  the  atmospheric  absorption  coefficient  (ABC)  as  a  function  of  frequency  is: 

ABC  =  (FREQ)(Fr  1  )(ALN)+(FT2)(FREQ)205  dB 

where 

FTl 
FT2 

(FREQ)(FT1) 

(FT2)(FREQ)205 

FUNCTION  ATKN(X,  Y,N,K,XI) 

This  function  is  a  general  AJTKEN  interpolation  function,  used  by  subroutine  ALPH  to 
compute  the  normalized  molecular  absorption  coefficient.  ATKN  was  obtained  from  the  ASD 
computer  center  library  (old  IBM  7094  library).  Since  this  is  a  common  interpolation  function 
defined  in  most  numerical  methods  texts,  no  additional  description  will  be  given.  The  function 
arguments  are  defined  in  comment  statements  in  the  listing. 

FUNCTION  ICV(R) 

This  function  is  used  in  numerous  routines  to  convert  variable  DR  from  floating  point  to 
integer.  The  resulting  integer  is  rounded  up  when  the  fractional  part  is  greater  than  or  equal  to 
0.5.  This  method  gives  more  consistent  and  conservative  noise  measure  data  than  would  be 
obtained  by  simple  truncation.  In  this  PC  version  of  OMEGA  10,  ICV  is  set  equal  to  the  MINT 
function.  Function  ICV  could  be  replaced  by  MINT  everywhere  in  the  program. 


=  (3.788785  x  10-3)(antilog[4.6833333  x  lO’^  x  TEMP]) 
=  (2.4931591  X  10-8)(antilog[6.33  x  IQ-^  x  TEMP]) 

=  maximum  molecular  absorption  coefficient. 

=  classical  absorption  coefficient. 
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SUBROUTINE  HEAD(IP) 


This  subroutine  is  called  from  subroutine  OUTH  to  print  the  header  block  on  output 
Pages  H,  I,  L  and  M.  The  IP  subroutine  argument  identifies  the  page  for  which  the  header  block 
is  requested.  These  header  blocks  are  108  columns  wide  (printer  columns  1 1  to  1 19).  Each 
header  block  contains  the  five  categories  of  data  information  described  below: 

( 1 )  The  identification  block  located  in  the  upper  right  comer  identifies  the  specific 
OMEGA  program,  the  aircraft  code,  the  operation  type  and  power  codes  (OPS),  the 
profile  version  code,  the  date  of  the  computer  ran,  and  the  page  number  of  that 
particular  computer  printout.  These  identification  data  enable  the  AL/OEBN  to 
locate  and  track  specific  results  in  our  data  bank  and  reconstruct  the  source  of  the 
data  and  all  processing  parameters  and  computations. 

(2)  The  table  title  printed  across  the  top  describe  the  measures  or  type  of  data  presented 
on  that  page. 

(3)  The  aircraft  block  identifies  the  specific  aircraft  (e.g.,  C-141). 

(4)  The  operation  block  presents  the  operation  power  description,  the  power  setting, 
and  the  airspeed. 

(5)  The  meteorology  block  gives  the  temperature  (°F),  relative  humidity  (percent),  and 
Delta  N  (dB).  Note  that  Delta  N  is  not  a  meteorology  related  constant. 

SUBROUTINE  IPA 

This  subroutine  is  called  from  the  MAIN  routine  to  initialize  the  ORD,  PD,  PB,  and  PM 
arrays  used  to  print  the  tab  plots  in  subroutines  OUTG  and  OUTJ.  The  content  of  these  arrays  is 
obvious  in  the  coding. 

SUBROUTINE  OUTG(DRAG, SOURCE) 

This  subroutine  is  called  from  the  MAIN  routine  to  print  the  plot  and  listing  of  the 
reference  SPL  spectrum  versus  1/3-octave  frequency  for  the  power  setting.  The  PNLA, 
PNLTA,  AL,  ALTA,  and  C  data  computed  from  this  mean  spectrum  and  THETA,  EPNLA, 
SELA  and  SELTA  from  the  reference  dataset  as  well  as  various  operation  and  identification 
parameters  are  also  printed  below  the  plot. 
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Method 


The  integer  SPL  data  are  computed  and  stored  in  array  ISRC.  The  maximum  annotated 
abscissa  scale  value  (MX)  is  the  nearest  multiple  of  ten  less  than  seven  plus  the  maximum  SPL. 
The  minimum  scale  value  (MN)  is  70  dB  less  than  MX.  The  abscissa  scale  is  defined  in  array 
ISC.  The  actual  plotted  maximum  and  minimum  are  two  greater  than  MX  and  two  less  than  MN. 

The  SPL  data  are  scaled,  set  up  in  the  print  array  (PP),  and  printed  in  the  label  150  loop. 
The  remainder  of  the  subroutine  completes  the  plot  and  prints  the  PNLA,  PNLTA,  etc.  data 
below  the  plot.  The  detailed  format  of  these  Page  G  data  can  best  be  observed  by  consulting  the 
sample  problem  in  Appendix  C.  There  should  be  sufficient  comment  statements  in  the 
subroutine  listing  to  identify  the  data  being  printed. 

SUBROUTINE  CDIST(IRD,RV) 

This  subroutine  is  called  from  the  MAIN  routine  to  compute  the  single  event  noise  profile 
data  (e.g.,  SELX  versus  distance  for  22  distances  from  200  to  25,(XX)  feet)  for  air-to-ground 
(IPROP=l)  or  ground-to-ground  (IPROP=2)  propagation  for  the  L^  reference  power  setting. 
These  data  are  also  adjusted  to  the  output  temperature  and  relative  humidity  and  to  the  program 
reference  (RV)  airspeed. 

The  subroutine  argument  IRD  is  the  index  of  the  standard  distance  which  is  within  one 
percent  of  the  reference  distance.  IRD  is  determined  in  this  subroutine  and  returned  to  the  MAIN 
routine. 

The  EA  data  statement  array  contains  the  excess  attenuation  data  for  frequency  bands  17 
through  40  and  distances  200  to  25000  feet.  These  data  are  presented  in  Table  3  which  was  taken 
from  AAMRL-TR-84-017.(^)  This  table  also  contains  the  data  for  bands  10  through  16  which 
are  the  same  as  for  band  17.  They  are  required  to  compute  the  ground-to-ground  propagation 
data. 

The  flowchart  in  Figure  3  provides  a  brief  description  of  the  data  computed  by  this 
subroutine. 

Extrapolated  SPL  Spectra 

The  first  segment  of  this  subroutine  (label  30  loop)  extrapolates  the  reference  SPL 
spectrum  to  each  of  the  22  profile  distances  defined  in  array  SX.  These  22  SPL  spectra  are  used 
to  compute  the  single  event  profile  data  for  air-to-ground  propagation.  The  extrapolation 
algorithm  is: 
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TABLE  3 

EXCESS  SOUND  ATTENUATION  MODEL  USED  IN  THE  OMEGA  10  AND  OMEGA  11  PROGRAMS  (in  dB) 
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SPLKu  =  SPLAj^j  -  ).  _ 20 log(5X,  /  HO)  +  DELN  dB 


where 


SPLXij 

SPLAlj 


ATNCj 


ATNRj 


SNTH 


=  the  sound  pressure  level  in  dB  for  the  distance  and  band. 

=  the  mean  sound  pressure  level  in  dB  for  the  reference  power  setting  and 
the  band. 


=  antilog 


7  +  22 


which  is  the  profile  distance  in  feet. 


THETAl  = 


DELN 


=  the  atmospheric  absorption  coefficients  in  dB  per  1000  feet  for  the  profile 
output  temperature  (ITEMP)  and  humidity  (IRHUM). 

=  the  atmospheric  absorption  coefficients  in  dB  per  1000  feet  for  standard  day 
conditions. 

=  (1000)  (SIN  [THETAl]). 

=  the  directivity  angle  in  degrees  for  the  power  setting. 

=  the  reference  minimum  slant  range  in  feet  for  the  power  setting. 

=  the  constant  dB  level  added  to  all  SPL  bands  for  the  power  setting. 


The  following  distance  function  used  later  in  this  subroutine  to  adjust  the  integrated 
single  event  data  for  changes  in  slant  range  is  also  computed  in  this  segment  of  the  program: 


D2Xi  =  6  log  (SXj/HO) 


Adjust  the  SPL  Spectra  for  Ground-to-Ground  Propagation 


The  SPL  spectra  computed  for  air-to-ground  propagation  are  adjusted  for  ground-to- 
ground  propagation  after  all  air-to-ground  data  are  computed  for  the  L^^  reference  power  setting: 


where 


EAi,i 


SPLXy  =  SPLXij -EAli 


=  the  distance  index  defined  for  the  22  distances. 

=  the  frequency  band  index  defined  for  frequency  bands  10  through  16. 

=  excess  attenuations  in  dB  for  the  distance  and  for  frequency  band  17. 


SPLXij  =  SPLXij -EAij 
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where  I  is  defined  for  all  22  distances  and  J  is  defined  for  bands  17  through  40.  The  excess 
attenuation  is  constant  for  frequency  bands  less  than  17  and  greater  than  36. 

A  constant  ground-to-ground  adjustment  factor  (DG)  is  also  subtracted  from  the  air-to- 
ground  SPL  for  all  distances  and  all  frequency  bands: 

SPLXij  =  SPLXi  j  -  DG 
This  DG  adjustment  factor  is  currently  0  dB. 

Single  Event  Noi.se  Data 

The  statement  label  150  loop  computes  all  single  event  noise  data  for  each  of  the  profile 
distances.  Subroutines  CPNL,  CPTC,  and  CAL  are  called  to  compute  the  perceived  noise  level 
(PNLX),  tone  correction  (PTC)  and  the  A-weighted  overall  sound  level  (ALX),  respectively,  for 
each  distance  spectrum.  If  PNLX  data  are  missing  for  distances  beyond  the  second  distance,  the 
PNLX  data  are  extrapolated  as  follows: 

PNLXi  =  (2)  (PNLXi. i)  -  PNLXi.2  PNdB 

The  tone-corrected  perceived  noise  level  (PNLTX)  and  tone-corrected  A-weighted 
overall  sound  level  (ALTX)  for  the  distance  are: 

PNLTXi  =  PNLXi  +  PTC  PNdB 

ALTXi  =  ALXi  +  PTC  dBA 

The  effective  perceived  noise  level  (EPNLX),  sound  exposure  level  (SELX),  and  tone- 


corrected  sound  exposure  level  (SELTX)  are: 

EPNLXi  = 

EPNLAl  +  PNLTXi  '  PNLTAl  +  D2Xi 

EPNdB 

SELXi  = 

SELAl  +  ALXi  -  ALAl  +  D2Xi 

dB 

SELTXj  = 

SELTAl  +  ALTXi  -ALTAl  +  D2Xi 

dB 

where  EPNLAl,  SELAl,  and  SELTAl  mean  data  from  the  reference  dataset  for  the  L^ 
power  setting.  The  PNLTAl,  ALAl,  “d  ALTAl  were  computed  from  the  reference  SPL 
spectrum  for  the  L*^  power  setting. 
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Smooth  Tone-Corrected  Data 

The  last  segment  of  this  subroutine  (statement  label  220  loop)  computes  the  smoothed 
ALTX,  PNLTX,  SELTX,  and  EPNLX  data  for  each  distance  (I)  and  adjusts  these  data  to  the 
program  reference  airspeed: 


ALTXi 

=  ALXi  +  (C2)  (D3) 

dBA 

PNLTXi 

=  PNLXi-J-(C2)(D3) 

PNdB 

SELTXi 

=  SELXi  +  (Cl)(D3) 

dB 

EPNLXi 

=  EPNLXi  +  (C1)(D3)  -  C3i 

EPNdB 

where 

Cl 

SELTXird  -  SELXird 

dB 

C2  = 

ALTXird  -  ALXird 

dBA 

IRD  = 

the  index  of  the  reference  distance 

C3i  = 

SELTXi  -  SELXi 

dB 

D3  =  1.0  for  distances  200  to  3,150  feet 

D3  =  (0.2)  (18-1)  for  distances  4,000  to  8,000  feet  (1=14  to  17). 

D3  =  0.0  for  distances  10,000  to  25,000  feet. 

The  SELTXj  used  to  compute  C3j  above  is  the  unsmoothed  SELTXj  defined  in  the  previous 
section;  that  is,  C3i  is  the  tone  correction  for  the  unsmoothed  SELTX  and  EPNLX  data. 

All  EPNLX,  SELTX,  and  SELX  data  for  the  L^h  reference  power  setting  are  adjusted  to 
the  program  reference  airspeed  and  stored  in  arrays  PRDA  or  PRDG  as  required  in  subroutine 
DELTA6.  The  four  remaining  single  event  measures,  PNLTX  through  ALX,  are  also  stored  in 
arrays  PRDA  and  PRDG,  but  without  the  airspeed  adjustment  which  does  not  apply  to  these 
nonintegrated  measures.  It  should  also  be  noted  here  that  all  seven  measures  are  stored  in  blank 
common  arrays  PNLTX  through  SELX  without  the  airspeed  adjustment. 


For  air-to-ground  propagation  (IPROP=l): 


PRDAi^L,!  = 

EPNLXi  -VFCT 

EPNdB 

PRDAi^L,2  = 

SELTXi  -  VFCT 

dB 

PRDAi,l,3  = 

SELXi  -  VFCT 

dB 

PRDAi,l,4  = 

PNLTXi 

PNdB 

PRDAi^L,5  = 

PNLXi 

PNdB 
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dBA 
dBA 

and  for  ground-to-ground  propagation  (IPROP=2): 


PRDAi,l,6  =  ALTXi 
PRDAilj  =  ALXi 


PRDGil,! 

=  EPNLX  I-  VFCT 

EPNdB 

PRDGil,2 

=  SELTXi  -  VFCT 

dB 

PRDGi,l,3 

=  SELXi  -  VFCT 

dB 

PRDGil,4 

=  PNLTXi 

PNdB 

PRDGi,l,5 

=  PNLXi 

PNdB 

PRDGil,6 

=  ALTXi 

dBA 

PRDGi,l,7 

=  ALXi 

dBA 

where  VFCT  adjusts  the  airspeed  from  the  reference  file  airspeed  for  the  power  setting  (IVl) 
to  the  program  reference  airspeed  (RV)  in  knots: 


( 

VFCr  =  10  log 

V 


RV  is  currently  defined  as  250  knots  in  the  main  routine.  These  single  event  measure  data, 
which  were  computed  from  the  single  event  reference  data  and  the  reference  SPL  spectrum 
extrapolated  to  each  of  the  22  profile  distances,  are  frequently  identified  as  "reference  file 
measure  data"  or  "reference  file  profile  data"  in  the  remainder  of  this  OMEGA  10 
documentation. 


SUBROUTINE  CPNL(I,ID) 

This  subroutine  is  called  from  the  MAIN  routine  and  subroutine  CDIST  to  compute  the 
perceived  noise  level  (PNL)  for  the  ph  spectrum  using  the  standard  method  described  in  several 
references.(^>^)  The  PNL  quantifies  the  relative  subjective  noisiness  of  different  sound  spectra 
and  is  widely  used  to  assess  the  annoyance  of  individual  sounds. The  ID  subroutine 
argument  determines  whether  the  PNL  data  are  computed  from  the  mean  SPL  spectrum  for  the 
power  setting  (SPLAj  j  for  ID=0)  or  the  extrapolated  SPL  spectrum  for  the  distance 
(SPLXijforID=l). 

Variables  BL  and  BBNH  are  the  indices  of  frequency  bands  17  and  40  (frequencies  50  to 
10,000  Hz)  as  required  by  the  PNL  algorithm.  FJ  is  a  PNL  weighting  factor  defined  for  the  1/3 
octave  band  data. 
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Method 

PNL  is  calculated  as  follows  from  each  1/3  octave  band  SPL  spectrum  (I): 

(1)  Convert  each  1/3  octave  band  SPL  (array  SPLA  or  SPLX)  for  frequency  band 
indices  EL  to  IBNH  to  perceived  noisiness  FNj  using  function  FNOY  which  defines 
the  noisiness  of  sound  in  noy  units  as  a  function  of  frequency  and  SPL.  Also  sum 
these  noy  values  (SUM)  and  determine  the  largest  value  (AMX). 

(2)  Determine  the  total  perceived  noisiness,  SUM,  as  follows: 

SUM  =  {FJ)[^{FNj)-AMX]  +  AMX  noys 

j 

where  FJ  =  0.15  for  1/3  octave  band  data. 

FNj  =  perceived  noisiness  value  for  frequency  index  J  from  Step  (1). 
AMX  =  number  of  noys  in  the  noisiest  band. 

(3)  Calculate  PNL: 

PNL  =  40+33.3  log  (SUM)  PNdB 

When  one  or  more  SPL  data  points  are  beyond  the  range  of  the  noy  algorithm,  PNL  is  set 
equal  to  9999. 

FUNCTION  FNOY(SPL,JJ) 

This  function  is  used  by  subroutine  CPNL  to  compute  the  perceived  noisiness  value  (in 
noys)  for  a  given  1/3  octave  band  sound  pressure  level  (SPL)  using  the  method  described  in  SAE 
ARP  865 A.(5)  Function  argument  JJ  is  the  frequency  band  index  (numeric  value  of  1  to  24) 
corresponding  to  the  24  frequency  bands  from  50  to  10,000  Hz  and  used  in  data  statement  arrays 
FL  and  FM. 

Method 

The  perceived  noisiness  value  FNOY,  in  noys,  for  a  particular  frequency  band  (JJ)  is 
related  to  the  band  sound  pressure  level,  SPL,  by  the  following  equations: 

(1)  For  FLjj^  1  <  SPL  <  FLjj  2 

FNOY  =  0.1  antilog  [(FMjj  i)(SPL  -  FLjj  j)]  noys 
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(2) 

For  FLjj_2  <  SPL  <  FLjj  3 

FNOY  =  antilog  [(FMjj  2)(SPL  -  FLjj  3)] 

noys 

(3) 

For  FLjj  3  <  SPL  <  FLjj  4 

FNOY  =  antilog  [(FMjj  3)(SPL  -  FLjj  3)] 

noys 

(4) 

ForFLjj^4<SPL<  150 

FNOY  =  antilog  [(FMjj,4)(SPL  -  FLjj^s)] 

noys 

(5) 

For  SPL  <  FLjj^  1 

FNOY  =  0.0 

noys 

(6) 

For  SPL  >  150,  FNOY  is  undefined  because  it  is  beyond  the  range  of  the  perceived 
noisiness  algorithm  (carried  as  5001). 

Data  statement  arrays  FL(24,5)  and  FM(24,4)  contain  the  Lj  to  L4  and  Mji  to  M4  versus 
frequency  data  listed  in  Table  4  (from  SAE  ARP  865A(5)). 

SUBROUTINE  CPTC(PTC,I,ID) 

This  subroutine  is  called  from  the  MAIN  routine  and  from  subroutine  CDIST  to  compute 
the  tone  correction  (PTC)  for  the  spectrum.  This  spectrum  will  be  from  array  SPLAj  j  for 
the  P^  power  setting  when  the  ID  subroutine  argument  is  zero  (CPTC  called  from  MAIN)  or 
from  array  SPLXj  j  for  the  ph  distance  for  ID=1  (CPTC  called  from  CDIST).  The  blank 
common  array  EPTCj  is  used  to  store  the  frequency  index  (J)  of  the  frequency  band  which 
determined  the  tone  correction  for  the  P^  spectrum. 

This  subroutine  uses  the  procedure  described  in  FAR  Part  36,  Section  36.3.  W  This 
procedure  requires  SPL  data  for  1/3  octave  frequency  bands  80  to  10,000  Hz;  however,  this 
subroutine  will  compute  PTC  when  SPL  data  greater  than  or  equal  to  20  dB  are  available  for  at 
least  10  consecutive  bands  within  this  1/3  octave  frequency  band  range  (only  SPL  data  greater 
than  or  equal  to  20  dB  are  used). 

Method 

Before  applying  the  tone  correction  procedure,  this  subroutine  determine  the  largest  SPL 
data  point  in  the  spectrum  and  also  the  number  of  consecutive  SPL  data  points  greater  than  20  dB 
on  each  side  of  this  peak  value  (statement  labels  370  to  490).  If  the  total  number  of  points  is  less 
than  ten,  tone  correction  is  undefined  (9999)  and  control  is  returned  to  the  calling  routine.  When 
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TABLE  4 

DATA  STORED  IN  ARRAYS  FL(24,5)  AND  FM(24,4)  IN  FUNCTION  FN0Y(5) 


Mj** 

LZ 

M2 

0.079520 

55 

0.058098 

0.068160 

51 

0.058098 

0.068160 

46 

0.052288 

0.059640 

42 

0.047534 

0.053013 

39 

0.043573 

0.053013 

m 

0.043573 

0.053013 

33 

0.040221 

0.053013 

30 

0.037349 

0.053013 

27 

0.034859 

0.053013 

25 

0.034859 

0.053013 

25 

0.034859 

0.053013 

25 

0.034859 

0.053013  i 

25 

0.034859 

25 

0.034859 

0.059640 

23 

0.034859 

0.053013 

21 

0.040221 

0.053013 

18 

0.037349 

0.047712 

15 

0.034859 

0.047712 

14 

0.034859 

0.053013 

14 

0.034859 

0.053013 

15 

0.034859 

0.068160 

17 

0.037349 

IHSHI 

23 

0.037349 

0.0596401 

29 

0.043573 

0.043478 


0.040570 


0.036831 


0.036831 


0.035336 


0.033333 


0.033333 


0.032051 


0.030675 


0.030103 


0.030103 


0.030103 


0.030103 


0.030103 


0.029960 


0.029960 


0.029960 


0.029960 


0.029960 


0.029960 


0.029960 


0.042285 


0.042285 


85.88 


87.32 


79.85 


79.76 


75.96 


73.96 


74.91 


94.63 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


44.29 


50.72 


M4 

L4 

0.030103 

52 

0.030103 

51 

0.030103 

49 

0.030103 

47 

0.030103 

46 

0.030103 

45 

0.030103 

43 

0.030103 

42 

0.030103 

41 

0.030103 

40 

0.030103 

40 

•  0.030103 

40 

0.030103 

40 

0.030103 

40 

0.030103 

38 

0.029960 

34 

0.029960 

32 

0.029960 

30 

0.029960 

29 

0.029960 

29 

0.029960 

30 

0.029960 

31 

0.029960 

34 

0.029960 

37 

*Lj  to  L4  data  are  stored  in  FL(24,5). 
**Mi  to  M4  data  are  stored  in  FM(24,4). 
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at  least  ten  points  are  available,  the  tone  correction  procedure  is  applied  over  redefined  frequency 
band  indices  IL2  to  IHl. 

The  remainder  of  this  subroutine  (statement  labels  5  to  220)  is  a  direct  application  with 
one  exception  of  the  ten  step  tone  correction  procedure  defined  in  FAR  Part  36,  Section  36.3. 
This  exception  to  Far  Part  36  is  explained  in  item  (4)  below.  The  final  tone  correction  is  returned 
in  variable  PTC.  The  following  comments  compare  the  notation  used  in  the  coding  and  in  the 
FAR  Part  36  description: 

(1)  The  frequency  band  index  i  in  Part  36  corresponds  to  J  in  the  coding.  Spectrum 
index  K  doesn't  apply  because  the  subroutine  operates  on  one  spectrum  per  call. 

(2)  Circled  SPL  data  are  denoted  by  variable  ICTj  equal  to  one. 

(3)  All  other  variable  correspondence  should  be  obvious. 

(4)  The  direct  application  of  Far  Part  36,  Section  36.3  to  compute  tone  correction  will 
occasionally  result  in  a  false  tone  at  a  frequency  greater  than  2500  Hz.  The  if 
statement  at  label  1 10  was  added  to  check  two  SPL  slopes  at  the  Fj  frequency.  This 
check  eliminates  most  of  the  false  tones  at  high  frequencies. 

SUBROUTINE  CAL(I,ID) 

This  subroutine  is  called  from  the  MAIN  routine  and  from  subroutine  CDIST  to  compute 
the  A-weighting  overall  sound  level  (AL)  for  the  spectrum.  The  subroutine  argument  ID 
determines  whether  the  AL  data  are  computed  from  the  mean  SPL  spectrum  for  the  power 
setting  or  the  extrapolated  SPL  spectrum  for  the  F^  distance.  ID  is  defined  the  same  as  in 
subroutines  CPTC  and  CPNL. 


Method 

AL  is  defined  as  follows: 


where 

J 

SPLij 


AWj 


AL,  =  101og[^antilog^ 


SPL,j+AWj^ 
10 


] 


dBA 


the  frequency  band  index;  J  is  defined  for  bands  10  to  40  (indices  1  to  31). 

the  sound  pressure  level  for  the  F^  spectrum  or  F^  distance  and  the  band 
(arrays  SPLX  or  SPLA). 

the  A-weighting  relative  response  in  dB  for  the  band  (see  Table  5). 
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TABLES 

WEIGHTING  FACTORS 
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If  the  AL  data  are  undefined  for  the  spectrum,  AL  is  set  equal  to  9999;  this  should  never 

occur  in  this  program. 

SUBROUTINE  OUTH(IRD,IPTC,SENX,LFLG,IPF) 

This  subroutine  is  called  from  the  MAIN  routine  to  print  the  SPL  spectra  versus  slant 
distance  on  Pages  H  and  L  and/or  the  single  event  noise  data  (seven  measures)  versus  slant 
distances  on  Pages  I  and  M.  The  air-to-ground  data  (IPROP=l)  are  printed  on  Pages  H  and  I  and 
the  ground-to-ground  data  (IPROP=2)  on  Pages  L  and  M.  Subroutine  argument  IRD  is  the  index 
of  the  profile  distance  which  is  within  one  percent  of  the  reference  distance  for  this  power 
setting.  The  remaining  arguments  are  defined  in  subroutine  comment  statements.  The  data 
format  and  page  layout  for  each  page  can  best  be  determined  by  consulting  the  sample  problem 
in  Appendix  C. 

Print  SPL  Spectra 

The  first  segment  of  this  subroutine  which  is  omitted  when  IPF>2  prints  the  SPL  spectra 
versus  distance  data  for  the  22  profile  distances.  Subroutine  HEAD  is  called  to  print  the  page 
header  block.  The  statement  label  50  loop  sets  up  the  integer  and  variable  format  arrays  and 
prints  each  of  the  22  spectra.  All  SPL  data  less  than  zero  are  omitted  from  the  printout  (blanked 
out).  The  frequency  band  w'hich  determines  the  tone  correction  in  the  reference  distance 
spectrum  is  flagged  with  the  ">"  symbol.  The  variable  format  is  required  to  accommodate  the 
blank  and  flagged  data. 

Print  Single  Event  Data 

The  remainder  of  this  subroutine  (from  label  50)  prints  the  single  event  versus  distance 
data  for  the  22  profile  distances.  The  sequence  of  these  seven  single  event  measures  in  array 
SENX  is  determined  by  the  EPF  program  flag.  For  IPF<2,  the  sequence  is  PNLTX,  PNLX, 
ALTX,  ALX,  EPNLX,  SELTX  and  SELX;  while  for  IPF>2,  the  sequence  is  EPNLX,  SELTX, 
SELX,  PNLTX,  PNLX,  ALTX  and  ALX.  Neither  sequence  is  the  same  as  the  print  sequence 
which  is  just  the  reverse  of  the  latter  case  (IPF>2).  Thus,  index  variable  JD  is  setup  to  identify 
the  single  event  data  in  array  SENX  in  the  required  print  sequence.  For  the  IPF>2  case,  variable 
JD  is  defined  by  the  ISEQD  data  statement  array  which  contains  the  indices  of  these  single  event 
measures  in  the  print  sequence.  These  data  are  first  written  to  the  DM  character  array  to  blank 
out  all  negative  single  event  data.  This  character  array  is  then  printed  with  format  statement 
2700. 
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At  label  55  subroutine  HEAD  is  called  to  print  the  Page  I  or  M  header  block.  The  label 
100  loop  prints  the  measure  data  for  each  distance.  The  label  80  loop  sets  up  the  vanable  print 
array  for  each  measure  for  the  distance. 

SUBROUTINE  OUTJ(IP,SENX,LFLG)  \ 

This  subroutine  is  called  from  the  MAIN  routine  to  print  the  single  event  measure  printer 
plots  for  air-to-ground  (IPROP=l,  Pages  J  and  K)  or  ground-to-ground  (IPROP=2,  Pages  N  and 
O)  data  for  the  Uh  power  setting.  The  value  of  the  subroutine  argument  IP  determines  the  single 
event  measures  plotted  on  each  page.  When  IP  is  one  or  four,  the  PNLTX,  PNLX,  ALTX,  and 
ALX  data  are  plotted  on  Pages  J  or  N,  and  when  IP  is  two  or  three,  the  EPNLX,  SELTX,  and 
SELX  data  are  plotted  on  Pages  K  and  O.  This  IP  program  flag  also  determines  the  sequence  in 
which  these  single  event  data  are  stored  in  array  SENX.  This  data  sequence  is  described  in  detail 
in  subroutine  OUTH  where  program  flag  EPF  is  equivalent  to  BP  in  this  routine.  Subroutine 
argument  LFLG  is  a  footnote  flag;  for  LFLG=1,  an  extrapolation  limit  flag  is  pnnted  at  the 
bottom  of  each  page.  The  plot  format  and  page  layout  can  best  be  determined  by  consulting  the 

sample  problem  in  Appendix  C. 

Method 

Each  printer  plot  contains  the  distance  scale  on  the  ordinate  and  the  dB  scale  on  the 
abscissa.  The  distance  scale  always  contains  the  same  21  profile  distances  from  200  to  25,000 
feet.  The  abscissa  scale  is  8.3  inches  wide  and  covers  an  83  dB  range.  The  maximum  abscissa 
scale  value  is  determined  from  the  measure  data  for  the  first  ten  distances  (smallest  distances). 
The  maximum  annotated  scale  value  is  the  largest  multiple  of  ten  which  is  less  than  the  largest 
measure  plus  eight.  The  largest  plotted  point  may  be  one  larger  than  this  multiple  of  ten.  All 
data  outside  the  abscissa  scale  range  are  deleted  from  the  plot. 

The  label  400  loop  sets  up  the  plot  array  and  the  grid  pattern  for  each  profile  distance. 

An  identification  block  is  printed  in  the  upper  left  comer  and  a  legend  block  is  pnnted  m  the 
lower  right  comer  of  the  grid.  These  blocks  do  not  normally  interfere  with  the  plotted  data 
points;  however,  they  do  make  the  coding  more  complex  and  thus  more  difficult  to  change. 

The  label  150  loop  scales  the  three  or  four  (J1  to  J2)  data  points  plotted  for  each  distance, 
sets  up  the  scaled  point  in  the  plot  array  with  the  correct  symbol,  and  saves  the  data  presently  in 
the  plot  array.  There  should  be  sufficient  comment  statements  in  this  subroutine  to  follow  the 
data  being  printed;  thus  no  additional  documentation  will  be  given  here. 
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SUBROUTINE  PPFDAT(PRDC,LFLG,COMD) 

This  subroutine  is  called  from  the  MAIN  routine  to  write  the  profile  datasets  on  Unit  3 
for  the  output  power  setting  (PSC).  The  air-to-ground  and  ground-to-ground  data  are  in 
array  PRDC.  LFLG  is  a  program  flag  which  when  equal  to  one  adds  an  additional  extrapolation 
limit  comment  line  to  each  profile  dataset.  The  format  and  content  of  these  profile  datasets  is 
given  in  Appendix  F.  There  should  be  sufficient  comment  statements  in  this  subroutine  to  follow 

the  data  being  written. 

SUBROUTINE  TITPG(IPR) 

This  subroutine  is  called  from  the  MAIN  routine  to  print  the  aircraft  title  page  which 
identifies  the  aircraft  being  analyzed  and  the  data  computed  in  the  analysis.  The  format  and 
content  of  this  page  can  best  be  determined  by  consulting  the  sample  problem  in  Appendix  C. 
There  should  be  sufficient  comment  statements  in  the  coding  to  follow  the  information  being 
printed.  IPR  is  the  program  print  flag. 

SUBROUTINE  SETUPD6(IREQ,N,NP,TYPE) 

This  subroutine  is  called  from  the  MAIN  routine  to  apply  the  interpolation/extrapolation 
rules  to  determine  the  reference  data  required  to  compute  the  profile  data  at  the  requested  output 
power  setting  (PSC).  The  indices  of  these  reference  data  required  for  each  PSC  are  stored  in 
array  IREQC  and  used  in  subroutine  DELTA6.  The  IREQ,  N,  and  NP  subroutine  arguments  and 
numerous  additional  variables  are  defined  in  the  comment  statements  at  the  beginning  of  this 
subroutine.  Argument  TYPE  is  a  character  array  used  to  store  the  interpolation  code  (I,  N  or  P) 
for  each  dataset  reference. 

In  the  OMEGA  10  version  7  update  in  May  1991,  the  method  used  to  interpolate  noise 
measure  data  at  non  reference  power  conditions  was  changed.  Formerly  a  complex  set  of  rules 
based  on  the  operation  power  code  (OPC)  was  used  to  setup  the  interpolation  equation.  These 
rules  were  abandoned  and  each  reference  dataset  was  assigned  an  interpolation  code  which  is 
used  to  setup  the  power  setting  interpolation  equation.  These  interpolation  codes  were  assigned 
to  each  reference  dataset  after  reviewing  all  the  available  data  for  a  specific  aircraft.  The  code  is 
stored  as  the  last  character  in  the  COMDECK  name  and  replaces  the  COMDECK  revision  code 
in  the  flight  reference  data  COMDECK  name  in  NOISEFILE. 

The  three  interpolation  codes  currently  used  in  NOISEFILE  and  OMEGA  10  are  I,  P  and 
N.  Datasets  coded  as  I  for  a  specific  aircraft  are  used  to  define  a  power  setting  versus  noise  level 
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point  to  point  linear  interpolation  function  which  is  used  to  interpolate  noise  level  for  an  output 
power  setting  whose  reference  OPCR  is  coded  as  I  and  whose  power  setting  is  not  the  same  as 
any  of  the  reference  file  power  settings.  The  P  interpolation  code  stands  for  parallel  inteipolation 
and  the  N  interpolation  code  stands  for  no  interpolation.  Parallel  interpolation  means  that  the 
coded  as  P  are  interpolated  along  a  slope  line  determined  by  the  point  to  point  I  interpolation 
function  but  with  the  slope  line  offset  so  it  passes  through  the  reference  power  versus  noise  level 
point  (OPCR).  Output  power  settings  requested  using  an  OPCR  with  type  code  N  must  match 
the  reference  power  setting  since  no  power  interpolation  is  permitted. 

Method 

The  OMEGA  10  setup  file  (code  sheet  input)  may  contain  up  to  twelve  output  power 
settings  (PSC)  for  each  aircraft  set.  Each  of  these  PSCs  must  have  a  reference  operation  power 
code  (OPCR)  identified  in  the  setup  file.  This  OPCR  must  be  in  the  reference  dataset  file  or  data 
for  the  PSC  are  omitted  from  the  computer  job.  This  OPCR  determines  the  operation  power 
description  assigned  to  the  corresponding  PSC  as  well  as  the  type  of  interpolation  (I,  N  or  P) 
used  to  compute  the  profile  data.  A  brief  overview  of  the  interpolation  procedure  is  described  as 
follows: 

(1)  If  the  OPCR  identifies  the  PSC  as  type  I,  then  this  subroutine  identifies  the  two  type 
I  reference  datasets  with  power  settings  which  bracket  the  given  PSC.  The  PSC  may 
be  interpolated  anywhere  along  the  type  I  power  versus  noise  level  function.  PSC 
extrapolations  beyond  the  high  or  low  end  points  use  the  last  two  type  I  reference 
data  points  and  are  limited  to  5  dB  beyond  the  end  point. 

(2)  If  the  OPCR  identifies  the  PSC  as  type  P,  then  this  subroutine  identifies  the  two  type 
I  reference  datasets  with  power  settings  which  bracket  the  reference  power  setting 
identified  by  the  OPCR.  This  pair  of  type  I  reference  points  are  used  to  determine 
the  slope  line  which  is  offset  to  pass  through  the  OPCR  noise  versus  power  setting 
point.  If  the  given  PSC  also  falls  within  the  power  setting  range  of  these  two  type  I 
reference  points,  the  noise  level  for  this  PSC  is  interpolated  along  the  slope  line 
which  is  offset  to  pass  through  the  OPCR  point.  If  the  PSC  power  setting  is  outside 
the  reference  points  which  bracket  the  OPCR,  then  the  noise  level  offset  between  the 
OPCR  and  this  slope  line  is  computed  and  applied  to  a  second  slope  line  determined 
by  a  pair  of  type  I  reference  points  which  bracket  the  given  PSC  power  setting.  This 
pair  of  type  I  reference  points  which  bracket  the  given  PSC  may  be  anywhere  along 
the  reference  power  versus  noise  level  function;  however,  the  extrapolation  from  the 
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noise  level  at  the  reference  power  setting  (OPCR)  and  at  the  reference  disf^ce  is 
limited  to  5  dB. 

(3)  If  the  OPCR  identifies  the  PSC  as  type  N,  then  this  PSC  must  exactly  match  the 
power  setting  (PS)  for  the  OPC  identified  by  this  OPCR.  If  PS  is  not  the  same  as 
PSC,  this  PSC  is  omitted  from  the  OMEGA  10  job  and  an  error  message  is  printed. 

This  subroutine  identifies  the  reference  file  OPC  power  settings  required  to  interpolate 
the  output  noise  profiles  for  each  output  PSC.  The  interpolation  is  done  later  in  subroutine 
DELTA6.  Array  IREQ(15)  is  used  as  a  flag  to  identify  the  reference  datasets  which  are  used  to 
compute  at  least  one  PSC.  If  IREQj  is  greater  than  zero,  this  reference  file  OPC  is  required 
and  all  reference  profile  data  must  be  computed.  Array  IREQC(7,15)  is  used  to  store  the  indices 
of  the  reference  OPC's  required  to  compute  the  profile  data  for  each  output  PSC.  A  brief 
description  of  this  IREQC  array  is  as  follows. 

(1)  If  IREQCi  j  is  zero  for  I  equal  one  to  three,  then  no  profile  data  are  computed  for 
this  Jth  PSC. 

(2)  If  IREQC3  j  is  zero  and  IREQCi  j  is  greater  than  zero,  then  no  interpolation  is 
required  and  IREQCi  j  contains  the  reference  file  index  of  the  required  profile  data. 

(3)  For  a  typical  interpolation  where  there  are  two  slope  end  points  and  the  reference 
point  from  which  the  interpolation  is  computed,  the  indices  of  the  two  slope  points 
are  stored  in  IREQCi  j  IREQC2,J  the  reference  point  is  stored  in 
IREQC3,j.  This  IREQC3  j  index  will  be  negative  when  the  5  dB  interpolation  limit 
must  be  checked  at  the  reference  distance.  For  type  I  interpolation,  index  IREQC3  j 
win  usually  be  the  same  as  either  the  IREQCi  j  IREQC2,J  index.  For  type  P 
interpolation,  the  three  indices  will  usually  be  different  since  the  type  P  reference 
point  is  usually  offset  from  the  type  I  slope  line. 

(4)  For  a  type  P  interpolation  which  requires  two  slope  lines,  the  first  slope  line  is  used 
to  determine  the  noise  level  offset  between  the  OPCR  point  and  the  type  I 
interpolation  slope  line  at  the  OPCR  power  setting.  The  second  slope  line  end 
points  bracket  the  output  PSC  for  this  OPCR;  this  slope  line  is  used  to  interpolate 
the  noise  levels  at  this  PSC  power  setting.  The  offset  from  the  first  slope  line  is 
added  to  the  interpolated  levels.  In  this  subroutine  the  indices  of  the  first  slope  line 
end  points  are  stored  in  IREQC4  j  and  IREQC5  j.  The  index  of  the  OPCR 
reference  point  is  in  IREQC5  j.  The  corresponding  indices  for  the  second  slope  line 
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are  saved  in  IREQC^  j  through  IREQC3  j.  The  IREQC3  j  and  IREQC5  j  reference 
point  indices  are  multiplied  by  minus  one  when  the  extrapolation  limit  must  be 
checked  in  subroutine  DELTA6.  The  IREQC7  j  or  seventh  array  element  is  used  to 
store  the  index  of  a  type  I  reference  power  setting  which  is  between  the  two  pairs  of 
slope  end  points.  This  fifth  type  I  point  is  not  used  to  compute  either  slope  line;  it 
is  used  in  the  interpolation  procedure  in  subroutine  DELTA6.  (Note  that  as  of 
February  10, 1993  no  NOISEFILE  reference  data  has  five  type  I  points  with  one  or 
more  type  P  points.  Thus  IREQC7  j  will  always  be  zero.) 

In  this  subroutine  the  loops  labeled  15  and  20  use  internal  read  statements  to  convert 
reference  (PS)  and  output  (PSC)  power  setting  values  from  character  to  real.  The  real  arrays  are 
PSIF  and  PSCF,  respectively. 

From  label  30  to  label  90  the  indices  of  the  KK  type  N  and  type  P  reference  data  are 
stored  in  the  first  KK  elements  of  array  NRI.  The  type  I  data  are  then  ranked  by  increasing 
power  setting  starting  with  rank  KK+1.  The  index  of  the  type  I  data  of  rank  K  is  stored  in  the 
K^  element  of  array  NRI. 

The  label  400  loop  sets  up  arrays  IREQ  and  IREQC  for  each  of  the  NP  output  power 
settings  (PSCF).  The  indices  setup  in  IREQC  are  used  in  subroutine  DELTA6  to  interpolate  the 
output  profile  data.  The  label  95  loop  identifies  the  index  (J)  of  the  reference  operation  power 
code  (OPC)  which  matches  the  operation  power  code  (OPCR)  assigned  to  the  output  power 
setting.  Variable  "J"  will  be  used  to  store  this  index  throughout  the  label  400  loop.  The  next 
section  of  code  from  label  100  to  label  105  checks  for  type  N  and  type  P  interpolation  and  also 
checks  the  input  (PSIF)  and  output  (PSCF)  power  settings  to  determine  if  they  are  the  same  in 
which  case  no  interpolation  is  required. 

The  label  110  loop  determines  the  type  I  power  setting  indices  which  bracket  the  output 
power  setting  for  type  I  or  bracket  the  input  power  setting  for  type  P  data.  The  first  three 
elements  of  IREQC  are  set  in  this  loop.  The  next  section  of  code  starting  at  label  115  checks  the 
type  P  data  to  determine  if  a  second  slope  line  in  required.  When  a  second  slope  line  is  required 
for  type  P,  elements  four  to  six  of  IREQC  are  set  equal  to  the  indices  which  bracket  the  reference 
PSIF  power  setting  and  a  second  pass  through  the  label  1 10  loop  defines  the  bracketing  indices 
around  the  output  PSCF  power  setting. 

The  section  of  code  containing  the  label  220  loop  is  executed  only  for  type  P 
interpolation  where  two  slope  lines  are  required.  The  type  I  power  setting  points  are  checked  to 
determine  if  there  is  an  unused  type  I  power  setting  value  between  the  two  slope  lines.  If  this  is 
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the  case,  the  index  of  this  point  is  saved  in  element  seven  of  IREQC.  This  index  will  be  needed 
in  subroutine  DELTA6  only  if  the  extrapolation  limit  is  exceeded  using  the  second  slope  line  for 
this  type  P  interpolation. 

The  section  of  code  from  labels  370  to  395  prints  error  messages  to  the  list  file  and  sets 
ERRFLG  to  TRUE.  The  label  500  to  510  code  is  executed  when  the  NP  equal  zero  option  is 
selected.  For  NP  equal  zero,  profile  data  are  computed  for  all  reference  file  power  settings  at  the 
reference  file  conditions  for  the  selected  aircraft.  The  first  N  elements  of  IREQ  are  set  to  one  and 
the  first  element  of  EREQC  for  each  of  the  N  power  settings  is  set  equal  to  the  index  of  that 
power  setting. 

SUBROUTINE  DELTA6(PRDI,PRDC,K,PSIF,PSCF, IREQC, LFLG,VFCT,LIM,EXTMX) 


This  subroutine  is  called  from  the  MAIN  routine  to  apply  the  interpolation  algorithms  in 
computing  the  final  flight  noise  PNLX,  PNLTX,  ALX,  ALTX,  SELX,  SELTX  and  EPNLX 
profile  data  at  the  requested  power  setting  (PSC)  and  airspeed.  All  the  subroutine  arguments  are 
defined  in  the  comment  section  at  the  beginning  of  this  subroutine  (source  listing).  The  term 
"reference  file  measure  data"  used  in  the  following  paragraphs  is  the  measure  data  (PNLX,  etc.) 
derived  from  the  SPL  reference  spectrum  after  it  is  extrapolated  to  each  of  the  22  profile 
distances  (see  subroutine  CDIST). 


Method 

The  following  linear  interpolation  (or  extrapolation)  function  (F)  is  defined  in  this 
subroutine  and  used  to  compute  the  single  event  noise  measures  at  each  profile  distance: 


where 


yc= 


y2-yi 

X2-XI 


j(XC-XB)  +  y5 


YC 

Y1  and  Y2 
YB 


XI  and  X2 

xc 

XB 


the  interpolated  or  extrapolated  measure  data  for  the  PSC  power  setting 
(array  PRDC); 

the  reference  file  measure  data  used  to  compute  the  slope  (array  PRDI); 

the  reference  file  measure  data  reference  point  defined  by  OPCR  (array 
PRDI); 

the  reference  file  power  setting  data  used  to  compute  the  slope  (array 
PSBF); 

the  requested  power  setting  (array  PSCF);  and 

the  reference  file  power  setting  defined  by  OPCR  (array  PSBF). 
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The  airspeed  adjustment  factor  (VFCT)  is  subtracted  from  this  YC  data  where 

VFCr  =  10  logf^ 

\  RV 

VXl  =  profile  output  airspeed  (knots)  for  PSC. 

RV  =  program  reference  airspeed  (presently  250  knots). 

Note  that  VFCT  is  always  zero  for  all  PNLX,  PNLTX,  ALX  and  ALTX  data  since  the  airspeed 
adjustment  only  applies  to  the  integrated  measures. 

The  linear  interpolation  algorithm  used  to  compute  the  limiting  power  setting  when  the 
noise  level  extrapolation  limit  (EXTMX)  is  exceeded  is: 

XC  =  XB  +  -  XI) 

where 

XC  =  the  new  PSC  power  setting  (PSCF)  corresponding  to  the  limiting 

measure  value  (YL);  and 

YL  =  the  limiting  measure  value  (YL=YB±EXTMX),  YL  is  EXT  in  the 

program. 

All  other  variables  are  defined  for  the  YC  equation  above. 

Array  EREQC  is  defined  in  detail  in  the  comments  at  the  beginning  of  this  subroutine  as 
well  as  in  the  subroutine  SETUPD6  documentation.  The  interpolation  procedure  is  also 
described  in  detail  in  subroutine  SETUPD6.  With  this  background  there  should  be  sufficient 
comment  statements  in  the  source  listing  to  follow  the  coding;  however,  the  following  outline  of 
this  subroutine  may  be  helpful.  The  principle  segments  of  this  subroutine  are: 

(1)  The  statement  label  40  loop  uses  the  slope  and  reference  points  to  interpolate  the 
measure  data  and  apply  the  airspeed  adjustment.  This  is  a  direct  application  of  the 
above  YC  equation. 

(2)  The  label  50  loop  is  used  to  compute  the  type  P  interpolation  where  two  slope  lines 
are  required.  First  the  noise  measure  offset  is  computed  at  the  reference  PSIF^g 
power  setting.  Next  the  interpolated  noise  measure  data  are  computed  at  the  output 
PSCF  power  setting. 
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(3)  The  label  80  loop  applies  the  airspeed  adjustment  to  the  given  reference  file  noise 
measure  data.  For  these  data  the  reference  and  output  power  settings  are  the  same 
and  no  interpolation  is  required. 

(4)  The  segment  of  code  in  the  IF  statement  following  label  150  checks  the  extrapolation 
limit  for  type  P  interpolation  where  two  slopes  are  required.  This  limit  is  checked 
only  at  the  1000  foot  reference  distance  and  only  for  one  noise  measure.  If  the 
computed  noise  measure  exceeds  the  extrapolation  limit  at  the  low  end,  program 
execution  is  transferred  to  label  180  where  a  new  output  power  setting  is  determined 
for  this  minimum  noise  measure  level.  If  the  limit  is  exceeded  at  the  high  end, 
program  execution  is  transferred  to  label  182  to  compute  the  new  power  setting  at 
the  upper  limit. 

(5)  The  next  section  of  code  starting  at  label  153  checks  the  extrapolation  limits  for  the 
single  slope  type  I  or  type  P  extrapolations.  If  a  limit  is  exceeded,  a  new  output 
PSCF  power  setting  is  computed  at  label  170. 

(6)  The  segments  of  code  following  labels  180  and  182  compute  new  limiting  output 
power  settings  for  cases  where  the  computed  noise  measure  was  below  the  lower 
limit  or  above  the  upper  limit,  respectively.  In  each  segment  the  new  PSCF  is 
checked  and  the  appropriate  slope  indices  are  assigned.  The  PSCF  is  recomputed  as 
needed  when  the  slope  line  changes. 

(7)  The  label  1 85  segment  is  executed  when  the  PSCF  computed  in  one  of  the  previous 
sections  (labels  180  and  182)  falls  within  the  first  slope  PSIF  range.  This  eliminated 
the  need  for  the  second  slope  line.  The  PSCF  is  recomputed  using  the  first  slope 
reference  points  and  the  BREQC  indices  are  redefined. 

(8)  The  final  segment  of  code  starting  at  label  180  sets  the  extrapolation  limit  flag  LFLG 
to  indicate  that  the  output  power  setting  was  recomputed.  This  code  also  formats  the 
computed  PSCF  into  the  PSC  character  array.  Program  execution  is  returned  to  label 
25  where  the  profile  data  are  computed. 

Note  that  all  code  after  label  150  deals  with  checking  the  extrapolation  limits  and 
recomputing  a  new  PSCF  where  the  limit  is  exceeded. 
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SUBROUTINE  SUMRY(IPU,COMD,EXTMX,N,NP, SOURCE, LFLG) 


This  subroutine  is  called  from  the  MAIN  routine  at  the  end  of  each  aircraft  analysis  to 
print  a  sunmiary  of  the  input  and  output  data  for  that  aircraft.  The  subroutine  arguments  are; 

(1)  IPU  is  the  profile  dataset  print  (on  Unit  3)  flag; 

(2)  COMD  contains  the  last  5  characters  of  the  reference  dataset  comdeck  name; 

(3)  EXTMX  is  the  maximum  extrapolation  permitted  by  the  program  (see  subroutine 
DELTA6); 

(4)  N  is  the  number  of  input  power  settings  (PS); 

(5)  NP  is  the  number  of  output  power  settings  (PSC); 

(6)  SOURCE  contains  the  date  of  the  original  reference  dataset  run;  and 

(7)  LFLG  is  a  flag  defined  for  each  output  power  setting  (see  subroutine  DELTA6). 

If  the  sample  problem  in  Appendix  C  is  used  as  a  guide,  there  should  be  sufficient 
comment  statements  in  this  subroutine  to  follow  the  data  being  printed.  The  summary  of  the 
output  data  includes  the  reference  file  comdeck  names  for  all  reference  data  used  to  interpolate 
the  profile  data  at  the  output  power  setting.  Information  in  arrays  IREQC  is  used  to  print  these 
comdeck  names. 

SUBROUTINE  RDFILE(NFFLYI,NFFLY2,*) 

This  subroutine  is  called  from  the  MAIN  routine  when  the  command  line  variable  ARCS 
in  the  MAIN  routines  is  NULL.  The  subroutine  arguments  NFFLYl  and  NFFLY2  are  two 
twenty-four  bjte  character  variables  used  to  store  the  NOISEFILE  database  filenames  for 
military  and  civil  aircraft,  respectively.  When  the  command  line  is  NULL,  these  two  filenames 
must  be  defined  in  the  MAIN  routine.  The  asterisk  in  the  argument  list  is  an  alternate  RETURN 
label  which  is  used  to  return  control  to  the  calling  routine  when  an  error  has  occurred. 

This  subroutine  performs  the  following  functions: 

(1)  Unit  1  is  opened  for  I/O  to  the  display; 

(2)  The  program  queries  the  operator  for  the  OMEGA  10  setup  filename  which  must  be 
a  standard  DOS  filename  with  a  maximum  of  twelve  characters.  This  file  is  opened 
on  Unit  5  by  the  program.  The  ERRFLG  is  set  to  TRUE  if  the  operator  does  not 
enter  an  exiting  filename. 
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(3)  The  program  queries  the  operator  for  the  profile  dataset  output  filename  which  must 
be  a  standard  DOS  filename  with  a  maximum  of  twelve  characters.  This  file  is 
opened  on  Unit  3  by  program. 

(4)  The  NOISEFILE  database  file  with  names  stored  in  NFFLYl  and  NFFLY2  are 
checked  to  be  sure  they  exist.  If  they  don't  exist  the  ERRFLG  is  set  to  TRUE  and  the 
job  is  terminated.  These  files  are  not  opened  at  this  time. 

(5)  The  list  file  is  always  opened  on  Unit  6  with  filename  OMEGA108.OUT.  If  this  file 
already  exists,  the  operator  is  queried  for  a  new  filename  and  given  the  option  to 
terminate  the  Job.  If  the  Job  is  terminated,  the  ERRFLG  is  set  to  TRUE. 

SUBROUTINE  CLFILE(NFFLY1,NFFLY2^GS*) 

This  subroutine  is  called  when  the  command  line  character  variable  ARGS  contains 
filename  information  which  must  be  decoded  into  the  I/O  filenames.  This  routine  separates  the 
command  line  string  which  may  be  up  to  127  character  long  into  the  I/O  filenames.  This  routine 
will  only  run  when  compiled  with  the  Lahey  compiler;  it  will  not  run  with  Microsoft  FORTRAN 
which  was  formerly  used  with  the  OMEGA  10  program. 

The  subroutine  arguments  NFFLYl  and  NFFLY2  are  two  twenty-four  byte  character 
variables  used  to  store  the  NOISEFILE  database  filenames  for  Military  and  civil  aircraft  data, 
respectively. 

The  command  line  must  contain  five  I/O  filenames  in  the  sequence  listed  here.  Each  file 
is  described  briefly  as  follows: 

(1)  The  setup  filename  is  decoded  and  opened  on  Unit  5; 

(2)  The  list  filename  is  decoded  and  opened  on  Unit  6; 

(3)  The  profile  dataset  filename  is  decoded  and  opened  on  Unit  3; 

(4)  The  first  NOISEFILE  database  filename  is  decoded  and  stored  in  variable  NFFLYl; 
and 

(5)  The  second  NOISEFILE  database  filename  is  decoded  and  stored  in  variable 
NFFLY2. 

Each  filename  must  be  separated  by  a  space  on  the  command  line.  The  first  byte  on  the 
command  line  must  not  be  a  space.  If  an  error  occurs  when  opening  or  checking  a  file,  the 
ERRFLG  is  set  to  TRUE  and  the  Job  is  terminated.  If  only  military  or  only  civil  data  are 
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available,  the  two  NOISEFILE  filenames  may  be  the  same;  however,  two  database  filenames 
must  be  entered  on  the  command  line.  The  two  database  files  are  never  both  open  at  the  same 
time. 

BLOCK  DATA 

The  BLOCK  DATA  routine  is  the  last  routine  in  this  program.  It  is  used  to  define 
numerous  variables  used  in  labeled  common  throughout  the  program.  See  the  program  listing  in 
Appendix  I  and  the  list  of  symbols  in  this  program  documentation  to  identify  the  variables 
defined  in  BLOCK  DATA. 
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GENERAL  OVERVIEW  OF  THE  OMEGA  11  PROGRAM 


The  OMEGA  1 1  program,  hereafter  referred  to  as  simply  the  "program",  is  designed  to 
compute  descriptions  of  the  ground  run-up  noise  of  an  aircraft  in  terms  of  tone-corrected 
perceived  noise  level  (PNLTX),  A-weighted  overall  sound  level  (ALX),  and  tone-corrected  A- 
weighted  overall  sound  level  (ALTX)  as  a  function  of  distance  to  the  aircraft,  aircraft  power 
setting  and  meteorological  conditions.  These  noise  measure  data  (profile  datasets)  are  computed 
for  aircraft  ground  run-up  reference  data  as  outlined  in  AMRL-TR-73-107.(0  They  are  required 
as  input  to  the  NOISEMAP  noise  exposure  forecast  program. 

To  compute  the  above  noise  measure  data  for  selected  operation  power  settings,  the 
program  inputs  all  reference  datasets  for  aircraft  ACC  from  the  NOISEFILE  database.  These 
reference  datasets  contain  sound  pressure  level  (SPL)  data  for  19  angles  normalized  to  a 
reference  distance  of  250  feet  and  to  standard  day  temperature  (59°F),  relative  humidity  (70%) 
and  barometric  pressure  (29.92  inches  Hg).  The  format  of  these  data  are  described  in  Appendix 
G. 


These  reference  spectra  are  then  extrapolated  to  each  of  the  22  standard  profile  distances 
(from  200  to  25,0(K3-feet)  at  the  requested  standard  or  non-standard  weather  conditions.  PNLTX, 
ALX  and  ALTX  noise  measures  are  determined  for  each  distance  and  angle  from  these 
extrapolated  SPL  spectra  for  the  reference  power  setting.  This  program  interpolates  between  two 
reference  file  power  settings  to  compute  the  noise  measure  data  for  each  power  setting  requested 
in  the  setup  file.  Noise  measure  data  can  be  computed  for  any  power  setting  within  the  range 
available  in  the  reference  file. 

A  brief  summary  of  the  program  operation  is  given  below: 

(1)  The  program  inputs  the  code  sheet  parameters  described  in  Appendix  B  and  reads  the 
reference  datasets  for  aircraft  ACC  from  the  reference  file. 

(2)  The  power  setting  data  are  ranked  and  the  indices  of  the  reference  datasets  required 
to  interpolate  each  output  power  setting  are  defined. 

(3)  The  cover  page  is  printed  when  DPR  is  greater  than  zero. 

(4)  The  aircraft  summary  page  is  always  printed. 

(5)  If  Delta  N  (DELN)  is  not  zero,  it  is  added  to  all  spectra  for  all  power  settings. 

(6)  If  print  flag  IPR  is  greater  than  zero,  all  reference  datasets  are  printed. 
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(7)  PNLTX,  ALX  and  ALTX  profile  data  for  each  requested  power  setting  (PSC)  are 
computed  in  sequence  from  low  to  high  PSC: 

(a)  If  the  PSC  is  the  same  as  one  of  the  reference  dataset  power  settings,  the  profile 
data  are  computed  and  no  interpolation  is  required. 

(b)  If  the  PSC  is  not  the  same  as  a  reference  dataset  power  setting,  the  profile  data 
are  computed  at  the  reference  conditions  for  the  nearest  reference  file  power 
settings  on  each  side  of  PSC  and  the  PSC  profile  data  are  linearly  interpolated. 
For  many  PSC'S,  the  reference  profile  data  for  at  least  one  dataset  will  have 
been  computed  for  the  previous  PSC  and  thus  stored  in  the  SENX  array. 

(8)  If  lEDIT  equal  zero,  the  PSC  profile  data  are  edited  to  select  the  10  angles  which 
best  represent  the  profile  dataset  at  the  reference  distance. 

(9)  For  BEDIT  greater  than  or  equal  to  zero,  all  PNLTX,  ALX  and/or  ALTX  profile  data 
requested  by  the  MEAS  program  flag  are  written  to  the  file  opened  on  Unit  2. 

(10)  For  IPR  greater  than  zero,  the  PNLX,  PNLTX,  ALX  and/or  ALTX  profile  data  are 
printed  as  requested  by  the  MEAS  flag.  A  plot  is  also  printed  for  PNLTX,  ALX  and 
ALTX  at  the  reference  distance. 

(1 1)  After  ail  data  are  analyzed  for  all  aircraft  or  when  an  error  occurs,  the  program 
terminates  with  a  call  to  EXIT  where  the  argument  is  one  when  an  error  has  occurred 
or  zero  when  there  is  no  error. 

The  content  and  format  of  the  above  mentioned  output  can  best  be  determined  by 
consulting  the  sample  problem  in  Appendix  D. 

The  CDC  FORTRAN  extended  (FORTRAN  IV)  computer  language  was  used  for  the 
original  version  of  this  program.  The  common  and  subroutine  features  of  the  language  were  used 
extensively  throughout  the  program  to  save  computer  time  and  memory. 

Over  the  past  ten  years,  this  program  was  converted  from  CDC  FORTRAN  to 
FORTRAN  77  for  the  Perkin-Elmer  and  later  to  Microsoft  and  Lahey  FORTRAN  for  the  PC. 
Version  1 1.3  is  currently  being  compiled  with  the  Lahey  compiler. 

The  following  sections  describe  the  detailed  tasks  accomplished  by  the  program.  It  is 
intended  to  document  the  procedures  within  each  subroutine  at  a  level  useful  to  either  a 
programmer  reading  this  while  working  with  the  code  or  a  reader  simply  interested  in  what 
happens  with  a  specific  subroutine.  The  algorithms  used  to  compute  the  noise  measures  are 
described  in  detail  in  the  individual  subroutines. 
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GENERAL  ORGANIZATION  OF  THE  OMEGA  II  PROGRAM 


The  general  organization  of  the  entire  program  is  shown  in  Figure  4.  The  diagram 
indicate  access  to  the  various  routines  rather  than  program  flow;  for  example,  MAIN  calls 
subroutine  CDIST  which  in  turn  calls  subroutine  CPNL,  CPTC  and  CAL  and  subroutine  CPNL 
calls  function  FNOY.  In  general,  program  flow  is  from  left  to  right  in  this  diagram. 

Using  Figure  4  as  a  guide,  this  section  summarizes  in  very  general  terms  the  functions 
performed  by  the  entire  program.  This  is  meant  to  serve  as  an  introduction  for  the  reader  to  the 
functions  of  the  individual  subroutines. 

The  control  routine,  MAIN,  reads  the  job  control  card  and  initializes  several  program  and 
test  variables.  The  required  input  and  output  files  are  opened  in  routines  CLFILE  and  RDFILE. 
Subroutine  TESTN  is  then  called  to  read  the  code  sheet  test  parameters,  initialize  numerous  test 
variables,  and  call  subroutine  ALPH  to  compute  the  atmospheric  absorption  data  for  non¬ 
standard  temperature  and  relative  humidity. 

Next  the  MAIN  routine  initializes  additional  test  parameters  and  calls  subroutine  RSPLN 
to  read  all  reference  datasets  for  aircraft  ACC  from  the  file  opened  on  Unit  7.  The  power  settings 
from  these  reference  datasets  as  well  as  the  output  power  setting  data  are  ranked  by  subroutine 
RANK  which  also  determines  the  reference  datasets  required  to  compute  the  profile  data  for  each 
output  power  setting. 

The  cover  (IPR>0)  and  summary  pages  are  printed  by  subroutines  TITPG  and  SUMRY, 
respectively.  Delta  N  (DELN)  is  added  to  all  reference  spectra  (for  DELNtK))  which  are  then 
printed  (IPR>0)  by  the  MAIN  routine. 

Subroutine  CDIST  is  called  to  compute  the  PNLX,  PNLTX,  ALX  and/or  ALTX  profile 
data  (as  requested  by  the  MEAS  flag)  for  the  reference  dataset  power  settings  required  to 
interpolate  the  output  power  settings  (PSC).  This  interpolation  between  power  settings  is 
performed  by  the  MAIN  routine.  Subroutine  CDIST  calls  subroutine  CPNL,  CPTC  and  CAL  to 
compute  the  PNLX,  tone  correction,  and  ALX  data.  These  interpolated  profile  data  are  edited 
(IEDIT=  0)  by  subroutine  EDIT  to  select  the  10  angles  which  best  describe  the  profile  data  at  the 
reference  distance.  EDIT  calls  subroutine  ERR  which  uses  linear  interpolation  between  the 
selected  angles  to  determine  how  well  they  represent  the  profile  function.  The  FINTP  linear 
interpolation  subroutine  is  used  by  subroutine  ERR.  EDIT  also  calls  subroutine  ITER  which 
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Figure  4.  General  Organization  of  the  OMEGA  1 1  Program. 


attempts  to  reduce  the  interpolation  errors  by  selecting  the  angles  with  the  largest  errors;  ITER 
also  calls  ERR. 

Subroutine  PPFDAT  writes  the  profile  datasets  to  the  file  opened  on  Unit  2  (for 
IEDIT>0)  and  writes  (IPR>0)  the  list  output  for  all  computed  measures  to  Unit  6..  Subroutine 
PLT  prints  (IPR>0)  a  plot  of  the  PNLTX,  ALX  and  ALTX  noise  level  versus  angle  data 
computed  for  the  reference  distance.  The  angles  for  which  profile  data  are  written  on  Unit  2  are 
listed  below  the  plot. 

After  computing  and  printing  the  profile  data  for  each  of  the  NP  power  settings,  control  is 
returned  to  label  10  in  the  MAIN  routine  and  the  program  repeats  the  above  for  each  aircraft  until 
an  end  of  file  is  read  from  the  INPUT  file. 
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DEFINITIONS  OF  SYMBOLS  AND  TERMINOLOGY 
USED  IN  THE  OMEGA  II  PROGRAM 


The  symbols  defined  here  are  used  in  this  report  and/or  in  the  OMEGA  1 1  program 
source  listing.  They  are  a  subset  of  the  complete  symbol  versus  reference  list.  Many  of  the 
symbols  used  in  this  program  are  really  dummy  variables  used  in  only  one  or  two  routines  and 
redefined  in  each  routine;  most  of  these  symbols  are  not  included  in  this  list  of  symbol 
definitions.  Symbols  which  are  arrays  will  be  listed  with  their  array  dimensions.  Variables  I,  J 
and  K  are  usually  but  not  always  used  as  array  subscripts  as  follows: 

(1)  The  subscript  "I"  is  a  running  index  associated  with  any  one  spectrum  (angle).  It  is 
also  frequently  used  as  a  "dummy"  index  to  initialize  variables. 

(2)  The  subscript  "J"  is  a  mnning  index  associated  with  any  one  band  in  the  set  of  octave 
or  1/3  octave  frequency  bands.  It  is  also  an  index  associated  with  a  specific  profile 
measure. 

(3)  The  subscript  "K"  is  a  running  index  associated  with  any  one  profile  distance. 
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ACC 

AL'rx:(  19,22,2) 

ALX(19,22,2) 

ARCS 

ASK 

ATN(31) 

ATN8(31) 

ATNC(31) 

AW(31) 

BLK 

CBLK 

C0MD(6) 

CRI 

CXD(19) 

DATE 

DATN(6) 

DELN 

DIST 

EA(22,24) 

ER(19,3) 


Three  character  aircraft  code  read  from  the  code  sheet  and  printed  on  all  output 
pages  and  in  all  output  datasets. 

Profile  tone-corrected  A-weighted  overall  sound  level  in  dBA  for  each  angle  and 
profile  distance  for  two  power  settings. 

Profile  A-weighted  overall  sound  level  in  dBA  for  each  angle  and  profile  distance 
for  two  power  settings. 

Character  variable  used  to  retrieve  the  command  line  argument  list. 

One  byte  character  variable  which  contains  an  asterisk.  ASK  is  used  to  check  for 
the  COMDECK  name  record.  Also  used  as  end  of  file  flag. 

Data  statement  array  containing  the  atmospheric  absorption  coefficients  in  dB  per 
1000  for  standard  day  temperature  and  relative  humidity. 

Atmospheric  absorption  coefficients  in  dB  per  1000  feet  for  OMEGA  1 1  profile 
output  temperature  and  relative  humidity. 

Atmospheric  absorption  coefficients  in  dB  per  1000  feet  for  OMEGA  1 1  reference 
input  temperature  and  relative  humidity  (usually  standard  day  conditions). 

A-weighting  coefficients  in  dBA. 

A  data  statement  variable  containing  a  blank  Hollerith  character  used  in  printing 
variable  format  data. 

Ten  byte  character  variable  which  is  filled  with  ten  spaces  or  blank  characters.  It 
is  used  throughout  the  program  to  blank  out  character  data. 

Character  variable  used  to  store  the  last  four  characters  of  each  reference  file 
COMDECK  name. 

Comdeck  revision  identifier  (see  code  sheet). 

Computed  tone  correction  for  each  angle  for  the  reference  distance  (dB). 

Date  of  computer  run  (see  code  sheet). 

Date  of  the  OMEGA  8  computer  run  which  generated  the  reference  dataset. 

Noise  adjustment  factor  added  to  the  reference  spectrum  (dB). 

Distance  in  feet  to  which  the  reference  file  data  in  NOISEFILE  are  normalized 
(presently  250  feet). 

Excess  atmospheric  attenuation  in  dB  for  bands  17  to  40  and  distances  200  to 
25000  feet  (subroutine  CDIST). 

Profile  dataset  angle  selection  error  data  for  each  angle  and  measure  (dB). 
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SYMBOL 

ERMAX 

FIMPR8 

FJ 

FL(24,5) 

FM(24,4) 

FMT(22) 

FMXER 

FREQ(31) 

FREQ3(24) 

FSPL(  19,3 1,6) 

IBNH 

EBNL 

IC 

ICC 

ID 

EDIT 

EC(6) 

ECC 

in(6) 

DFH 


DESCRIPTION 

Maximum  acceptable  angle  selection  error  in  dB;  no  attempt  is  made  to  improve 
the  angle  selection  for  errors  less  than  ERMAX. 

Characteristic  impedance  ratio  using  reference  and  profile  output  temperature  and 
barometric  pressure. 

Constant  used  in  the  perceived  noise  level  computations;  FJ=0.15  for  1/3  octave 
band  data. 

Data  statement  array  used  in  noy  computations  (functions  FNOY)  containing  the 
band  sound  pressure  levels  in  dB  given  in  Table  4. 

Data  statement  array  used  in  noy  computations  (function  FNOY)  containing  the 
reciprocals  of  the  slopes  given  in  Table  4. 

Variable  format  array  used  to  print  profile  data  in  subroutine  PPFDAT. 

Maximum  angle  selection  error  permitted  without  an  error  message  being  printed 
(see  code  sheet). 

Data  statement  array  containing  the  frequency  values  in  Hz  in  character  format  for 
printing. 

Geometric  mean  and  lower  limiting  frequencies  required  to  compute  atmospheric 
absorption  coefficients  for  1/3  octave  band  data. 

Normalized  reference  SPL  in  dB  for  each  angle,  frequency  band,  and  operation 
power  code. 

Largest  band  number  index  (31  corresponds  to  band  40). 

Initial  band  number  index  (1  corresponds  to  band  10). 

Index  which  is  usually  associated  with  output  power  setting  (PSC)  data. 

Integer  variable  containing  the  power  setting  index  for  the  specific  power  setting 
being  analyzed. 

Integer  variable  used  to  store  the  integer  value  of  the  standard  reference  distance. 

Program  flag  which  controls  the  qu2intity  of  profile  data  written  to  the  file  opened 
on  Unit  2  (see  code  sheet). 

Program  flag  used  to  flag  special  case  output  power  setting  data  (see  code  sheet). 
Index  used  to  count  the  number  of  EC's  >  0. 

Program  flag  used  to  flag  special  case  (Afterburner,  etc.)  reference  file  power 
setting  data. 

Index  used  to  count  the  number  of  El's  >  0. 
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SYMBOL 

IH 

IH8 

IHH 

n 

IL 

IPAGE 

IPR 

IPRCK(6) 

IR(19) 

IRD 

IREQ(2,6) 

IT 

ITS 

IVER 

L 

M 

MEAS(3) 

MM 

N 

NC 

NFRUNl 

NFRUN2 


DESCRIPTION 

Index  of  the  largest  frequency  band  defined  for  this  test  (IH=3 1  for  band  40). 
Profile  output  relative  humidity  (code  sheet  input). 

Reference  relative  humidity  (70%). 

Index  which  is  frequently  associated  with  the  computation  of  the  two  sets  of 
profile  data  used  to  interpolate  profile  data  at  the  requested  power  settings. 

Index  of  the  lowest  frequency  band  defined  for  this  test  (IL=1  for  band  10). 

Page  number  counter. 

Program  print  control  flag  (see  code  sheet). 

Integer  variable  used  to  flag  the  reference  data  required  to  compute  the  requested 
output. 

Array  used  to  print  integer  values  of  the  reference  SPL  data. 

Profile  distance  index  which  corresponds  to  the  reference  minimum  slant  range 
(presently  IRD=2). 

Index  of  the  one  or  two  reference  datasets  required  to  compute  the  profile  data  for 
each  requested  power  setting. 

Reference  temperature  in  degrees  F  (always  59°F). 

Profile  output  temperature  in  degrees  F  (see  code  sheet). 

Program  version  code. 

Index  frequently  associated  with  data  for  a  specific  reference  file  power  setting. 

Program  flag  indicating  whether  data  are  for  octave  (M=l)  or  1/3  octave  (M=3) 
band  center  frequency  (always  3  in  this  program). 

Program  job  control  variable  which  flags  the  profile  measures  to  be  computed  (see 
code  sheet). 

Increment  of  the  frequency  index  (one  for  1/3  octave  band  data). 

Number  of  reference  dataset  operation  power  codes  read  from  the  reference  file 
(Unit  7)  for  aircraft  ACC. 

Number  of  angles  for  which  SPL  spectra  are  defined  in  the  reference  dataset 
(always  19  for  the  present  reference  datasets). 

Twenty-four  byte  character  variable  used  to  store  the  first  reference  filename. 
Twenty  four  byte  character  variable  used  to  store  the  second  reference  filename. 
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SYMBOL 


DESCRIPTION 


NN 

NP 

NPM 

NR(17,3) 

NRC(6) 

NRI(6) 

0PC(6) 

0PCC(6) 

OPCDM 

OPCSP(14) 

OPD(2,6) 

OPDl,  OPD2 

PI 

P8 

PNLTX(19,22,2) 

PNLX(19,22,2) 

PS(6,6) 

PSC(6) 

PSCF(6) 


Maximum  number  of  reference  file  operation  power  codes  permitted  for  each 
aircraft  (NN=6). 

Number  of  output  operation  power  codes  (power  settings)  to  be  computed  for  this 
aircraft. 

maximum  number  of  profile  operation  power  codes  permitted  for  each  aircraft 
(NPM=6). 

Change  in  slope  rank  data  computed  in  the  angle  selection  routines  for  angles  10  to 
170  degrees  for  each  measure.  Angles  with  rank  greater  than  nine  are  included  in 
the  final  profile  dataset. 

NRC(K)  contains  the  index  of  the  output  power  setting  data  (in  arrays  PSC  and 
PSCF)  whose  power  setting  rank  is  K. 

NRI(K)  contains  the  index  of  the  reference  power  setting  data  (in  arrays  PS  and 
PSIF)  whose  power  setting  rank  is  K. 

Reference  operation  power  codes. 

Output  operation  power  codes. 

Two  byte  character  variable  containing  the  output  operation  power  code. 

Special  case  operation  power  codes  for  which  no  inteipolation  is  permitted.  This 
array  must  be  updated  when  new  codes  are  added. 

Power  description  data  for  each  reference  operation  power  code  (20  characters). 

Power  description  data  defined  for  output  power  settings  when  the  reference  and 
output  power  settings  are  the  same. 

Reference  barometric  pressure  in  inches  Hg  (always  29.92  inches  Hg). 

Profile  output  barometric  pressure  in  inches  Hg. 

Profile  tone-corrected  perceived  noise  level  in  PNdB  for  each  angle  and  distance 
for  two  power  settings. 

Profile  perceived  noise  level  in  PNdB  for  each  angle  and  distance  for  two  power 
settings. 

Power  setting  data  for  each  reference  file  operation  power  code  (power  setting 
value  and  units  in  character  format). 

Power  setting  data  for  each  profile  output  operation  power  code  in  character 
format  (see  code  sheet). 

Numeric  form  of  PSC(6)  data  defined  above. 
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SYMBOL  DESCRIPTION 

PSIF(6)  Numeric  form  of  the  first  power  setting  for  each  reference  operation  power  code  in 

PS(6,6)  above;  e.g.,  PSEFQ  =PS(1, 1). 

PSU  Power  setting  units  for  the  profile  output  data. 

PTC  Tone  correction  in  dB. 


PV 

RMS(3) 

RUN(6) 

RUNC(6) 

SENX(19,22,12) 

SPLX(  19,31) 
SX(22) 

T12(2) 

TEST(6) 

TT(6.6) 

TYPE(3) 


Profile  version  code  (see  code  sheet). 

Root  mean  square  of  the  angle  selection  error  for  each  measure. 

Two  character  run  number  from  each  reference  dataset. 

Data  statement  array  containing  run  numbers  01  to  06  in  character  format.  These 
run  numbers  are  assigned  to  the  output  data  in  power  setting  sequence. 

Array  equivalent  to  arrays  PNLX  through  SENXD  in  blank  common  where 
SENXD(  19,22,4)  contains  the  interpolated  profile  data  for  the  measures. 

Profile  sound  pressure,  level  in  dB  for  each  angle  (spectrum). 

Distance  data  in  feet  for  the  22  profile  distances. 

Data  statement  array  used  to  set  up  the  variable  formats  for  the  profile  data  tab 
output. 

Ten  character  test  number  from  each  reference  dataset. 

Two  lines  of  25  characters  from  the  reference  dataset  describing  the  noise  source. 
The  two  lines  from  the  first  reference  dataset  are  printed  in  the  Noise 
Source/Subject  block  on  all  output  pages. 

Five  byte  character  variable  containing  the  PNLT,  AL  and  ALT  noise  measure 
type  names  for  printing. 
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DETAILED  DESCRIPTION  OF  THE  OMEGA  11  PROGRAM 


This  section  discusses  the  MAIN  routine  and  each  subroutine  in  the  OMEGA  1 1  program. 
Procedures  within  most  routines  are  documented  at  a  level  useful  to  a  programmer  reading  this  while 
working  with  the  code  or  a  reader  interested  in  what  happens  within  a  specific  subroutine.  Most  routines 
contains  numerous  comments  which  should  be  very  helpful  in  following  the  code. 

The  program  algorithms  and  I/O  are  discussed  in  the  routines  in  which  they  are  coded.  The 
program  code  sheet  in  Appendix  B  and  the  sample  problem  in  Appendix  D  are  referenced  to  simplify  the 
description  of  the  input  and  output.  The  more  complex  subroutines  are  supplemented  by  flowcharts 
drawn  from  the  point  of  view  of  function  performed  rather  than  block  instructions. 

COMMON  VARIABLES 

Extensive  use  is  made  of  common  in  the  program  for  communications  between  the  various 
routines.  Many  of  the  storage  locations  in  blank  common  are  used  in  different  ways  or  with  different 
variable  names  throughout  the  program.  This  was  done  to  reduce  the  memory  required  to  run  the 
program  on  the  old  CDC  computers.  Several  of  the  large  arrays  are  included  in  blank  common  rather 
than  labeled  common  because  on  the  CDC  6600  computer  a  large  blank  common  reduces  the  total 
memory  required  to  load  and  execute  the  program.  The  variables  used  in  labeled  coirunon  are  usually  of 
a  similar  type  and/or  used  in  many  of  the  same  subroutines. 

The  variables  assigned  to  blank  and  labeled  common  in  the  MAIN  routine  and  the  total 
common  length  are  listed  in  Table  6.  The  subroutines  in  which  the  labeled  common  are  used  are 
listed  in  Table  7.  A1  blank  and  labeled  common  are  included  in  the  MAIN  routine.  All  common 
variables  are  defined  in  the  complete  list  of  symbols  at  the  beginning  of  the  report.  The  blank 
and  labeled  common  are  described  in  the  following  paragraphs. 

CHAR  Common 

All  the  variables  in  this  labeled  common  are  character  variables.  Most  of  the  variables 
were  formerly  part  of  HEADC  labeled  common  before  all  character  variables  were  put  in  this 
CHAR  labeled  conunon  when  the  program  was  converted  to  run  on  a  PC.  Most  of  these 
variables  are  used  in  the  list  and  COMDECK  output  to  identify  the  data.  There  are  456  bytes  in 
this  labeled  common. 
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TABLE  6 

MAIN  ROUTINE  VARIABLES  IN  BLANK  AND  LABELED  COMMON 


Length  in 
Bytes 


LABELED  COMMON 


Blank  Common 

ATTC 

HE  ADC 

CHAR 

ERRORS 

M 

ATNC(31) 

DELN 

ACC 

ERRFLG 

MM 

ATNR(31) 

FIMPR8 

ASK 

IL 

SX(22) 

IC 

BLK 

m 

ICC 

CBLK 

NC 

IFC(6) 

COMD(6) 

L 

IFCC 

CRI 

N 

in(6) 

DATE 

ID 

IFE 

DATN(6) 

Disr 

niH 

OPC(6) 

ms 

OPCC(6) 

FSPI.(19.31-6) 

IPAGB 

OPCDM 

Sr'LX(  19,22) 

IVT^R 

OPD(2,6) 

PNLX(  19,22,2'! 

rr 

OPDl 

PNLTX(  19,22,2) 

IT8 

OPD2 

ALX(  19,22,2) 

NP 

PSU 

ALTX(  19,22,2) 

PI 

PS  (6,6) 

SENXD(  19.22,4) 

P8 

PSC(6) 

rR(19) 

psm(6) 

PV 

IDPRCK(6) 

PSCr(6) 

RUN(6) 

DM\(202) 

NRC(6) 

RUNC(6) 

TEST(6) 

TT(6,6) 

36,S28 

336 

ISO 

456 

4 
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TABLE  7 

SUBROUTINES  CONTAINING  THE  LABELED  COMMON 


ATTC 

CHAR 

CDIST 

BLOCK  DATA 

TESTN 

HEADS 

PLT 

PPFDAT 

RANK 

RSPLN 

SUMRY 

TESTN 

TITPG 

HEADC 

ERRORS 

BLOCK  DATA 

CLFILE 

HEADS 

RANK 

PPFDAT 

RDFILE 

RANK 

RSPLN 

RSPLN 

TESTN 

SUMRY 

TESTN 

TITPG 
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ERRORS  Common 

This  common  contains  the  ERRFLG  which  is  used  to  set  the  exit  code.  ERRFLG  is  set  to 
TRUE  when  there  is  an  error. 

Blank  Common 

The  36,828  bytes  used  by  blank  common  in  the  MAIN  routine  are  the  maximum  required 
in  any  routine  in  the  program.  Almost  all  routines  use  some  blank  common  but  only  a  few 
require  the  36,828  bytes.  The  variable  names  assigned  to  blank  common  vary  throughout  the 
program.  Most  variables  are  defined  to  communicate  with  several  subroutines  and  then 
redefined  for  the  next  series  of  routines.  In  several  cases  variables  are  changed  or  equivalenced 
to  simplify  the  computations  and/or  output  routines. 

ATTNC  Common 

These  variables  are  used  to  extrapolate  data  from  a  given  set  of  distance  and  weather 
conditions  to  a  new  set  of  conditions.  For  example,  from  reference  distance  and  standard  day 
conditions  to  the  profile  distances  at  the  profile  temperature  and  humidity.  These  336  bytes  are 
defined  the  same  throughout  the  program. 

HE  A  DC  Common 

These  numeric  variables  are  used  to  store  the  weather  and  power  setting  data  required  to 
perform  the  interpolation  at  the  output  conditions.  Some  of  these  data  are  printed  in  the  page 
header  blocks  and  the  output  COMDECK  comment  cards.  The  180  bytes  in  this  labeled 
common  are  defined  the  same  throughout  the  program. 

MAIN  PROGRAM 

MARI  is  the  executive  routine  for  the  entire  OMEGA  1 1  program.  Its  principal  funcUon 
is  to  call  the  subroutines  required  to  read  the  input  data,  perform  the  data  analysis,  and  write  the 
list  and  data  file  outputs.  However,  the  program  also  reads  the  program  control  parameters, 
initializes  numerous  test  variables,  interpolates  the  profile  noise  levels  as  a  function  of  power 
setting,  and  prints  part  of  the  output.  The  MAIN  routine  is  discussed  in  the  following  paragraphs 
using  the  program  listing  in  Appendix  J,  the  flowchart  in  Figure  5,  and  the  sample  problem  in 
Appendix  D  as  guides. 

Data  Statement  Arrays 

The  FREQ  data  statement  array  contains  frequency  data  (Hz)  in  character  format  which  is 
printed  in  the  frequency  versus  SPL  listings  of  the  reference  datasets. 
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Figure  5.  Flowchart  for  OMEGA  11  MAIN  Routine. 
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Initialize  Test  Variable 


The  first  segment  of  this  routine  performs  the  following  variable  initialization  steps; 

(1)  The  Lahey  FORTRAN  routine  GETCL  is  called  to  get  the  command  line 
argument  which  is  stored  in  character  variable  ARCS.  If  ARCS  is  NULL  tlien 
subroutine  RDFILE  is  called  to  open  and/or  check  the  I/O  files.  If  ARCS  is  not 
NULL  then  subroutine  CLFILE  is  called  to  get  the  I/O  file  names  from  ARCS 
and  then  open  and/or  check  the  files. 

(2)  The  SX  standard  profile  distance  array  is  evaluated: 

SX,  =  antilog 

where  I  is  the  distance  index  running  from  1  to  22. 

(3)  The  code  sheet  job  control  record  is  read  once  per  execution  of  the  program. 

(4)  Subroutine  TESTN  is  called  to  read  all  test  (aircraft)  code  sheet  parameters  and 
initialize  numerous  other  test  variables  based  on  these  code  sheet  input  data  (e.g., 
atmospheric  absorption  data). 

Read  Reference  Data 

The  ptogram  calls  subroutine  RSPLN  to  read  the  reference  data  for  aircraft  ACC  from  the 
file  opened  on  Unit  7.  This  Unit  7  file  may  contain  data  for  numerous  aircraft,  however,  all  data 
for  aircraft  ACC  must  be  back  to  back.  All  reference  data  (for  ACC)  for  normal  power  settings 
are  stored  in  the  program;  the  special  case  (afterburner,  etc.)  power  setting  data  are  stored  only 
when  at  least  one  special  case  output  power  setting  is  requested.  Note  that  array  OPCSP  in 
subroutine  RSPL  must  be  updated  when  new  special  case  OPC's  are  added  to  NOISEFELE. 

Setup  Interpolation  Index  Array 

Subroutine  RANK  is  called  to  rank  from  low  to  high  the  input  and  output  power  setting 
data.  These  ranked  power  setting  data  are  used  to  determine  the  reference  data  required  to 
interpolate  the  profile  data  for  each  output  power  setting.  The  indices  of  these  reference  data  are 
stored  in  array  EREQ. 

Subroutine  TITPG  prints  the  cover  page  when  flag  DPR  is  greater  than  zero.  Subroutine 
SUMRY  uses  the  abo^'e  IREQ  array  data  to  print  the  aircraft  summary  page  which  is  the  only  list 
output  when  IPR  is  less  than  one. 
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Adjust  and  Print  Reference  Data 

The  following  constant  is  computed  as  required  by  subroutine  CDIST  later  in  the 
program: 


D1  =  10  log  (FIMPR8)  +  20  log  (DIST) 


where 


FIMPR8  =  the  characteristic  impedance  ratio  computed  in  subroutine  TESTN. 

DIST  =  the  reference  distance  in  feet. 

The  statement  label  55  loop  checks  array  IREQ  to  determine  which  reference  data  are 
actually  required  for  this  aircraft  analysis.  If  Delta  N  is  not  equal  to  zero,  it  is  added  to  all 
spectra  in  these  required  datasets  which  are  then  printed  in  the  label  150  loop  when  flag  IPR  is 
greater  than  zero.  Subroutine  HEADS  prints  each  reference  page  header  block.  The  overall  SPL 
is  also  computed  and  printed  below  the  last  SPL  band.  The  overall  SPL  is  defined  as  follows: 

Overall  SPL  =  10  lo 

where 

FSPLj  j  L  “  reference  band  SPL  in  dB  at  distance  DIST  and  for  angle  index  I,  band 
index  J,  and  power  setting  index  L. 

Single  Event  Profile  Data 

The  statement  label  500  loop  computes  the  single  event  profile  noise  data  for  each  of  the  NP 
power  settings  (PSC).  Each  of  the  NP  power  settings  will  require  profile  data  from  one  reference  dataset 
when  the  PSC  is  the  same  as  a  reference  power  setting  or  from  two  reference  datasets  when  the  PSC 
must  be  interpolated.  The  dataset  indices  required  for  each  PSC  are  stored  in  array  IREQ.  Subroutine 
CDIST  computes  the  profile  data  requested  by  the  MEAS  flag  (PNLTX,  ALX  and/or  ALTX)  for  each 
reference  dataset. 

Note  that  profile  data  are  computed  and  stored  for  no  more  than  two  reference  datasets  at  any  one 
time;  however  array  IREQ  is  set  up  and  checked  to  avoid  duplicate  computations.  This  limit  of  two  was 
set  to  keep  the  program  memory  requirements  within  reasonable  limits  for  the  old  CDC  computers. 

When  no  power  setting  interpolation  is  required,  the  final  profile  data  are  stored  in  array 
SENX  (by  subroutine  CDIST)  as  indicated  in  the  program  comment  statements  (see  label  370  and 
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380).  When  the  PSC  are  not  the  same  as  a  reference  file  power  setting,  the  profile  data  are 
interpolated  for  each  angle  (I)  and  distance  (K): 


SENXj  f^  j 


SENX/  —  SENXj  jj 
PSIF,,-PSIF,, 


[PSCF,cc  -  PSIF,,  ]  +  SENX,_,,jj 


PNdB  or  dBA 


where 

SENXi,K,J 

SENXi,K,J1 


SENXi,k,JJ 


the  final  profile  data  in  PNdB  or  dBA  for  the  angle,  distance  and 
measure  index  J  (J  ranges  from  9  to  12). 

the  profile  data  for  the  Jl^^  measure  index  from  the  second  (L2) 
reference  dataset  for  I  and  K. 

the  profile  data  for  the  measure  index  from  the  first  (LI)  reference 
dataset  for  I  and  K. 


PSIFl2  “  power  setting  from  the  second  (L2)  reference  dataset. 

PSEFli  =  power  setting  from  the  first  (LI)  reference  dataset. 

PSCFj(2(2  =  the  output  profile  power  setting  for  which  the  profile  data  are 
interpolated. 


Exiit  and  Write  Profile  Data 

When  flag  EEDIT  is  zero,  subroutine  EDfT  is  called  to  select  the  ten  angles  which  best 
describe  each  profile  measure  at  the  reference  distance.  The  ten  selected  angles  are  angles  0  and 
180  degrees  plus  the  eight  angles  from  10  to  170  degrees  for  which  the  value  of  array  NRj  j  is 
greater  than  nine,  where  I  is  the  angle  index  (1=1  to  17)  and  J  is  the  measure  index  (J=l  to  3  for 
PNLTX,  ALX  and  ALTX). 


Subroutine  PPFDAT  is  called  to  write  the  final  profile  data  for  the  requested  measures 
(MEAS>0)  to  the  file  opened  on  Unit  2  and  to  print  the  profile  data  list  output.  The  quantity  of 
data  written  to  Unit  2  depends  on  the  value  of  the  EEDIT  flag: 

(1)  For  IEDrr<0,  all  Unit  2  output  are  omitted. 

(2)  For  IEDIT=0,  write  the  10  angles  selected  by  subroutine  EDIT. 

(3)  For  EEDn'=l,  write  all  19  angles  of  profile  data. 

A  description  of  the  content  and  format  of  these  profile  datasets  is  given  in  Appendix  H. 
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The  profile  data  listings  are  printed  when  program  flag  BPR  is  greater  than  zero.  They 
contain  data  for  all  19  angles  and  22  distances;  one  page  is  printed  for  each  measure  (listing 
includes  the  PNLX  data  when  PNLTX  is  computed). 

Subroutine  PLT  prints  (IPR>0)  a  plot  of  PNLTX,  ALX  and  ALTX  noise  level  versus 
angle  for  the  reference  distance.  This  plot  is  very  helpful  when  selecting  the  ten  angles 
(manually)  or  when  checking  the  angles  selected  by  the  edit  routine.  The  angles  written  to  Unit  2 
for  each  measure  are  printed  below  this  plot. 

At  this  point,  control  is  returned  to  label  10  to  begin  the  next  aircraft  analysis  or  to 
terminate  the  job.  An  OMEGA  1 1.3  job  will  terminate  with  a  call  to  EXTr(l)  if  an  error  has 
occurred  (ERRFLG  =  TRUE)  or  a  call  to  EXIT(O)  when  the  job  terminates  without  an  error 
(ERRFLG  =  FALSE).. 

SUBROUTINE  TESTN(NPM) 

This  subroutine  is  called  from  the  MAIN  routine  to  input  the  aircraft  code  sheet 
parameters.  The  subroutine  also  initializes  numerous  test  variables  and  sets  code  sheet  default 
parameters.  There  should  be  sufficient  comment  statements  in  the  listing  to  document  most  of 
the  coding. 

Method 

The  first  segment  of  the  subroutine  (to  label  10)  reads  the  code  sheet  test  parameters  and 
sets  the  program  default  values.  The  computer  job  is  terminated  when  an  end  of  file  is  read  from 
Units.  The  code  sheet  parameters  and  default  conditions  are  described  in  Appendix  B.  The 
reference  temperature  (ITS)  is  converted  from  °F  to  °C  as  required  to  compute  the  impedance 
data. 


The  impedance  ratio  for  reference  and  profile  output  conditions  is; 

FIMPR2,  =  I 


■  273  + TM  ■ 

1/2 

'PS' 

.273  +  TMS. 

-PI. 

where 


TM  =  reference  temperature  (15°C) 

TM8  =  profile  output  temperature  in  °C 

PI  =  reference  barometric  pressure  (29.92  in  Hg). 

PS  =  profile  output  barometric  pressure  in  inches  Hg. 
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Subroutine  ALPH  is  called  to  compute  the  atmospheric  absorption  data  in  dB  per  1000 
feet  for  non-standard  day  profile  conditions  (ATN8).  For  standard  day  conditions  (59°F  and 
70%),  these  absorption  data  are  stored  in  the  ATN  data  statement  array. 

SUBROUTINE  ALPH(REL,  TEMP^BC) 

This  subroutine  is  called  from  subroutine  TESTN  to  compute  the  octave  or  1/3  octave 
band  atmospheric  absorption  data  in  dB  per  1000  feet  for  non-standard  day  temperature  and 
relative  humidity.  This  subroutine  is  the  same  as  subroutine  ALPH  in  the  OMEGA  10  program. 

FUNCTION  ATKN(X,  Y,N,K,XI) 

This  function  is  a  general  ATTKEN  interpolation  function,  used  by  subroutine  ALPH  to 
compute  the  normalized  molecular  absorption  coefficient.  ATKN  was  obtained  from  the  ASD 
computer  center  library  (old  IBM  7094  library).  Since  this  is  a  common  interpolation  function 
defined  in  most  numerical  methods  texts,  no  additional  description  will  be  given. 

SUBROUTINE  HEADS(IPH) 

This  subroutine  is  called  from  the  MAIN  routine  and  from  subroutines  CDIST  and  PLT 
to  print  page  header  blocks  at  the  top  of  all  output  pages.  The  header  blocks  are  1 12  characters 
wide  for  the  reference  sound  pressure  level  output  (label  2000  formats)  and  126  characters  wide 
for  all  profile  output  pages  (label  3000  formats).  The  BPH  subroutine  argument  determines  the 
specific  page  header  block  printed  for  each  call. 

The  content  and  format  of  the  header  blocks  can  best  be  observed  by  consulting  the 
sample  problem  in  Appendix  D.  This  subroutine  should  contain  sufficient  comment  statements  to 
identify  the  data  being  printed.  It  may  also  be  helpful  to  consult  the  subroutine  HEAD 
documentation  in  the  OMEGA  10  program;  the  five  categories  of  data  described  there  also  apply 
here. 

SUBROUTINE  RSPLN(NNJERR) 

This  subroutine  is  called  from  the  MAIN  routine  to  input  one  reference  dataset  from  the 
file  opened  on  Unit  7  for  each  normal  operation  power  code  (OPC)  in  the  reference  file 
(NOISEFILE )  for  aircraft  ACC.  When  the  Unit  7  file  contains  more  than  one  reference  dataset 
for  one  or  more  OPCS,  only  the  last  dataset  is  stored  by  this  routine.  The  special  case  operation 
power  code  data  (Afterburner,  Wet  or  With  Jets)  are  read  only  when  at  least  one  special  case  was 
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requested  on  the  code  sheet  (IFCC>0).  The  special  case  operation  power  codes  are  stored  in 
array  OPCSP  which  must  be  updated  when  new  codes  are  added  to  the  NOISEFTLE  database. 
NN  is  the  maximum  number  of  datasets  for  which  storage  has  been  allocated  in  the  program 
(NN=6).  The  lERR  argument  is  a  program  error  flag  which  returns  codes  of  1  to  4  for  different 
type  input  errors. 

Band  10  to  16  were  added  to  all  reference  data  in  1992.  The  SPL  levels  are  zero  for  these 
seven  bands  for  all  current  runup  data.  The  label  150  loop  checks  these  seven  bands  for  all 
spectra  for  the  current  aircraft.  Variable  IL  is  set  equal  to  the  band  index  of  the  first  non  zero 
band. 


The  flowchart  in  Figure  6  and  the  comment  statements  in  the  listing  should  make  it  easy 
to  follow  the  data  being  read;  the  format  and  content  of  these  ground  run-up  reference  datasets 
are  also  given  in  Appendix  G. 

FUNCTION  ICV(R) 

This  function  is  used  in  the  MAIN  and  PLT  routines  to  convert  variable  R  to  an  integer. 

R  is  rounded  up  when  the  fractional  part  is  greater  than  or  equal  to  0.5.  Function  ICV  could  be 
replaced  by  function  NIOT  which  is  a  standard  FORTRAN  77  function. 

SUBROUTINE  CDIST(IRD,DI,II) 

This  subroutine  is  called  from  the  MAIN  routine  to  compute  the  single  event  noise  profile 
data  (e.g.,  PNLTX  versus  distance  for  22  distances  from  200  to  25,000  feet)  from  the  reference 
SPL  spectra  stored  in  array  FSPL. 

The  subroutine  argument  IRD  is  the  index  of  the  standard  distance  which  is  within  one 
percent  of  the  reference  distance.  IRD  is  defined  in  this  subroutine  and  returned  to  the  MAIN 
routine.  For  the  present  standard  reference  distance  DIST  of  250  feet,  IRD  will  always  be  two. 

D1  is  a  constant  computed  in  the  MAIN  routine  and  defined  in  the  next  section,  n  is  a  power 
setting  related  index  (see  MAIN  routine). 

The  EA  data  statement  array  contains  the  excess  attenuation  data  for  frequencies  50  Hz  to 
10,000  Hz  and  distances  200  to  25,000  feet.  These  data  are  presented  in  Table  3  which  was 
taken  from  AAMRL-TR-84-0 17.(3)  jhis  table  also  contains  the  data  for  bands  10  through  16 
which  are  the  same  as  for  band  17. 
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Figure  6.  Flowchart  for  Subroutine  RSPLN. 


Extrapolated  SPL  Data 


The  first  segment  of  this  subroutine  (statement  label  135  loop)  computes  the  SPL 
spectrum  extrapolated  to  the  standard  profile  distance; 


SPUC,,,  =  FSPL,,„  F  D1  20Iog(5X, )  ^  ) 


1000 


1000 


dB 


where 

SPLXij 


FSPLij^l 

D1 


the  calculated  band  SPL  in  dB  at  the  profile  distance  from  the  source  for 
angle  index  I  and  frequency  index  J. 

reference  band  SPL  in  dB  at  distance  DIST  and  for  angle  index  I,  band 
index  J,  and  power  setting  index  L. 

defined  in  the  MAIN  routine 


FIMPR8 

DIST 

EAD 


SXk 

ATN8j 


ATNCj 


=  10  log  (FIMPR8)  +  20  log  (DIST) 

=  the  ratio  of  profile  and  reference  characteristic  impedance  (see  subroutine 

TESTN). 

=  reference  distance  in  feet  from  the  source  (250  feet). 

=  excess  atmospheric  attenuation  of  sound  in  dB  over  distance  SXj^  for 
fi-equency  band  index  J. 

the  standard  profile  distance  in  feet  (defined  in  the  MAIN). 

sound  absorption  coefficient  in  dB  per  1000  feet  for  frequency  index  J  and 
profile  output  temperature  and  relative  humidity. 

same  as  above  except  for  standard  day  conditions  (reference) 


Frequency  Weighted  Measures 

The  statement  label  250  loop  controls  the  computation  of  the  perceived  noise  level 
(PNLX),  A-weighted  overall  sound  level  (ALX)  and  tone  correction  (CXD  —  reference  distance 
only)  for  each  spectrum  (angle)  for  the  profile  distance  and  II^  power  setting.  These  data 
are  computed  only  when  requested  by  the  MEAS  program  flag.  Subroutine  CPNL  is  called  to 
compute  the  PNLX  data  for  the  spectrum.  If  PNLX  data  are  missing  beyond  the  second 
distance,  PNLX  is  extrapolated  linearly  from  the  last  two  good  points: 


=  (2)(PNLX, )  -  PNLX,_^_^„ 


PNdB 


The  A-weighted  overall  sound  level  for  the  I*  spectrum  are  computed  by  subroutine 
CAL.  These  ALX  data  are  the  final  A-weighted  profile  data.  The  tone  correction  data  are 
computed  for  the  spectrum  and  reference  distance  (ERD)  by  subroutine  CPTC. 


>in2le  Event  Profile  Data 


The  PNLTX  and  ALTX  profile  data  are  computed  by  adding  a  smoothed  tone  correction 
to  the  PNLX  and  ALX  data  (label  400  loop): 


P^LTK,.k.„  =  PNU:,_,j,+  {Cl){CXD, )  PNdB 


where 

CXDj  =  the  tone  correction  in  dB  for  the  ph  spectrum  at  the  reference  distance. 

Cl  =  1 .0  for  distances  200  to  3, 1 50  feet. 

=  (0.2)(18-K)  for  distances  4,000  to  8,000  feet  (K=14  to  17). 

=  0.0  for  distances  10,000  to  25,000  feet. 


SUBROUTINE  CPNL(FJ,I,II) 


This  subroutine  is  called  from  subroutine  CDIST  to  compute  the  perceived  noise  level 
(PNLX)  for  the  ph  spectrum  (angle),  distance  and  Iph  power  setting.  Function  FNOY  is 
used  to  compute  the  noy  data  for  each  SPL  data  point.  Subroutine  argument  FJ  is  0.15  for  this 
1/3  octave  band  data.  The  PNLX  algorithms  applied  in  this  subroutine  are  the  same  as  applied  in 

subroutine  CPNL  in  the  OMEGA  10  program;  thus  no  additional  documentation  will  be  given 
here. 


FUNCTION  FNOY(SPL,JJ) 

This  function  is  used  in  subroutine  CPNL.  It  is  the  same  as  function  FNOY  described  in 
the  OMEGA  10  program  documentation. 


SUBROUTINE  CPTC(PTC,I) 

This  subroutine  is  called  from  subroutine  CDIST  to  compute  the  tone  correction  (PTC) 
for  the  Fh  spectrum  (angle)  and  the  reference  distance.  This  subroutine  is  the  same  as  subroutine 
CPTC  described  in  the  OMEGA  10  program  documentation. 
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SUBROUTINE  CAL(I,  II) 


This  subroutine  is  called  from  subroutine  CDIST  to  compute  the  A-weighted  overall 
sound  level  (ALX)  for  the  1*^^  spectrum,  the  distance,  and  the  power  setting.  The  ALX 
algorithm  is  the  same  as  described  for  subroutine  CAL  in  the  OMEGA  10  documentation.  Many 
of  the  variable  names  are  different  in  this  subroutine;  however,  they  should  be  defined  in 
sufficient  detail  in  comment  statements  in  the  listing. 

SUBROUTINE  PPFDAT(JIJ2JI, LI, L2,IPR,IEDIT) 

This  subroutine  is  called  from  the  MAIN  routine  to  write  (IEDIT>-1)  the  PNLTX,  ALX, 
and/or  ALTX  single  event  profile  data  to  the  file  opened  on  Unit  2  and  to  print  (IPR>0)  the 
PNLX,  PNLTX,  ALX  and/or  ALTX  list  output  to  the  OUTPUT  file.  For  both  types  of  output, 
only  measures  requested  by  the  MEAS  flag  are  printed.  To  simplify  the  coding  and  reduce  the 
number  of  write  statements,  these  profile  data  are  stored  in  the  blank  common  array  SENX. 

The  subroutine  arguments  are  defined  as  follows: 

(1)  J1  and  J2  are  the  first  and  last  indices  of  the  SENX  profile  data  and  JI  is  the 
increment  of  Jl  to  J2. 

(2)  LI  and  L2  are  the  indices  of  the  normalized  data  used  to  compute  (interpolate)  the 
profile  data  for  the  requested  power  condition. 

(3)  IPR  is  the  program  list  print  flag. 

(4)  lEDIT  is  the  profile  flag  which  controls  the  quantity  of  data  printed  to  the  file  on 
Unit  2. 

Write  Profile  Datasets  to  Unit  2 

The  number  of  angles  of  profile  data  written  to  Unit  2  by  the  label  100  loop  depends  on 
the  value  of  the  EEDIT  program  flag.  All  19  angles  or  the  10  angles  selected  by  the  EDIT  routine 
are  written  for  BEDIT  equal  to  one  and  zero,  respectively.  A  complete  description  of  the  content 
and  format  of  the  10  angle  ground  run-up  profile  datasets  is  given  in  Appendix  H.  When  all 
angles  are  written,  the  final  10  angles  which  most  accurately  describe  each  noise  profile  must  be 
selected  from  these  19  angles.  The  plot  of  the  profile  data  for  the  reference  distance,  printed  by 
subroutine  PLT  on  output  Page  J,  is  included  to  aid  in  the  selection  of  these  10  angles. 
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Print  Sin2le  Event  Data 


The  remainder  of  this  subroutine  (label  400  loop)  prints  the  four  single  event  measures 
(PNLX,  PNLTX,  ALX,  and  ALTX)  on  Pages  D  through  G,  respectively.  These  data  are  printed 
from  the  SENX  array  using  the  Jl,  J2,  and  J1  subroutine  arguments  defined  above.  Subroutine 
HEADS  is  called  to  print  each  page  header  block.  All  single  event  data  less  than  zero  are 
blanked  out  in  the  printout.  The  content  and  format  of  these  pages  can  best  be  observed  by 
consulting  the  sample  problem  in  Appendix  D. 

SUBROUTINE  TITPG 

This  subroutine  is  called  from  the  MAIN  routine  to  print  the  cover  or  title  page  for  each 
test.  This  cover  page  provides  the  following  information: 

(1)  The  aircraft  name  and  code  for  which  the  data  were  measured. 

(2)  The  program  used  in  the  test  analysis. 

(3)  The  date  of  the  computer  run. 

(4)  The  types  of  data  computed  and  printed  for  the  test. 

The  content  and  format  of  this  title  page  can  best  be  understood  by  consulting  the  sample 
problem  in  Appendix  D.  With  this  sample  problem  title  page  as  a  guide,  there  should  be  no 
problem  in  following  the  coding. 

SUBROUTINE  PLT(IRD,JJI,JJ2JJI) 

This  subroutine  is  called  from  the  MAIN  routine  to  print  a  plot  of  PNLTX,  ALX,  and 
ALTX  versus  angle.  These  data  are  for  the  reference  distance  (IRD  =  250  feet)  from  the  profile 
datasets.  The  content  and  format  of  this  plot  can  best  be  understood  by  consulting  Page  J  in  the 
sample  problem  in  Appendix  D.  This  plot  simplifies  the  selection  or  checking  of  the  10  angles 
which  most  accurately  describe  the  profile  data. 

The  JJl,  JJ2  and  JH  subroutine  arguments  are  the  indices  and  the  index  increment  of 
these  profile  measures  in  array  SENX.  IRD  is  the  index  of  the  reference  distance  in  array  SENX. 

Method 

The  first  segment  of  this  subroutine  initializes  the  P  plot  array,  calls  subroutine  HEADS 
to  print  the  page  header  block,  and  prints  the  plot  symbol  identification  line,  the  top  border  line 
and  the  first  grid  Une  below  the  header  block.  The  label  25  loop  determines  the  maximum 
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PNLTX  value  which  is  then  used  to  set  up  the  abscissa  scale  values.  The  maximum  annotated 
scale  value  is  the  first  multiple  of  ten  greater  than  this  maximum  PNLTX.  The  minimum  value 
is  100  less  than  the  maximum.  The  actual  minimum  and  maximum  values  are  two  less  and  two 
greater  than  these  annotated  scale  values. 

The  label  200  loop  sets  up  and  prints  the  plot  for  each  of  the  19  angles.  The  ordinate 
annotation  and  title  and  grid  pattern  are  determined  from  the  angle  index  I.  The  label  120  loop 
scales  the  PNLTX,  ALX  and  ALTX  data  for  the  angle  and  sets  up  the  corresponding  symbol 
in  the  plot  array.  The  data  are  then  printed  and  the  plot  array  is  reinitialized  with  the  data  stored 
in  array  S AV.  There  should  be  sufficient  comments  in  the  coding  to  follow  the  setup  and 
printing  of  the  data. 

SUBROUTINE  RANK(IREQ,IERR) 

This  subroutine  is  called  from  the  MAIN  routine  to  determine  the  reference  data  (one  or 
two  datasets)  required  to  compute  the  profile  output  for  each  requested  output  power  setting 
(PSC  or  PSCF).  The  indices  of  these  reference  data  are  stored  in  array  EREQ  for  each  power 
setting.  Subroutine  argument  lERR  is  a  program  error  flag  which  is  returned  greater  than  zero 
when  errors  occurred  in  this  subroutine. 

Method’ 

Internal  reads  are  used  to  convert  the  input  (PS)  and  output  (PSC)  power  setting  data 
from  character  format  to  floating  point.  The  floating  point  data  are  required  for  ranking  and 
interpolating  the  data. 

Next  this  subroutine  ranks  all  normal  reference  file  power  setting  data  (PSIF)  and  all 
normal  requested  output  power  setting  data  (PSCF).  The  special  case  data  which  may  not  be 
interpolated  (Afterburner,  Wet  and  With  Jets)  are  not  ranked.  These  ranked  power  setting  data 
are  then  used  to  determine  which  reference  data  are  required  to  interpolate  profile  data  for  each 
(IC)  output  power  setting.  Indices  of  these  reference  data  are  defined  in  array  IREQ  for  each 
output  power  setting.  If  profile  data  for  the  same  reference  power  setting  are  required  to 
interpolate  profile  data  for  two  or  more  consecutive  output  power  settings,  the  index  of  this 
reference  power  setting  is  stored  in  consecutive  columns  in  the  same  row  in  array  IREQ. 

Profile  data  at  these  reference  power  settings  are  later  extrapolated  from  the  reference 
SPL  spectra  (see  MAIN  routine)  as  requested  in  array  IREQ;  thus,  only  those  profile  reference 
data  required  to  interpolate  data  at  the  profile  output  power  settings  are  extrapolated  by  the 
program.  Also,  to  conserve  storage,  reference  profile  data  storage  is  limited  to  the  two  power 
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settings  required  to  compute  each  output  profile  power  setting;  however,  IREQ  is  checked  to 
avoid  duplicate  reference  profile  computations. 

There  should  be  sufficient  comments  in  the  coding  to  follow  the  setup  of  the  IREQ  index 

array. 

SUBROUTINE  SUMRY(IREQ,IEDIT,FMXER) 

This  subroutine  is  called  from  the  MAIN  routine  to  print  the  OMEGA  1 1  summary  page 
which  lists  job  identification  parameters  as  well  as  a  summary  of  the  input  and  output  data.  This 
page  is  the  only  list  output  when  the  program  print  flag,  IPR,  is  less  than  one.  The  subroutine 
arguments  are  defined  as  follows: 

(1)  IREQ  is  an  array  containing  the  indices  of  the  reference  datasets  used  to  interpolate 
the  profile  data. 

(2)  EDIT  is  a  program  flag  which  controls  the  editing  of  the  profile  data. 

(3)  FMXER  is  the  maximum  error  permitted  in  the  EDIT  subroutine. 

If  the  sample  problem  summary  page  is  used  as  a  guide,  the  coding  in  this  subroutine 
should  be  easy  to  follow. 

SUBROUTINE  EDIT(IRDJ1J2,JI^CC,PSC,PSU, FMXER) 

This  subroutine  is  called  from  the  MAIN  routine  to  select  the  ten  angles  for  each  measure 
(PNLTX,  ALX  and  ALTX)  which  best  define  the  angle  versus  noise  level  data  for  that  measure 
at  the  reference  distance  (250  feet).  Actually  only  eight  angles  are  selected  by  this  routine,  since 
angles  0  and  180  degrees  are  always  included.  Reference  distance  profile  data  for  all  three 
measures  at  all  19  angles  are  used  to  compute  the  three  sets  of  ten  angles,  even  though  only  one 
or  two  of  the  measures  are  required  for  the  data  analysis  (see  MEAS  program  flag). 

The  subroutine  arguments  are  defined  as  follows: 

(1)  IRD  is  the  index  of  the  reference  distance. 

(2)  J1  and  J2  are  the  indices  of  the  PNLTX  and  ALTX  data  in  array  SENXD  and  JI  is 
the  index  increment  (ALX  data  is  in  index  J 1  +  JI). 

(3)  ACC  is  aircraft  code,  PSC  is  the  power  setting  for  which  the  data  are  being 
computed,  and  PSU  are  the  power  setting  units. 
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(4)  FMXER  is  the  maximum  error  permitted  in  this  EDIT  routine  (default 
FMXER=5.0). 


Method 

The  following  outline  describes  the  method  used  to  obtain  the  best  set  of  ten  angles: 

(1)  Compute  the  18  angle  to  angle  changes  in  PNdB  or  dBA  for  each  measure  (slopes) 
and  store  in  array  SL. 

(2)  Compute  the  17  changes  in  slope  (DSL)  for  each  measure  and  assign  to  angles  10, 
20, ...,  170  degrees.  Store  these  data  in  array  DSL. 

(3)  Rank  the  absolute  values  of  these  17  changes  in  slope  (DSL)  for  each  measure 
(label  50  loop). 

(4)  The  eight  angles  with  the  largest  rank  (i.e.,  largest  change  in  slope;  rank  10  to  17) 
plus  angles  0°  and  180°  will  be  the  initial  guess  for  the  desired  10  angles  for  each 
measure.  There  are  three  sets  of  10  angles,  one  for  each  measure.  The  rank  data  are 
stored  in  array  NR. 

(5)  Using  the  selected  10  angles  and  linear  interpolation,  compute  the  difference 
between  the  original  and  interpolated  data  for  each  angle;  i.e.,  compute  the  error. 
The  error  for  the  10  selected  angles  will  always  be  zero;  thus,  the  RMS  of  the  error 
is  computed  only  for  the  nine  remaining  angles: 

9 

These  error  data  are  computed  by  subroutine  ERR. 

(6)  Compute  the  error  and  RMS  error  described  in  step  5  above  for  each  measure  three 
times  using  the  three  sets  of  angles  (step  4)  for  each  measure.  For  example,  for 
PNLTX  compute  the  RMS  of  the  error  using  the  angles  selected  for  PNLTX,  ALX 
and  ALTX.  For  each  measure,  save  the  angles,  error  data  and  RMS  for  the  smallest 
RMS  error  (see  statement  label  150  and  2(X)  loops). 

(7)  Determine  the  angle  with  the  largest  error  using  the  error  data  in  step  6  above  (for 
one  measure  at  a  time).  If  this  error  is  greater  than  ERMAX  (ERMAX=1.49  dB), 
select  this  angle  as  one  of  the  eight  to  be  chosen  and,  one  at  a  time,  delete  each  of 
the  eight  angles  previously  selected  computing  the  error  and  RMS  error  for  each  of 
the  eight  sets.  Determine  the  best  of  the  eight  RMS  values  and,  if  better  than  the 
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RMS  in  step  6  above,  revise  the  angle  set  for  this  measure.  If  no  improvement  is 
found,  repeat  the  above  using  the  second  largest  error  greater  than  ERMAX.  This 
step  is  performed  by  subroutine  ITER  called  from  this  routine. 

(8)  Repeat  step  7  above  a  maximum  of  five  times  (in  subroutine  ITER)  or  until  all  error 
are  less  than  ERMAX  (maximum  of  five  passes  including  the  second,  third  etc., 
largest  error  passes). 

(9)  Repeat  steps  7  and  8  above  for  each  of  the  three  measures  (label  200  loop). 

(10)  At  this  point  we  have  revised  sets  of  angles  (compared  with  step  6  above)  for  a 
measure  only  if  one  or  more  errors  were  greater  than  ERMAX  for  that  measure.  If 
any  angle  sets  have  been  revised,  repeat  steps  6  through  9  above  using  the  new 
angle  sets  instead  of  those  from  step  4;  otherwise,  use  angle  sets  selected  in  step  6 
as  final  values  for  each  measure.  The  repeat  of  steps  6  through  9  are  computed  by 
the  statement  label  220  and  230  loops  which  also  use  subroutines  ERR  and  ITER. 

(11)  Return  the  angle  set  for  each  requested  measure  to  the  MAIN  routine  in  array  NR. 
The  angle  set  consists  of  those  angles  whose  rank  is  from  10  to  17  in  array  NR. 

Print  error  messages  if  the  interpolation  errors  using  these  final  selected  angles 
exceed  FMXER  (default  FMXER=5.0). 

SUBROUTINE  ERR(J,JA,NRA,RMSA,ERA,JM) 

This  subroutine  is  called  from  subroutines  EDIT  and  ITER  to  compute  the  error  data  as 
described  in  step  5  in  subroutine  EDIT.  The  subroutine  arguments  are: 

(1)  J  is  the  measure  index  of  the  original  data. 

(2)  JA  is  the  measure  index  in  the  rank  (NRA),  error  (ERA)  and  RMS  error  (RMS A) 
arrays.  JA  indicates  the  angle  set  being  used  to  interpolate  the  above  J  measure 
data. 

(3)  NRA,  ERA  and  RMSA  are  defined  in  (2)  above. 

(4)  JM  is  the  variable  dimension  of  arrays  NRA,  ERA  and  RMSA. 

Method 

Subroutine  FINTP  is  called  to  linearly  interpolate  the  PNdB  or  dBA  levels  of  the  nine 
angles  with  rank  less  than  ten.  The  error  data  are  then  computed  as  follows: 

ERAi  JA  =  SENXjj  -  DBCj 
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where 


RMSAj^ 


ERAj  =  the  interpolation  error  in  PNdB  or  dB  A  for  the  angle  and  measure  index 
JA. 

SENXj  j  =  the  computed  measure  data  in  PNdB  or  dB  A  for  angle  I  and  measure  J. 

DBCj  =  the  interpolated  measure  value  in  PNdB  or  dBA  for  angle  I. 

RMS  Aja  =  the  RMS  error  for  measure  index  JA. 

The  above  error  data  (ERAj  j/^)  are  zero  for  the  eight  angles  with  rank  greater  than  nine; 
thus  the  RMS  error  are  computed  by  dividing  the  sum  of  the  error  data  by  nine. 

SUBROUTINE  FINTP(AG,DBC,X,Y) 

This  is  a  linear  interpolation  subroutine  called  by  subroutine  ERR  to  interpolate  the 
measure  data  for  angles  with  rank  less  than  ten.  Arrays  X  and  Y  contain  the  angle  and 
corresponding  noise  level  data  (PNdB  or  dBA)  for  angles  0  and  180  degrees  plus  the  eight  angles 
with  rank  greater  than  nine.  Array  AG  contains  the  19  angle  values  in  ten  degree  increments 
from  0  to  180  degrees.  Array  DBG  contains  the  noise  level  computed  by  this  subroutine  for  the 
19  angles. 

Method 

For  angles  defined  in  array  X  (rank  greater  than  9),  DBCj  is  set  equal  to  the  given  noise 

level: 

DBCj  =  Yl  PNdB  or  dBA 

where 

DBCj  =  noise  level  for  the  angle. 

Yl  =  the  given  noise  level  corresponding  to  angle  Xl  which  is  the  same  as  angle 

AGj. 

For  angles  which  must  be  interpolated  the  linear  interpolation  function  is: 

^  Yf  —Y  \ 

DBCj  PNdBordBA 

V  X/,  -  At,  J 

where 
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DBCj 

=  the  interpolated  noise  level  for  the  angle. 

XI 

=  angle  AGj 

L 

=  the  index  of  the  first  angle  in  array  X  greater  than  XI. 

LI 

=  L-1. 

XL 

=  first  angle  greater  than  XI. 

Xli 

=  first  angle  less  than  XI. 

yl 

=  given  noise  level  corresponding  to  angle  Xl- 

Yu 

=  given  noise  level  corresponding  to  angle  X^j. 

SUBROUTINE  ITER(J,JMN,ICK,ERMAX,NRD) 


This  subroutine  is  called  from  subroutine  EDIT  to  perform  the  interations  described  in 
steps  7  and  8  in  subroutine  EDIT.  The  subroutine  arguments  are; 

(1)  J  is  the  index  of  the  original  measure  data  and  of  the  rank  (NRD)  and  RMS  error 
(RMS)  data. 

(2)  JMN  is  the  index  of  the  minimum  error  set  in  array  ER. 

(3)  ICK  is  an  error  flag  which  is  one  when  an  error  greater  than  ERMAX  was  found  in 
this  routine.  This  usually  results  in  a  change  in  rank  (angles  selected). 

(4)  ERMAX  is  a  dB  level  defined  in  subroutine  EDIT.  An  attempt  is  made  to  improve 
the  angle  selection  for  interpolation  errors  greater  than  or  equal  to  this  level. 

(5)  NRD  is  a  dummy  rank  array. 


SUBROUTINE  RDFILE(NFRUNI,NFRUN2*) 

This  subroutine  is  called  from  the  MAIN  routine  when  the  command  line  variable  ARGS 
in  the  MAIN  routine  is  NULL.  The  subroutine  arguments  NFRUNl  and  NFRUN2  are  two 
twenty-four  byte  character  variables  used  to  store  the  NOISEFILE  database  filenames  for 
military  runup  data.  When  the  command  line  is  NULL,  these  two  filenames  must  be  defined  in 
the  main  routine.  The  asterisk  in  the  argument  list  is  an  alternate  RETURN  label  which  is  used 
to  return  control  to  the  calling  routine  when  an  error  has  occurred. 

This  subroutine  performs  the  following  functions: 

(1)  Unit  1  is  opened  for  I/O  to  the  display; 
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(2)  The  program  queries  the  operator  for  the  OMEGA  1 1  setup  filename  which  must  be 
a  standard  DOS  filename  with  a  maximum  of  twelve  characters.  This  file  is  opened 
on  Unit  5  by  the  program.  The  ERRFLG  is  set  to  TRUE  if  the  operator  does  not 
enter  an  existing  filename. 

(3)  The  program  queries  the  operator  for  the  profile  dataset  output  filename  which  must 
be  a  standard  DOS  filename  with  a  maximum  of  twelve  characters.  This  file  is 
opened  on  Unit  2  by  program. 

(4)  The  NOISEFILE  database  files  with  names  stored  in  NFRUNl  and  NFRUN2  are 
checked  to  be  sure  they  exist.  If  they  don't  exist  the  ERRFLG  is  set  to  TRUE  and  the 
job  is  terminated.  These  files  are  not  opened  at  this  time. 

(5)  The  list  file  is  always  opened  on  Unit  6  with  filename  OMEGA  1  LOUT.  If  this  file 
already  exists,  the  operator  is  queried  for  a  new  filename  and  given  the  option  to 
terminate  the  job.  If  the  job  is  terminated,  the  ERRFLG  is  set  to  TRUE. 

SUBROUTINE  CLFILE(NFRUNI,NFRUN2^GS, *) 

This  subroutine  is  called  when  the  command  line  character  variable  ARGS  contains 
filename  information  which  must  be  decoded  into  the  I/O  filenames.  This  routine  separates  the 
command  line  string  which  may  be  up  to  127  character  long  into  the  I/O  filenames.  This  routine 
will  only  run  when  compiled  with  the  Lahey  compiler;  it  will  not  run  with  Microsoft  FORTRAN 
which  was  formerly  used  with  the  OMEGA  1 1  program. 

The  subroutine  arguments  NFRUNl  and  NFRUN2  are  two  twenty-four  byte  character 
variables  used  to  store  the  NOISEFILE  database  filenames  for  military  runup  data. 

The  command  line  must  contain  five  I/O  filenames  in  the  sequence  listed  here.  Each  file 
is  described  briefly  as  follows: 

(1)  The  setup  filename  is  decoded  and  opened  on  Unit  5. 

(2)  The  list  filename  is  decoded  and  opened  on  Unit  6. 

(3)  The  profile  dataset  filename  is  decoded  and  opened  on  Unit  3. 

(4)  The  first  NOISEFILE  database  filename  is  decoded  and  stored  in  variable 
NFRUNl. 

(5)  The  second  NOISEFILE  database  filename  is  decoded  and  stored  in  variable 
NFRUN2. 
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Each  filename  must  be  separated  by  a  space  on  the  command  line.  The  first  byte  on  the 
command  line  must  not  be  a  space.  If  an  error  occurs  when  opening  or  checking  a  file,  the 
ERRFLG  is  set  to  TRUE  and  the  job  is  terminated.  If  only  one  NOISEFILE  database  is 
available  and  if  all  data  to  be  analyzed  are  in  this  database,  then  this  database  filename  may  be 
entered  in  both  filename  fields  on  the  command  line.  The  two  database  files  are  never  open  at 
the  same  time. 

BLOCK  DATA 

The  BLOCK  DATA  routine  is  the  last  routine  in  this  program.  It  is  used  to  define 
numerous  variables  used  in  labeled  common  throughout  the  program.  See  the  program  listing  in 
Appendix  J  and  the  list  of  symbols  in  this  program  documentation  to  identify  the  variables 
defined  in  BLOCK  DATA. 
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APPENDIXA 

OMEGA  10  CODE  SHEETS  AND  SETUP  PROCEDURE 


This  appendix  contains  the  standard  procedure  for  setup  and  execution  of  the  OMEGA  10 
program  which  includes  the  OMEGA  10  program  code  sheets  and  a  detailed  description  of  each 
code  sheet  parameter  and  alphanumeric  data  field. 
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OMEGA  10  PROGRAM  CODE  SHEET 


I.  JOB  CONTROL  RECORD  (One  Per  JohV 

Col.  1  _  _  b  _  _  _b  _  _  DATE;  eg  10  JUN  93 

1 1  _____  _  DATN;  Data  in  DAMOYR  form,  e.g.,  100693 
1 8  _  IPR  [0  --  no  print] 

20 _  MEAS(l)  EPNL  [0  for  rPR=0;  1  for  IPR=1] 

22 _  MEAS(2)  ALM  [0  for  IPR=0;  1  for  1PR=I] 

24 _  MEAS(3)  SEL  [1  for  IPR=0;  1  for  IPR=1] 

26 _  IPU  [1  for  IPR=0;  0  for  IPR=1] 


II.  OUTPUT  PARAMETERS  FOR  EACH  AIRCRAFT  t2  or  3  Recordsl: 
Record  #1 


Col.  1  _  _  _  ACC 

4  _  _  _  ITEMP  [59”F] 

7  _  _  _  IRHUM  [70%] 

11  _  PV  [w] 

13  _  CRI  [0] 

15  _  _  _  _  _  DELN  [0.0] 

20  _  _  NP 

22 _  _  _  _  _  _  PSU  (left  justify) 

Record  #2  fProfile  Output  Power  Datal 

Col.  1  _  _  _  _  _PSC#1  Col.  6_  _  _  VX  Col.  9_  _  OPCR  Col.  11  _  _  OPCC  [81] 
13  _  _  _  _  _PSC#2  18  _  _  _  VX  21  _  _  OPCR  23  _  _  OPCC  [82] 

25  _  _  _  _  _  PSC  #3  30  _  _  _  VX  33  _  _  OPCR  35  _  _  OPCC  [83] 

37  _  _  _  _  _  PSC#4  42  _  _  _  VX  45  _  _  OPCR  47  _  _  OPCC  [84] 

49  _  _  _  _  _  PSC#5  54  _  _  _  VX  57  _  _  OPCR  59  _  _  OPCC  [85] 

61  _  _  _  _  _  PSC#6  66  _  _  _  VX  69  _  _  OPCR  71  _  _  OPCC  [86] 

Record  #3  (Profile  Output  Power  Data  Conti,;  required  only  when  NP>61 

Col.  1  _____  PSC  #  7  Col.  6  _  _  _  VX  Col.  9  _  _  OPCR  Col.  11  _  _  OPCC  [87] 

13  _  _  _  _  _  PSC#  8  18  _  _  _  VX  21  _  _  OPCR  23  _  _  OPCC  [88] 

25  _____  PSC  #  9  30  _  _  _  VX  33  _  _  OPCR  35  _  _  OPCC  [89] 

37  _  _  _  _  _  PSC  #10  42  _  _  _  VX  45  _  _  OPCR  47  _  _  OPCC  [90] 

49  _  _  _  _  _  PSC#11  54  _  _  _  VX  57  _  _  OPCR  59  _  _  OPCC  [91] 

61  _  _  _  _  _  PSC  #12  66  _  _  _  VX  69  _  _  OPCR  71  _  _  OPCC  [92] 


(Right  Justify  "PSC"  Data) 


(1 )  Repeat  section  II  for  each  aircraft  set  in  the  job. 

(2)  [  ]  —  Program  default  values  for  the  above  parameters. 

(3)  See  Standard  Procedure  for  Setup  and  Execution  of  the  OMEGA  10  Program  for  above  parameter 
definitions. 
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STANDARD  PROCEDURE  FOR  SETUP 
AND 

EXECUTION  OF  THE  OMEGA  10  PROGRAM 


1)  The  OMEGA  10  setup  deck  must  contain  the  following: 

(a)  One  job  control  record: 

(b)  One,  two  or  three  output  parameter  records  for  each  set  of  profile  data  for  each 
aircraft. 

The  parameters  required  for  each  record  are  listed  on  the  OMEGA  10  code  sheet  and 

described  in  detail  in  Steps  2  and  3  below. 

2)  Fill  in  the  OMEGA  10  code  sheet  job  control  record  parameters  where  default  conditions 

do  not  apply  (code  sheet  Item  I).  The  required  parameters  are  defined  as  follows: 

(a)  The  DATE  in  columns  1  to  9  will  be  printed  on  all  output  pages  and  in  the  first 
comment  line  in  each  profile  dataset  (9  alphanumeric  characters). 

(b)  The  date  in  the  day,  month,  year  form  (DATN)  in  columns  1 1  to  1 6  is  printed  as 
part  of  the  run  identification  on  all  the  plot  pages  (G,  J,  K,  N,  and  O).  It  is  not  used 
m  the  no-print  mode  and  thus  may  be  left  blank  on  the  code  sheet. 

(c)  The  value  of  IPR  (integer)  sets  the  print  (IPR=1)  or  no-print  (IPR=0  or  blank)  mode 
of  the  OMEGA  10  job.  In  the  no  print  mode  (default),  only  error  messages  and  the 
summary  page  are  printed  on  the  OUTPUT  file  (Unit  6)  and  only  the  measure(s) 
requested  on  the  code  sheet  is  computed.  In  the  print  mode,  all  seven  measures  are 
computed  and  all  output  are  printed  on  the  OUTPUT  file. 

(d)  The  MEAS(l),  MEAS(2)  and  MEAS(3)  parameters  (integers)  determine  which 
noise  measures  (EPNL,  ALM  and  SEE  respectively)  are  computed  by  the  program. 
In  the  print  mode  (IPR=1),  all  three  are  always  set  equal  to  one  by  the  program 
because  all  three  measures  are  required;  thus  they  may  be  left  blank  on  the  code 
sheet.  In  the  no-print  mode  (IPR=0  or  blank),  measure  data  will  be  computed  only 
for  the  one  or  more  MEAS  parameters  which  are  greater  than  zero  (default  is 
MEAS(3)=1  and  SEE  is  computed).  Note  that  the  PNE,  PNET,  SEET  and  AET 
noise  measures  are  always  computed  and  printed  in  the  print  mode;  however  they 
cannot  be  written  to  the  profile  dataset  file  opened  on  Unit  3. 

(e)  IPU  (integer)  controls  the  printing  of  the  profile  datasets  to  the  file  opened  on  Unit 
3.  For  IPU  equal  blank  or  zero  (default  for  IPR=1),  no  data  are  written  to  Unit  3 

For  IPU  equal  to  one  (default  for  IPR=0  or  blank),  all  requested  profile  data  are 
written  to  Unit  3. 
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Fill  in  the  OMEGA  10  code  sheet  output  parameters  for  each  aircraft  where  default 

parameters  do  not  apply  (code  sheet  Item  II).  One,  two,  or  three  records  are  required  for 

each  aircraft  set  depending  on  the  value  of  NP  defined  below.  The  required  parameters 

are  defined  as  follows: 

(a)  ACC  must  be  a  three  character  numeric  aircraft  code.  ACC  is  part  of  the  profile 
dataset  COMDECK  name  and  part  of  the  output  profile  identification  required  by 
the  NOISEMAP  program.  It  is  also  used  to  search  the  reference  file  for  the  required 
reference  data. 

(b)  H  EMP  and  IRUUM  (integers)  are  the  output  temperature  (°F)  and  relative 
humidity  (%),  respectively.  The  default  values  listed  on  the  code  sheet  are  for 
standard  day  conditions. 

(c)  PV  is  an  alphanumeric  profile  version  code  (usually  alpha)  which  is  printed  on  all 
output  and  is  the  second  last  character  in  the  profile  COMDECK  name.  One 
function  of  PV  would  be  to  make  COMDECK  names  unique  when  the  saune  aircraft 
data  are  run  for  different  weather  conditions  (NOISEMAP  does  not  read  the  PV 
code). 

(d)  CRI  is  an  alphanumeric  COMDECK  revision  identifier  (usually  integer)  which  is 
printed  on  the  summary  page  and  is  the  last  character  in  the  profile  COMDECK 
name.  CRI  is  designed  to  make  COMDECK  names  unique  when  several  revisions 
of  the  same  data  are  computed.  The  CRI  default  value  will  normally  apply  (not 
read  by  the  NOISEMAP  program). 

(e)  DELN  (or  "DELTA  N")  is  the  dB  level  to  be  added  to  all  reference  data  for  this 
aircraft  (floating  point  or  right  justified  integer).  It  is  printed  on  the  summary 
profile  data  pages  and  is  typically  used  to  adjust  data  for  multiple  engines. 

(f)  NP  is  the  number  of  power  settings  (PSCs)  to  be  computed  for  this  ACC.  If  NP  is 
zero  or  blank,  profile  data  are  computed  for  all  operation  power  codes  (for  ACC)  in 
the  reference  file  (max  of  15)  with  the  reference  file  airspeed  and  power  settings; 
also  the  output  operation  power  codes  (OPCC)  are  the  same  as  the  operation  power 
codes  in  the  reference  file  (that  is,  for  NP=0,  records  two  and  three  are  not  read  by 
the  program).  The  NP  entered  in  the  setup  file  must  be  an  integer  less  than  or  equal 
to  12. 

(g)  PSU  is  the  power  setting  unit  which  applies  to  the  PSC  data.  PSU  is  also  used  to 
select  the  reference  file  power  setting  data;  thus  it  must  exactly  match  (including 
blanks)  one  of  the  units  in  the  reference  file  or  the  ACC  data  will  be  terminated. 

PSU  must  be  left  justified  in  the  6  character  field  (alphanumeric).  For  the  NP=0 
option,  the  PSU  field  may  be  left  blank  in  which  case  the  data  for  die  first  power 
setting  will  be  used  to  compute  the  profile  data. 

(h)  The  PSCs  are  the  power  setting  values  for  which  profile  data  are  computed.  The 
PSC  data  may  be  listed  as  integer  or  floating  point  in  any  desired  sequence,  but  in 
either  case  must  be  right  justified  in  the  five  character  field.  Note  that  each 


reference  dataset  contains  the  source  power  setting  in  one,  two,  or  three  different 
units  in  the  third  comment  record.  The  PSU  parameter  described  in  item  3g  above 
specifies  which  of  these  power  settings  the  OMEGA  10  program  will  use  to 
inteipolate  the  profile  data.  The  PSCs  are  printed  in  the  third  comment  record  in 
each  profile  dataset,  on  all  profile  output  pages,  and  on  the  summary  page. 

(i)  The  VX's  are  the  airspeed  values  (right  justified  integers)  in  knots  to  which  the 
profile  data  are  adjusted.  They  are  printed  in  the  first  comment  line  in  each  profile 
dataset,  on  all  profile  data  output  pages,  and  on  the  final  summary  page. 

(j)  The  OPCR's  (numeric  with  leading  zero)  are  reference  file  operation  power  codes 
which  determine  the  operation  power  descriptions  of  the  profile  data  and  the 
reference  points  (SEL,  ALM,  or  EPNL  versus  power  setting)  from  which  the 
measure  data  for  the  PSC  power  settings  are  interpolated  or  extrapolated  (the  slope 
line  passes  through  this  reference  point).  The  COMDECK  names  of  the  reference 
file  datasets  for  this  reference  point  and  for  the  two  slope  points  are  printed  on  the 
summary  page. 

(k)  The  OPCCs  are  two  character  numeric  operation  power  codes  assigned  to  the 
output  profile  data.  The  OPCC  are  part  or  the  COMDECK  name  and  are  also  part 
of  the  profile  identification  used  by  the  NOISEMAP  program.  Default  values  are 

8 1  to  92.  For  NP  equal  to  zero,  the  OPCCs  are  set  equal  to  be  reference  file  OPCS. 

4)  The  following  comments  apply  to  the  PSC  and  OPCR  data: 

(a)  Each  reference  dataset  is  labeled  with  an  interpolation  type  code  of  I,  N  or  P.  This 
interpolation  type  code  which  is  included  as  the  last  character  in  the  COMDECK 
name  is  described  in  the  reference  dataset  format  description  in  Appendix  E.  For  all 
OPCR's  with  interpolation  type  code  N  the  PSC  power  setting  must  be  the  same  as 
the  reference  file  power  setting  or  the  PSC  data  will  be  deleted  from  the  run;  that  is, 
no  extrapolation  or  interpolation  is  permitted  for  these  OPCR's. 

(b)  When  the  PSC  power  setting  is  the  same  as  the  reference  file  power  setting 
designated  by  OPCR,  then  only  the  OPCR  reference  file  data  are  required  to 
compute  the  PSC  profile  data.  For  all  other  cases  where  PSC  is  not  the  same,  at 
least  one  additional  reference  file  dataset  is  required.  It  is  assumed  here  that  the 
individual  completing  the  code  sheet  is  familiar  with  the  interpolation/extrapolation 
requirements  and  limitations. 

5)  Repeat  Step  3  above  for  each  set  of  profile  data  for  each  aircraft  in  the  job.  If  two  or 
more  sets  of  profile  data  are  required  for  the  same  aircraft  because  of  additional  PSCs  or 
different  DELN  or  weather  data,  the  program  uses  the  reference  data  read  for  the  previous 
set  and  thus  avoids  searching  the  reference  file  for  the  same  data.  Since  NOISEMAP 
uses  an  ID  composed  of  the  ACC  and  OPCC  codes,  computer  jobs  run  for  NOISEMAP 
input  can  not  use  the  default  OPCC  for  multiple  sets  for  the  same  ACC;  default  OPCCs 
will  not  provide  unique  ID's  for  NOISEMAP.  For  jobs  unrelated  to  NOISEMAP,  PV  or 
CRI  codes  can  be  used  to  make  the  COMDECK  names  unique. 
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The  program  will  read  the  reference  data  from  the  file  opened  on  Unit  7.  These  data  must 
be  on  disk  in  the  format  described  in  Appendix  E.  The  program  rewinds  the  Unit  7  file 
before  searching  for  the  aircraft  (ACC)  data;  thus,  aircraft  sequence  is  not  important  in 
the  OMEGA  10  job  or  on  the  Unit  7  file  (all  data  for  each  aircraft  must  be  back  to  back  in 
the  Unit  7  file).  The  program  reads  and  stores  all  reference  data  for  all  operation  power 
codes  (MAX  of  15)  available  in  the  reference  file  for  aircraft  ACC.  If  the  reference  file 
contains  more  than  15  operation  power  codes  for  aircraft  ACC,  a  warning  message  is 
printed  and  only  the  first  15  are  read  and  stored  by  the  program. 

Execution  of  the  OMEGA  10  program  requries  approximately  150K  of  memory. 

After  execution  of  the  computer  job,  the  following  data  are  available: 

(a)  The  summary  pages  and  all  other  list  output  are  in  the  default  OMEGA108.OUT 
output  file  or  in  the  file  opened  for  list  output  on  Unit  6. 

(b)  The  profile  data  are  in  the  file  opened  on  Unit  3.  There  is  no  default  file  name. 

Data  for  aircraft  code  ACC  will  be  deleted  from  the  computer  job  and  an  error  message 
printed  when  one  or  more  of  the  following  problems  occur  (ERRFLG  is  set  to  TRUE): 

(a)  The  PSU  from  the  code  sheet  doesn't  match  any  of  the  available  power  setting  units 
for  ACC  in  the  reference  file. 

(b)  There  is  an  error  in  the  operation  power  code,  aircraft  code,  or  record  sequence  in 
one  or  more  reference  file  data  records. 

(e)  The  reference  file  minimum  slant  range  is  not  within  1%  of  a  standard  profile 
distance. 

(d)  No  reference  data  were  found  for  this  aircraft  in  the  reference  file. 

Data  for  operation  power  code  OPCC  will  be  deleted  from  the  aircraft  set  and  an  error 
message  printed  when  one  or  more  of  the  following  occur  (ERRFLG  is  set  to  TRUE): 

(a)  There  is  insufficient  data  for  extrapolation  or  interpolation  to  the  PSC  power 
setting. 

(b)  The  reference  file  and  the  PSC  power  settings  are  not  equal  as  required  for  this 
OPCR. 

(c)  The  reference  operation  power  code  (OPCR)  was  not  found  in  the  reference  file. 

Warning  messages  are  printed  when  one  or  more  of  the  following  occur  (ERRFLG  is  set 
to  TRUE): 

(a)  The  reference  file  contains  more  than  15  datasets;  only  the  first  15  are  read  for  this 
aircraft. 


(b)  The  reference  file  minimum  slant  range  is  not  equal  to  1000  feet  for  all  aircraft  and 
250  feet  for  some  helicopters  as  assumed  by  the  program  in  subroutines  DELTA6 
and  PPFDAT.  The  extrapolation  limit  check  is  always  computed  at  1000  feet  even 
for  these  helicopter  data.  This  will  probably  be  corrected/updated  when  the 
helicopter  procedure  is  finalized. 

12)  The  entire  computer  job  is  terminated  when  an  end  of  file  is  read  from  the  input  file.  This 
is  the  normal  job  termination. 
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APPENDIX  B 

OMEGA  11  CODE  SHEETS  AND  SETUP  PROCEDURE 


This  appendix  contains  the  standard  procedure  for  setup  and  execution  of  the  OMEGA  1 1 
program  which  includes  the  OMEGA  1 1  program  code  sheets  and  a  detailed  description  of  each 
code  sheet  parameter  and  alphanumeric  data  field. 
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OMEGA  11  PROGRAM  CODE  SHEET 


L  JOB  CQ>rrROL  RECORD  (One  Per  Job): 

Col.  DATE; eg  10 JUN 93 

12  _ IPR[0] 

13  _  DEDITIO) 

16 _  MEAS(l)  PNLT  [0] 

18 _  MEAS(2)  AL  [0] 

20  _  MEAS(3)  ALT  [0] 

21  _  _  _  _  _  FMXER  [5.0  dB) 


II.  OUTPUT  PARAMETERS  FOR  EACH  AIRCRAFT  (2  Records!: 
Recordll 


Col.  1  _  _  _  ACC 

8  _  _  __  IT8[59'‘F] 

11  _  _  _ _  P8[29.92  IN  HG] 

18  _  _  _  IH8[70%] 

22 _ PV[W] 

24 _ CR[0] 

26  _  _  _  _  _  DELN[0.0] 


Record  #2  (Profile  Output  Power  Data^ 

Col.  1  _  NP 

2  _____  _  PSU  (Left  Justify) 


PSC#1  Col. 

21  _____  PSC#2 
31  _  _  _  _  _  PSC#3 

41  _  _  _  _  _  PSC#4 

51  _  _  _  _  _  PSC#5 

61  _  _  _  _  _  PSC#6 


(Right  Justify  ”PSC”  and  weather  data) 

SECOND  AIRCRAFT 
Record  #1 

Col.  1  _  _  _  ACC 

8  _  _  _  IT8[59°F] 

11  _  _  _ _  P8[29.92  IN  HG) 

18  _  _  _  IH8[70%) 

22 _ PV[W] 

24  _  CR[0] 

22  _  _  _  _  _  DELN[0.0) 


17 

_  EFC 

Col.  19  _ 

_  OPCC[91] 

27 

IFC 

29  _ 

_  OPCC[92] 

37 

_  IFC 

39  _ 

_  OPCC[93) 

47 

_  IFC 

49  _ 

_  OPCC[94] 

57 

IFC 

59  _ 

_  OPCCI95) 

67 

_  IFC 

69  _ 

_  OPCC[96] 

Rccord_#2 

Col.  1  _  NP 

2  _  _  _ 
11  —  _  - 
21  _  _  . 
31  _  _  „ 
41  _  _  . 
51  _  _  . 
61  _  _  . 


PSU  (Left  Justify) 


PSC#1 

Col.  17 

_  IFC 

Col.  19  _ 

_  OPCC[911 

PSC  #2 

27 

_  IFC 

29  _ 

_  OPCC[92) 

PSC  #3 

37 

_  IFC 

39  _ 

_  OPCC[931 

PSC  #4 

47 

_  IFC 

49  _ 

_  OPCC[94] 

PSC  #5 

57 

_  IFC 

59  _ 

_  OPCC[951 

PSC  #5 

67 

IFC 

69  _ 

_  OPCC[96] 

(1)  Repeat  section  11,  records  1  and  2  for  each  aircraft  set  in  the  job. 

(2)  [  ]  —  Program  default  values  for  the  above  parameters. 

(3)  See  Standard  Procedure  for  Setup  and  Execution  of  the  OMEGA  1 1  Program  for  parameter  definitions. 
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STANDARD  PROCEDURE  FOR  SETUP 
AND 

EXECUTION  OF  THE  OMEGA  11  PROGRAM 


1)  The  OMEGA  1 1  setup  deck  must  contain  the  following: 

(a)  One  job  control  record; 

(b)  Two  output  parameter  records  for  each  set  of  profile  data  for  each  aircraft. 

The  parameters  required  for  each  record  are  listed  on  the  OMEGA  1 1  code  sheet  and 
described  in  detail  in  Steps  2  and  3  below. 

2)  Fill  in  the  OMEGA  1 1  code  sheet  job  control  parameters  where  default  conditions  do  not 
apply  (code  sheet  Item  I).  The  required  parameters  are  defined  as  follows: 

(a)  The  DATE  in  columns  1  to  9  will  be  printed  on  all  output  pages  and  in  the  first 
comment  line  in  each  profile  dataset  (9  alphanumeric  characters). 

(b)  The  value  of  IPR  determines  the  quantity  of  list  output  written  to  the  file  opened  on 
Unit  6.  If  IPR  is  zero  or  blank,  only  the  summary  page  and  program  error  messages 
are  printed  for  each  aircraft.  If  IPR  is  greater  than  zero,  all  list  output  are  printed 
(IPR  must  be  an  integer). 

(c)  The  value  of  BEDIT  determines  the  quantity  of  profile  dataset  output  written  to  the 
file  opened  on  Unit  2.  K  EDIT  is  zero  or  blank,  the  profile  datasets  for  each 
requested  measure  contain  data  for  the  10  angles  which  best  define  the  profile  data 
(using  linear  interpolation)  from  0°  to  180°.  If  EDIT  is  greater  than  zero,  profile 
data  for  all  19  angles  are  written  to  Unit  2;  for  EDIT  less  than  zero,  all  profile  data 
are  omitted  from  Unit  2  (EDIT  must  be  an  integer). 

(d)  The  MEAS(l),  MEAS(2)  and  MEAS(3)  parameters  determine  which  noise  measures 
(PNLT,  AL  and  ALT  respectively)  are  computed  by  the  program.  If  all  three  are 
zero  or  blank,  all  measures  are  computed;  otherwise,  only  measures  corresponding 
to  MEAS  greater  than  zero  are  computed  (all  MEAS  must  be  integer). 

(e)  The  value  of  FMXER  is  the  maximum  linear  interpolation  error  permitted  in  the 
angle  selection  routine.  For  errors  greater  than  FMXER,  error  messages  are  printed 
but  the  aircraft  analyses  are  not  terminated.  Typically,  the  largest  error  in  a 
relatively  smooth  profile  is  1 .0  to  1 .5  dB  or  less.  FMXER  which  applies  only  when 
EDIT  is  zero  or  blank  must  be  listed  as  a  floating  point  number  or  right  justified 
integer  (default  is  5.0  dB). 

3)  Fill  in  the  OMEGA  1 1  code  sheet  output  parameters  for  each  aircraft  where  default 
parameters  do  not  apply  (code  sheet  Item  II).  The  required  parameters  are  defined  as 
follows: 
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(a)  ACC  must  be  a  three  character  numeric  aircraft  code.  ACC  is  part  of  the  profile 
dataset  COMDECK  name  and  part  of  the  output  profile  identification  required  by 
the  NOISEMAP  program.  It  is  also  used  to  search  the  reference  file  for  the  required 
reference  data  (see  Step  4  below  regarding  multiple  sets  for  the  same  aircraft  code). 

(b)  ITS,  PS  and  IHS  are  the  output  temperature  (°F),  barometric  pressure  (in  Hg)  and 
relative  humidity  (%)  respectively.  The  default  values  listed  on  the  code  sheet  are 
for  standard  day  conditions.  ITS  and  IHS  must  be  integers  and  PS  must  be  floating 
point. 

(c)  PV  is  an  alphanumeric  profile  version  code  (usually  alpha)  which  is  printed  on  all 
list  output  and  is  the  second  last  character  in  the  profile  COMDECK  name.  One 
function  of  PV  would  be  to  make  COMDECK  names  unique  when  the  same  aircraft 
data  are  run  for  different  weather  conditions  (NOISEMAP  does  not  read  the  PV 
code). 

(d)  CRI  is  an  alphanumeric  COMDECK  revision  identifier  (usually  integer)  which  is 
printed  on  the  summary  page  and  is  the  last  character  in  the  profile  COMDECK 
name.  CRI  is  designed  to  make  COMDECK  names  unique  when  several  revisions 
of  the  same  data  are  stored  in  a  database.  The  CRI  default  value  will  normally 
apply  (not  read  by  the  NOISEMAP  program). 

(e)  DELN  (or  "DELTA  N")  is  the  dB  level  to  be  added  to  all  reference  data  for  this 
aircraft  (floating  point  or  right  justified  integer).  It  is  printed  on  the  summary  and 
reference  data  pages  and  is  typically  used  to  adjust  data  for  multiple  engines. 

(f)  NP  is  the  number  of  power  settings  (PSCs)  to  be  computed  for  this  ACC.  If  NP  is 
zero  or  blank,  profile  data  are  computed  for  all  operation  power  codes  (for  ACC)  in 
the  reference  file;  also  the  output  operation  power  codes  (OPCC)  are  the  same  as 
the  operation  power  codes  in  the  reference  file  (that  is,  for  NP=0,  record  2  columns 
2  to  70  are  not  read  by  the  program).  NP  must  be  an  integer  less  than  or  equal  to  6. 

(g)  PSU  is  the  power  setting  unit  which  applies  to  the  PSC  data.  PSU  must  exactly 
match  (including  blanks)  one  of  the  units  in  the  reference  file  or  the  ACC  data  will 
be  terminated.  PSU  must  be  left  justified  in  the  6  character  field  (alphanumeric). 
Note  that  the  program  checks  ail  available  power  setting  data  in  the  reference  file. 

(h)  The  PSC’s  are  the  power  setting  values  for  which  profile  data  are  computed.  Ail 
PSC  values  must  lie  within  the  extremes  available  in  the  reference  file.  For 
uniformity,  list  special  case  PSC  data  (IFC=1)  first  followed  by  the  normal  PSC  data 
(IFC=0);  each  in  sequence  from  low  to  high  (not  required  by  the  program).  The 
PSC  data  may  be  listed  as  integer  or  floating  point,  but  in  either  case  must  be  right 
justified  in  the  five  character  field.  Note  that  each  reference  dataset  contains  the 
source  power  setting  in  one,  two,  or  three  different  units  in  the  third  comment 
record.  Any  available  power  setting  (value  and  units)  may  be  used  to  interpolate 
the  profile  data. 
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(i)  IFC  is  a  program  flag  used  to  separate  the  AFTERBURNER,  WET,  WITH  JETS 
and  helicopter  special  case  data  from  the  normal  power  data.  IFC  must  be  zero  or 
blank  for  normal  data  and  one  (1)  for  special  case  data.  For  special  case  data 
(IFC=1),  no  interpolation  is  permitted  and  thus  PSC  must  be  exactly  the  same  as  the 
power  setting  in  the  normalized  file.  There  are  currently  14  special  case  power 
settings  which  are  identified  by  their  operation  power  codes.  The  14  speical  case 
codes  are:  01, 02, 03, 10, 35,  38, 42, 49,  60,  61,  62, 63, 64  and  65. 

(j)  The  OPCC  are  two  character  numeric  operation  power  codes  assigned  to  the  output 
profile  data.  The  OPCC  are  part  of  the  COMDECK  name  and  are  also  part  of  the 
profile  identification  used  by  the  NOISEMAP  program.  Default  values  are  91  to 
96. 

4)  Repeat  Step  3  above  for  each  set  of  profile  data  for  each  aircraft  in  the  job.  If  two  or 
more  sets  of  profile  data  are  required  for  the  same  aircraft  because  of  additional  PSC's  or 
different  DELN  or  weather  data,  ACC  should  be  set  equal  to  "***";  however,  all  other 
parameter  must  be  defined  or  default  conditions  apply.  When  ACC=***,  the  program 
uses  the  reference  data  read  for  the  previous  set  and  thus  avoids  searching  the  reference 
file  for  the  same  data.  Since  NOISEMAP  uses  an  ID  composed  of  the  ACC  and  OPCC 
codes,  computer  jobs  run  for  NOISEMAP  input  can  not  use  the  default  OPCC  for 
multiple  sets  for  the  same  ACC;  default  OPCC's  will  not  provide  unique  ID's  for 
NOISEMAP.  For  jobs  unrelated  to  NOISEMAP,  PV  or  CRI  codes  can  be  used  to  make 
the  COMDECK  names  unique. 

5)  The  program  will  read  the  reference  data  from  the  file  opened  on  Unit  7.  These  data  must 
be  on  disk  in  the  format  described  in  Appendix  G.  The  program  rewinds  the  Unit  7  file 
before  searching  for  the  aircraft  (ACC)  data;  thus,  aircraft  sequence  is  not  important  in 
the  OMEGA  1 1  job  or  in  the  dataset  file.  However,  all  data  for  each  aircraft  must  be 
back  to  back  in  the  reference  dataset.  The  program  reads  and  stores  all  reference  data  for 
all  operation  power  codes  (MAX  of  6)  available  on  the  reference  file  for  aircraft  ACC.  If 
this  reference  file  contains  more  than  six  operation  power  codes  for  aircraft  ACC,  all 
computations  for  that  aircraft  are  terminated  (one  exception:  if  all  IFC=0,  special  case 
data  are  not  stored).  The  first  Test  number  and  Noise  Source/Subject  description  read 
from  the  reference  file  (Unit  7)  are  the  Test  and  Noise  Source/Subject  description  used 
for  this  set  of  profile  data  for  this  aircraft.  Run  number  is  assigned  from  01  to  06  in  the 
sequence  of  the  PSC  output. 

6)  After  execution  of  the  computer  job,  the  following  data  are  available: 

(a)  The  summary  pages  and  all  other  computer  printouts  are  on  file  OMEGAl  l.OUT  or 
the  file  opened  on  Unit  6  as  the  list  output  file. 
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(b)  The  profile  data  are  in  the  file  opened  on  Unit  2.  There  is  no  default  profile  dataset 
file. 

7)  Execution  of  the  OMEGA  1 1  program  requires  approximately  150K  of  memory. 

8)  Data  for  aircraft  code  ACC  will  be  deleted  from  the  computer  job  and  an  error  message 

printed  when  one  or  more  of  the  following  problems  occur  (ERRFLG  is  set  to  TRUE): 

(a)  No  reference  data  are  found  for  ACC  in  the  reference  file  opened  on  Unit  7. 

(b)  The  Unit  7  file  contains  data  for  more  than  the  maximum  (6)  number  of  operation 
power  codes  for  ACC. 

(c)  Error  in  reference  dataset  input  caused  by  missing  or  extra  data  records  (angle  ID 
errors). 

(d)  The  PSU  from  the  code  sheet  doesn't  match  any  of  the  power  setting  units  for  ACC 
in  the  reference  file. 

(e)  The  requested  special  case  power  setting  (PSC  for  IFC=1)  was  not  found  in  the 
reference  file. 

(f)  One  or  more  requested  power  settings  (PSCs)  are  outside  the  range  of  the  power 
setting  data  available  in  the  reference  file. 

9)  The  entire  computer  job  will  be  terminated  when  either  of  the  following  occur  (ERRFLG 

is  set  to  TRUE): 

(a)  The  number  of  power  settings  requested  (NP)  is  greater  than  the  maximum  number 
permitted  (NPM=6).  An  error  message  will  be  printed. 

(b)  An  end  of  file  is  read  from  this  input  file.  This  is  the  normal  job  termination. 
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APPENDICES  C  THROUGH  K  MAY  BE  OBTAINED  BY  REQUEST  TO 


AL/OEBN 

2610  SEVENTH  ST. 

WRIGHT-PATTERSON  AFB  OH  45433-7901 
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